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Abstract 

A new gauge theory of gravity is presented. The theory is constructed 
in a flat background spacetime and employs gauge fields to ensure that 
all relations between physical quantities are independent of the positions 
and orientations of the matter fields. In this manner all properties of 
the background spacetime are removed from physics, and what remains 
are a set of 'intrinsic' relations between physical fields. For a wide range 
of phenomena, including all present experimental tests, the theory repro- 
duces the predictions of general relativity. Differences do emerge, however, 
through the first-order nature of the equations and the global properties of 
the gauge fields, and through the relationship with quantum theory. The 
properties of the gravitational gauge fields are derived from both classi- 
cal and quantum viewpoints. Field equations are then derived from an 
action principle, and consistency with the minimal coupling procedure se- 
lects an action that is unique up to the possible inclusion of a cosmological 
constant. This in turn singles out a unique form of spin-torsion interac- 
tion. A new method for solving the field equations is outlined and applied 
to the case of a time-dependent, spherically-symmetric perfect fluid. A 
gauge is found which reduces the physics to a set of essentially Newtonian 
equations. These equations are then applied to the study of cosmology, 
and to the formation and properties of black holes. Insistence on find- 
ing global solutions, together with the first-order nature of the equations, 
leads to a new understanding of the role played by time reversal. This al- 
ters the physical picture of the properties of a horizon around a black hole. 
The existence of global solutions enables one to discuss the properties of 
field lines inside the horizon due to a point charge held outside it. The 
Dirac equation is studied in a black hole background and provides a quick 
derivation of the Hawking temperature. Some applications to cosmology 
are also discussed, and a study of the Dirac equation in a cosmological 
background reveals that the only models consistent with homogeneity are 
spatially flat. It is emphasised throughout that the description of gravity 
in terms of gauge fields, rather than spacetime geometry, leads to many 
simple and powerful physical insights. The language of 'geometric alge- 
bra' best expresses the physical and mathematical content of the theory 
and is employed throughout. Methods for translating the equations into 
other languages (tensor and spinor calculus) are given in appendices. 
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1 Introduction 



In modern theoretical physics particle interactions are described by gauge theo- 
ries. These theories are constructed by demanding that symmetries in the laws 
of physics should be local, rather than global, in character. The clearest ex- 
positions of this principle are contained in quantum theory, where one initially 
constructs a Lagrangian containing a global symmetry. In order to promote this 
to a local symmetry, the derivatives appearing in the Lagrangian arc modified so 
that they are unchanged in form by local transformations. This is achieved by 
the introduction of fields with certain transformation properties ('gauge fields'), 
and these fields are then responsible for inter-particle forces. The manner in 
which the gauge fields couple to matter is determined by the 'minimal coupling' 
procedure, in which partial (or directional) derivatives are replaced by covariant 
derivatives. This is the general framework that has been applied so successfully 
in the construction of the standard model of particle physics, which accounts 
for the strong, weak and electromagnetic forces. 

But what of gravity: can general relativity be formulated as a gauge theory? 
This question has troubled physicists for many years [1, 2, 3]. The first work 
that recovered features of general relativity from a gauging argument was due to 
Kibble [2], who elaborated on an earlier, unsuccessful attempt by Utiyama [1]. 
Kibble used the 10-component Poincare group of passive infinitesimal coordi- 
nate transformations (consisting of four translations and six rotations) as the 
global symmetry group. By gauging this group and constructing a suitable 
Lagrangian density for the gauge fields, Kibble arrived at a set of gravitational 
field equations — though not the Einstein equations. In fact, Kibble arrived 
at a slightly more general theory, known as a 'spin-torsion' theory. The neces- 
sary modifications to Einstein's theory to include torsion were first suggested 
by Cartan [4], who identified torsion as a possible physical field. The connec- 
tion between quantum spin and torsion was made later [2, 5, 6], once it had 
become clear that the stress-energy tensor for a massive fermion field must be 
asymmetric [7, 8]. Spin-torsion theories are sometimes referred to as Einstein- 
Cartan-Kibble-Sciama theories. Kibble's use of passive transformations was 
criticised by Hehl et al. [9], who reproduced Kibble's derivation from the stand- 
point of active transformations of the matter fields. Hehl et al. also arrived at a 
spin-torsion theory, and it is now generally accepted that torsion is an inevitable 
feature of a gauge theory based on the Poincare group. 

The work of Hehl et al. [9] raises a further issue. In their gauge theory deriva- 
tion Hehl et al. arc clear that ' coordinates and frames are regarded as fixed once 
and for all, while the matter fields are replaced by fields that have been rotated 
or translated!. It follows that the derivation can only affect the properties of the 
matter fields, and not the properties of spacetime itself. Yet, once the gauge 
fields have been introduced, the authors identify these fields as determining the 
curvature and torsion of a Riemann-Cartan spacetime. This is possible only if it 
is assumed from the outset that one is working in a Riemann-Cartan spacetime, 
and not in flat Minkowski spacetime. But the idea that spacetime is curved is 
one of the cornerstone principles of general relativity. That this feature must 
be introduced a priori, and is not derivable from the gauge theory argument, is 
highly undesirable it shows that the principle of local gauge invariance must 
be supplemented with further assumptions before general relativity is recovered. 
The conclusions are clear: classical general relativity must be modified by the 
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introduction of a spin-torsion interaction if it is to be viewed as a gauge theory, 
and the gauge principle alone fails to provide a conceptual framework for general 
relativity as a theory of gravity. 

In this paper we propose an alternative theory of gravity that is derived from 
gauge principles alone. These gauge fields are functions of position in a single 
Minkowski vector space. But here we immediately hit a profound difficulty. 
Paramctcrising points with vectors implies a notion of a Newtonian 'absolute 
space' (or spacetime) and one of the aims of general relativity was to banish 
this idea. So can we possibly retain the idea of representing points with vectors 
without introducing a notion of absolute space? The answer to this is yes 
— we must construct a theory in which points are parameterised by vectors, 
but the physical relations between fields are independent of where the fields 
are placed in this vector space. We must therefore be free to move the fields 
around the vector space in an arbitrary manner, without in any way affecting the 
physical predictions. In this way our abstract Minkowski vector space will play 
an entirely passive role in physics, and what will remain are a set of 'intrinsic' 
relations between spacetime fields at the same point. Yet, once we have chosen 
a particular paramcterisation of points with vectors, we will be free to exploit 
the vector space structure to the full, secure in the knowledge that any physical 
prediction arrived at is ultimately independent of the paramcterisation. 

The theory we aim to construct is therefore one that is invariant under arbi- 
trary field displacements. It is here that we make contact with gauge theories, 
because the necessary modification to the directional derivatives requires the 
introduction of a gauge field. But the field required is not of the type usually 
obtained when constructing gauge theories based on Lie-group symmetries. The 
gauge field coupling is of an altogether different, though very natural, character. 
However, this does not alter the fact that the theory constructed here is a gauge 
theory in the broader sense of being invariant under a group of transformations. 
The treatment presented here is very different from that of Kibble [2] and Hchl 
et al. [9]. These authors only considered infinitesimal translations, whereas we 
are able to treat arbitrary finite field displacements. This is essential to our aim 
of constructing a theory that is independent of the means by which the positions 
of fields are parameterised by vectors. 

Once we have introduced the required 'position-gauge' field, a further space- 
time symmetry remains. Spacetime fields are not simply scalars, but also consist 
of vectors and tensors. Suppose that two spacetime vector fields are equated 
at some position. If both fields are then rotated at a point, the same intrinsic 
physical relation is obtained. We therefore expect that all physical relations 
should be invariant under local rotations of the matter fields, as well as dis- 
placements. This is necessary if we are to achieve complete freedom from the 
properties of the underlying vector space — we cannot think of the vectors rep- 
resenting physical quantities as having direction defined relative to some fixed 
vectors in Minkowski spacetime, but are only permitted to consider relations 
between matter fields. Achieving invariancc under local rotations introduces a 
further gauge field, though now we are in the familiar territory of Yang-Mills 
type interactions (albeit employing a non-compact Lie group). 

There are many ways in which the gauge theory presented here offers both 
real and potential advantages over traditional general relativity. As our theory 
is a genuine gauge theory, the status of physical predictions is always unam- 
biguous — any physical prediction must be extracted from the theory in a 
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gauge-invariant manner. Furthermore, our approach is much closer to the con- 
ventional theories of particle physics, which should ease the path to a quantum 
theory. One further point is that discarding all notions of a curved spacetime 
makes the theory conceptually much simpler than general relativity. For exam- 
ple, there is no need to deal with topics such as differentiable manifolds, tangent 
spaces or fibre bundles [10]. 

The theory developed here is presented in the language of 'geometric alge- 
bra' [11, 12]. Any physical theory can be formulated in a number of different 
mathematical languages, but physicists usually settle on a language which they 
feel represents the 'optimal' choice. For quantum field theory this has become 
the language of abstract operator commutation relations, and for general rela- 
tivity it is Ricmannian geometry. For our gauge theory of gravity there seems 
little doubt that geometric algebra is the optimal language available in which 
to formulate the theory. Indeed, it was partly the desire to apply this language 
to gravitation theory that led to the development of the present theory. (This 
should not be taken to imply that geometric algebra cannot be applied to stan- 
dard general relativity — it certainly can [11, 13, 14, 15]. It has also been 
used to elaborate on Utiyama's approach [14].) To us, the use of geometric 
algebra is as central to the theory of gravity presented here as tensor calculus 
and Ricmannian geometry were to Einstein's development of general relativity. 
It is the language that most clearly exposes the structure of the theory The 
equations take their simplest form when expressed in geometric algebra, and 
all reference to coordinates and frames is removed, achieving a clean separation 
between physical effects and coordinate artefacts. Furthermore, the geometric 
algebra development of the theory is entirely self-contained. All problems can 
be treated without ever having to introduce concepts from other languages, such 
as differential forms or the Newman-Penrose formalism. 

We realise, however, that the use of an unfamiliar language may deter some 
readers from exploring the main physical content of our theory — which is of 
course independent of the language chosen to express it. We have therefore 
endeavoured to keep the mathematical content of the main text to a minimum 
level, and have included appendices describing methods for translating our equa- 
tions into the more familiar languages of tensor and spinor calculus. In addition, 
many of the final equations required for applications are simple scalar equations. 
The role of geometric algebra is simply to provide the most efficient and trans- 
parent derivation of these equations. It is our hope that physicists will find 
geometric algebra a simpler and more natural language than that of differential 
geometry and tensor calculus. 

This paper starts with an introduction to geometric algebra and its space- 
time version - the spacetime algebra. We then turn to the gauging arguments 
outlined above and find mathematical expressions of the underlying principles. 
This leads to the introduction of two gauge fields. At this point the discussion 
is made concrete by turning to the Dirac action integral. The Dirac action is 
formulated in such a way that internal phase rotations and spacetime rotations 
take equivalent forms. Gauge fields are then minimally coupled to the Dirac 
field to enforce invariance under local displacements and both spacetime and 
phase rotations. We then turn to the construction of a Lagrangian density 
for the gravitational gauge fields. This leads to a surprising conclusion. The 
demand that the gravitational action be consistent with the derivation of the 
minimally-coupled Dirac equation restricts us to a single action integral. The 
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only freedom that remains is the possible inclusion of a cosmological constant, 
which cannot be ruled out on theoretical grounds alone. The result of this work 
is a set of field equations that are completely independent of how we choose to 
label the positions of fields with a vector x. The resulting theory is conceptu- 
ally simple and easier to calculate with than the metric-based theory of general 
relativity. We call this theory ' gauge theory gravity' (GTG). Having derived 
the field equations, we turn to a discussion of measurements, the equivalence 
principle and the Newtonian limit in GTG. We end Part I with a discussion of 
symmetries, invariants and conservation laws. 

In Part II we turn to applications, concentrating mainly on time-dependent 
spherically-symmetric systems. We start by studying perfect fluids and derive 
a simple set of first-order equations that describe a wide range of physical phe- 
nomena. The method of derivation of these equations is new and offers many 
advantages over conventional techniques. The equations are then studied in 
the contexts of black holes, collapsing matter and cosmology. We show how 
a gauge can be chosen that affords a clear, global picture of the properties of 
these systems. Indeed, in many cases one can apply simple, almost Newtonian, 
reasoning to understand the physics. For some of these applications the pre- 
dictions of GTG and general relativity arc identical, and these cases include 
all present experimental tests of general relativity However, on matters such 
as the role of horizons and topology, the two theories differ. For example, we 
show that the black-hole solutions admitted in GTG fall into two distinct time- 
asymmetric gauge sectors, and that one of these is picked out uniquely by the 
formation process. This is quite different to general relativity, which admits eter- 
nal time-reverse symmetric solutions. In discussing differences between GTG 
and general relativity, it is not always clear what the correct general relativistic 
viewpoint is. We should therefore be explicit in stating that what we intend 
when we talk about general relativity is the full, modern formulation of the 
subject as expounded by, for example, Hawking & Ellis [16] and D'Inverno [17]. 
This includes ideas such as worm-holes, exotic topologies and distinct 'universes' 
connected by black holes [18, 19]. 

After studying some solutions for the gravitational fields we turn to the 
properties of electromagnetic and Dirac fields in gravitational backgrounds. For 
example, we give field configurations for a charge held at rest outside a black 
hole. We show how these field lines extend smoothly across the horizon, and 
that the origin behaves as a polarisation charge. This solution demonstrates 
how the global properties of the gravitational fields are relevant to physics out- 
side the horizon, a fact that is supported by an analysis of the Dirac equation in 
a black-hole background. This analysis also provides a quick, though physically 
questionable, derivation of a particle production rate described by a Fermi-Dirac 
distribution with the correct Hawking temperature. We end with a discussion of 
the implications of gauge-theory gravity for cosmology. A study of the Maxwell 
and Dirac equations in a cosmological background reveals a number of surpris- 
ing features. In particular, it is shown that a non-spatially-flat universe does not 
appear homogeneous to Dirac fields — fermionic matter would be able to detect 
the 'centre' of the universe if k ^ 0. Thus the only homogeneous cosmological 
models consistent with GTG are those that are spatially flat. (This does not 
rule out spatially-flat universes with a non-zero cosmological constant.) A con- 
cluding section summarises the philosophy behind our approach, and outlines 
some future areas of research. 
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2 An outline of geometric algebra 



There are many reasons for preferring geometric algebra to other languages em- 
ployed in mathematical physics. It is the most powerful and efficient language 
for handling rotations and boosts; it generalises the role of complex numbers 
in two dimensions, and quaternions in three dimensions, to a scheme that effi- 
ciently handles rotations in arbitrary dimensions. It also exploits the advantages 
of labelling points with vectors more fully than either tensor calculus or differ- 
ential forms, both of which were designed with a view to applications in the 
intrinsic geometry of curved spaces. In addition, geometric algebra affords an 
entirely real formulation of the Dirac equation [20, 21], eliminating the need 
for complex numbers. The advantage of the real formulation is that internal 
phase rotations and spacetime rotations are handled in an identical manner in 
a single unifying framework. A wide class of physical theories have now been 
successfully formulated in terms of geometric algebra. These include classical 
mechanics [22, 23, 24], relativistic dynamics [25], Dirac theory [20, 21, 26, 27], 
electromagnetism and electrodynamics [12, 26, 28], as well as a number of other 
areas of modern mathematical physics [29, 30, 31, 32, 33]. In every case, geomet- 
ric algebra has offered demonstrable advantages over other techniques and has 
provided novel insights and unifications between disparate branches of physics 
and mathematics. 

This section is intended to give only a brief introduction to the ideas and 
applications of geometric algebra. A fuller introduction, including a number of 
results relevant to this paper, is set out in the series of papers [12, 21, 26, 31] 
written by the present authors. Elsewhere, the books by Doran & Lasenby [34], 
Hestenes [14, 22] and Hestenes & Sobczyk [11] cover the subject in detail. A 
number of other helpful introductory articles can be found, including those by 
Hestenes [35, 36], Void [24, 28], and Doran & Lasenby [37]. The conference 
proceedings [38, 39, 40] also contain some interesting and useful papers. 

Geometric algebra arose from Clifford's attempts to generalise Hamilton's 
quaternion algebra into a language for vectors in arbitrary dimensions [41]. 
Clifford discovered that both complex numbers and quaternions are special cases 
of an algebraic framework in which vectors are equipped with a single associative 
product that is distributive over addition 4 . With vectors represented by lower- 
case Roman letters (a, 6), Clifford's 'geometric product' is written simply as 
ab. A key feature of the geometric product is that the square of any vector is a 
scalar. Now, rearranging the expansion 

(a + b) 2 = (a + b){a + b) = a 2 + {ab + ba) + b 2 (2.1) 

to give 

ab + ba = [a + bf -a 2 -b 2 , (2.2) 

where the right-hand side of (2.2) is a sum of squares and by assumption a 
scalar, we see that the symmetric part of the geometric product of two vectors 
is also a scalar. We write this 'inner' or 'dot' product between vectors as 

a-b= \{ab + ba). (2.3) 

4 The same generalisation was also found by Grassmann [42], independently and somewhat 
before Clifford's work. This is one of many reasons for preferring Clifford's name {'geometric 
algebra') over the more usual 'Clifford algebra'. 
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The remaining antisymmetric part of the the geometric product represents the 
directed area swept out by displacing a along b. This is the 'outer' or 'exterior' 
product introduced by Grassmann [42] and familiar to all who have studied 
the language of differential forms. The outer product of two vectors is called a 
bivector and is written with a wedge: 

aAb= \{ab-ba). (2.4) 

On combining (2.3) and (2.4) we find that the geometric product has been 
decomposed into the sum of a scalar and a bivector part, 

ab = a-b + af\b. (2.5) 

The innovative feature of Clifford's product (2.5) lies in its mixing of two dif- 
ferent types of object: scalars and bivectors. This is not problematic, because 
the addition implied by (2.5) is precisely that which is used when a real number 
is added to an imaginary number to form a complex number. But why might 
we want to add these two geometrically distinct objects? The answer emerges 
from considering reflections and rotations. Suppose that the vector a is reflected 
in the (hyper)plane perpendicular to the unit vector n. The result is the new 
vector 

a — 2(a-n)n = a — (an + na)n = —nan. (2-6) 

The utility of the geometric algebra form of the resultant vector, —nan, becomes 
clear when a second reflection is performed. If this second reflection is in the 
hypcrplane perpendicular to the unit vector to, then the combined effect is 

a i — ► mnanm. (2-7) 

But the combined effect of two reflections is a rotation so, defining the geometric 
product mn as the scalar-plus-bivector quantity f2, we see that rotations are 
represented by 

a i-> RaR. (2.8) 

Here the quantity R = nm is called the 'reverse' of R and is obtained by 
reversing the order of all geometric products between vectors: 

(o6...c)~ =c...ba. (2.9) 

The object R is called a rotor. Rotors can be written as an even (geometric) 
product of unit vectors, and satisfy the relation RR = 1. The representation 
of rotations in the form (2.8) has many advantages over tensor techniques. By 
defining cos6* = m ■ n we can write 

R = mn = exp(^^r9/2), (2.10) 
\\mAn\ I 

which relates the rotor R directly to the plane in which the rotation takes place. 
Equation (2.10) generalises to arbitrary dimensions the representation of planar 
rotations afforded by complex numbers. This generalisation provides a good 
example of how the full geometric product, and the implied sum of objects 
of different types, can enter geometry at a very basic level. The fact that 
equation (2.10) encapsulates a simple geometric relation should also dispel the 
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notion that Clifford algebras are somehow intrinsically 'quantum' in origin. The 
derivation of (2.8) has assumed nothing about the signature of the space being 
employed, so that the formula applies equally to boosts as well as rotations. 
The two-sided formula for a rotation (2.8) will also turn out to be compatible 
with the manner in which observables are constructed from Dirac spinors, and 
this is important for the gauge theory of rotations of the Dirac equation which 
follows. 

Forming further geometric products of vectors produces the entire geometric 
algebra. General elements are called 'multivectors' and these decompose into 
sums of elements of different grades (scalars are grade zero, vectors grade one, 
bivectors grade two etc.). Multivectors in which all elements have the same grade 
are termed homogeneous and are usually written as A r to show that A contains 
only grade-r components. Multivectors inherit an associative product, and the 
geometric product of a grade-r multivector A r with a grade-s multivector B s 
decomposes into 

A r B s = (AB) r+s + (AB) r+s _ 2 ... + (AB) , r _ s | , (2.11) 

where the symbol (M) r denotes the projection onto the grade-r component of 
M. The projection onto the grade-0 (scalar) component of M is written (M). 
The '•' and 'A' symbols are retained for the lowest-grade and highest-grade terms 
of the series (2.11), so that 

A r -B s = (AB)| r _ s | (2.12) 
ArhB s = (AB) r+S) (2.13) 

which are called the interior and exterior products respectively. The exterior 
product is associative, and satisfies the symmetry property 

A r AB s = {-l) rs B s AA r . (2.14) 

Two further products can also be defined from the geometric product. These 
are the scalar product 

A*B=(AB) (2.15) 

and the commutator product 

AxB = \{AB - BA). (2.16) 

The scalar product (2.15) is commutative and satisfies the cyclic reordering 
property 

(A...BC) = (CA...B). (2.17) 

The scalar product (2.15) and the interior product (2.12) coincide when acting 
on two homogeneous multivectors of the same grade. The associativity of the 
geometric product ensures that the commutator product (2.16) satisfies the 
Jacobi identity 

Ax (B xC) + B x (C x A) + C x (Ax B) = 0. (2.18) 

Finally we introduce some further conventions. Throughout we employ the 
operator ordering convention that, in the absence of brackets, inner, outer, 
commutator and scalar products take precedence over geometric products. Thus 
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a -be means (a-b)c, not a -(be). This convention helps to eliminate unwieldy 
numbers of brackets. Summation convention is employed throughout, except 
for indices that denote the grade of a multivector, which are not summed over. 
Natural units (H = c = 4ireo = G = 1) are used except where explicitly stated. 
Throughout we refer to a Lorentz transformation (i.e. a spatial rotation and/or 
boost) simply as a 'rotation'. 



2.1 The spacetime algebra 

Of central importance to this paper is the geometric algebra of spacetime, the 
spacetime algebra [14]. To describe the spacetime algebra it is helpful to intro- 
duce a set of four orthonormal basis vectors {7 M }, /U = . . . 3, satisfying 

7/x-7f = «V = diag(+ ). (2.19) 

The vectors {7^} satisfy the same algebraic relations as Dirac's 7-matrices, but 
they now form a set of four independent basis vectors for spacetime, not four 
components of a single vector in an internal 'spin-space'. The relation between 
Dirac's matrix algebra and the spacetime algebra is described in more detail 
in Appendix A, which gives a direct translation of the Dirac equation into its 
spacetime algebra form. 

A frame of timelike bivectors {07.}, k = 1 . . . 3 is defined by 

<r k = 7fc7o, (2.20) 

and forms an orthonormal frame of vectors in the space relative to the 70 direc- 
tion. The algebraic properties of the {07} are the same as those of the Pauli 
spin matrices, but they now represent an orthonormal frame of vectors in space, 
and not three components of a vector in spin-space. The highest-grade element 
(or 'pseudoscalar') is denoted by I and is defined as: 

I = 7o7i7273 = a- 1 a- 2 cr i . (2-21) 

The pseudoscalar satisfies I 2 = — 1. Since we are in a space of even dimension, I 
aniicommutcs with odd-grade elements, and commutes only with even-grade el- 
ements. With these definitions, a basis for the 16-dimcnsional spacetime algebra 
is provided by 

1 {7^} Wk, I<r k } {Ij^} I - 2 22 

1 scalar 4 vectors 6 bivectors 4 trivectors 1 pseudoscalar. 

Geometric significance is attached to the above relations as follows. An 
observer's rest frame is characterised by a future-pointing timelike (unit) vector. 
If this is chosen to be the 70 direction then the 70-vector determines a map 
between spacetime vectors a = 0^7^ and the even subalgebra of the spacetime 
algebra via 

«7o = a + a, (2.23) 

where 

a =a-7o (2.24) 
a = aA7 . (2.25) 
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The 'relative vector' a can be decomposed in the {<Tk} frame and represents a 
spatial vector as seen by an observer in the 70-frame. Since a vector appears to 
an observer as a line segment existing for a period of time, it is natural that what 
an observer perceives as a vector should be represented by a spacetime bivector. 
Equation (2.23) embodies this idea, and shows that the algebraic properties of 
vectors in relative space are determined entirely by the properties of the fully 
rclativistic spacetime algebra. 

The split of the six spacetime bivectors into relative vectors and relative 
bivectors is a frame-dependent operation different observers sec different 
relative spaces. This fact is clearly illustrated with the Faraday bivector F. 
The 'space-time split' [23] of F into the 70-system is made by separating F into 
parts that anticommute and commute with 70. Thus 

F = E + IB, (2.26) 

where 

E= \{F- l0 F l0 ) (2.27) 

= i(F + 7 F 7o ). (2.28) 

Both E and B are spatial vectors in the 70-frame, and IB is a spatial bivector. 
Equation (2.26) decomposes F into separate electric and magnetic fields, and 
the explicit appearance of 70 in the formulae for E and B shows how this split 
is observer-dependent. 

The identification of the algebra of three-dimensional space with the even 
subalgebra of the spacetime algebra necessitates a convention that articulates 
smoothly between the two algebras. For this reason relative (or spatial) vectors 
in the 70-system are written in bold type to record the fact that they are in fact 
spacetime bivectors. This distinguishes them from spacetime vectors, which are 
left in normal type. Further conventions are introduced where necessary. 

2.2 Geometric calculus 

Many of the derivations in this paper employ the vector and multivector deriva- 
tives [11, 31]. Before defining these, however, we need some simple results for 
vector frames. Suppose that the set {ek} form a vector frame. The reciprocal 
frame is determined by [11] 

e J = (-l) j - 1 eiAe 2 A- ■ -AejA- ■ -AenE' 1 (2.29) 

where 

E = eiAe 2 A- • -Ae„ (2.30) 

and the check on ej denotes that this term is missing from the expression. The 
{efc} and {e k } frames are related by 

e r e h = b). (2.31) 

An arbitrary multivector B can be decomposed in terms of the {e k } frame into 

B= B i ... j e i A---Ae i (2.32) 

i<-<3 
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where 



Bi...j =B-(e i A---Ae i ). (2.33) 

Suppose now that the multivector F is an arbitrary function of some multi- 
vector argument X, F = F(X). The derivative of F with respect to X in the 
A direction is defined by 

A, dx F { X)=^ X + ^- F ^. (2.34) 
From this the multivector derivative dx is defined by 

d x = e i A---Ae j {e j A---Ae i )*d x . (2.35) 

%<■■■<] 

This definition shows how the multivector derivative dx inherits the multivector 
properties of its argument X, as well as a calculus from equation (2.34). 

Most of the properties of the multivector derivative follow from the result 
that 

dx(XA) =P X {A), (2.36) 

where Px{A) is the projection of A onto the grades contained in X. Leibniz' rule 
is then used to build up results for more complicated functions (see Appendix B). 
The multivector derivative acts on the next object to its right unless brackets are 
present; for example in the expression dxAB the dx acts only on A, but in the 
expression dx(AB) the dx acts on both A and B. If the dx is intended to act 
only on B then this is written as dxAB, the overdot denoting the multivector 
on which the derivative acts. As an illustration, Leibniz' rule can be written in 
the form 

dx{AB) = d x AB + d x AB. (2.37) 

The overdot notation neatly encodes the fact that, since dx is a multivector, it 
does not necessarily commute with other multivectors and often acts on func- 
tions to which it is not adjacent. 

We frequently make use of the derivative with respect to a general vector to 
perform a range of linear algebra operations. For such operations the following 
results are useful: 

d a a-A r = rA r (2.38) 
d a aAA r = (n- r)A r (2.39) 
d a A r a = (-l) r (n - 2r)A r , (2.40) 

where n is the dimension of the space (n = 4 for all the applications considered 
here). Further results are given in appendix B. 

The derivative with respect to spacetime position x is called the vector 
derivative, and is of particular importance. Suppose that we introduce an arbi- 
trary set of spacetime coordinates x 11 . These define a coordinate frame 

«, - ^ (2-41) 
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The reciprocal frame is denoted e^, and from this we define the vector derivative 
Vby 

V = = e"0„. (2.42) 

(The importance of the vector derivative means that it is sensible to give it 
a unique symbol.) The vector derivative inherits the algebraic properties of a 
vector in the spacetime algebra. The usefulness of the geometric product for 
the vector derivative is illustrated by electromagnetism. In tensor notation, 
Maxwell's equations are 

dpF^ = J u , d [a F^ = 0. (2.43) 

These have the spacetime algebra equivalents [14] 

V-F=J, VAF = 0. (2.44) 

But wc can utilise the geometric product to combine these into the single equa- 
tion 

VF = J. (2.45) 

The great advantage of the V operator is that it possesses an inverse, so a first- 
order propagator theory can be developed for it [11, 26]. This is not possible 
for the separate V- and VA operators. 



2.3 Linear algebra 

Geometric algebra offers many advantages over tensor calculus in developing 
the theory of linear functions [11, 34, 36]. A linear function mapping vectors to 
vectors is written in a sans-serif font, f(a), where a is an arbitrary arguments 
and plays the role of a placeholder. The argument of a linear function is usually 
assumed to be a constant vector, unless stated otherwise. The adjoint function 
is written with an overbar, f(a), and is defined such that 

o-f(6) = f(o)-6. (2.46) 

It follows that 

f(o) =d b {f(b)a) = e|t <f(&)e"). (2.47) 

We will frequently employ the derivative with respect to the vectors a, b etc. to 
perform algebraic manipulations of linear functions, as in equation (2.47). The 
advantage is that all manipulations are then frame-free. Of course, the d a and 
a vectors can easily be replaced by the sum over a set of frame vectors and their 
reciprocals, if desired. 

A symmetric function is one for which f(a) = f(a). For such functions 

d Af(a) = d a Ad b (af(b)) = f(b)Ad b . (2.48) 

It follows that for symmetric functions 

a Af(a)=O, (2.49) 

which is equivalent to the statement that f(a) = f(a). 
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Linear functions extend to act on multivcctors via 



f(aA6A. . .Ac) = f(a)Af(6) . . .Af(c), (2.50) 

so that f is now a grade-preserving linear function mapping multivcctors to 
multivcctors. In particular, since the pseudoscalar / is unique up to a scale 
factor, we can define 

dct(f) = f(Z)/" 1 . (2.51) 

Viewed as linear functions over the entire geometric algebra, f and f are 
related by the fundamental formulae 

A r -f(B a ) = f(f(A r )-B s ) r<s 

f(A r )-B s = f(A r -f(B s )) r>s, [ ' ' 

which are easily derived [11, Chapter 3]. The formulae for the inverse functions 
are found as special cases of (2.52), 

f-\A) = det(f)- 1 f(AI)I- 1 

- i 11 (2.53) 

r 1 {A) = dct(f)~ 1 r 1 f(iA). 

A number of further results for linear functions are contained in Appendix B. 
These include a coordinate-free formulation of the derivative with respect to a 
linear function, which proves to be very useful in deriving stress-energy tensors 
from action integrals. 



3 Gauge principles for gravitation 

In this section we identify the dynamical variables that will describe gravita- 
tional interactions. We start by reviewing the arguments outlined in the intro- 
duction. The basic idea is that all physical relations should have the generic form 
a(x) = b(x), where a and b are spacetime fields representing physical quantities, 
and x is the spacetime position vector. An equality such as this can certainly 
correspond to a clear physical statement. But, considered as a relation between 
fields, the physical relationship expressed by this statement is completely inde- 
pendent of where we choose to think of x as lying in spacetime. In particular, we 
can associate each position x with some new position x' = f(x) and rewrite the 
relation as a(x') = b(x'), and the equation still has precisely the same content. 
(A proviso, which will gain significance later, is that the map f(x) should be 
non-singular and cover all of spacetime.) 

A similar argument applies to rotations. The intrinsic content of a relation 
such as a(x) = b(x) at a given point Xq is unchanged if we rotate each of a and 
b by the same amount. That is, the equation Ra(xo)R = Rb(xo)R has the same 
physical content as the equation a(xo) = b(xo)- For example scalar product 
relations, from which we can derive angles, are unaffected by this change. These 
arguments apply to any physical relation between any type of multivector field. 
The principles underlying gauge theory gravity can therefore be summarised as 
follows: 

(i) The physical content of a field equation must be invariant under arbitrary 
local displacements of the fields. (This is called position-gauge invariance.) 
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(ii) The physical content of a field equation must be invariant under arbitrary 
local rotations of the fields. (This is called rotation-gauge invariance.) 

In this theory predictions for all measurable quantities, including distances 
and angles, must be derived from gauge-invariant relations between the field 
quantities themselves, not from the properties of the background spacetime. On 
the other hand, quantities that depend on a choice of 'gauge' are not predicted 
absolutely and cannot be defined operationally. 

It is necessary to indicate how this approach differs from the one adopted 
in gauge theories of the Poincare group. (This is a point on which we have 
been confused in the past [43].) Poincare transformations for a multivector field 
M (x) are defined by 

MO) h-> M' = RM(x')R (3.1) 

where 

x'=RxR + t, (3.2) 

R is a constant rotor and t is a constant vector. Transformations of this type mix 
displacements and rotations, and any attempt at a local gauging of this space- 
time symmetry fails to decouple the two [2, 9]. Furthermore, the fact that the 
rotations described by Poincare transformations include the displacement (3.2) 
(with t = 0) means that the rotations discussed under point (ii) above are not 
contained in the Poincare group. 

As a final introductory point, while the mapping of fields onto spacetime 
positions is arbitrary, the fields themselves must be well-defined. The fields 
cannot be singular except at a few special points. Furthermore, any remapping 
of the fields in the must be one-to-one, else we would cut out some region 
of physical significance. In later sections we will see that general relativity 
allows operations in which regions of spacetime are removed. These are achieved 
through the use of singular coordinate transformations and are the origin of some 
potential differences between GTG and general relativity. 

3.1 The position-gauge field 

We now examine the consequences of the local symmetries we have just dis- 
cussed. As in all gauge theories we must study the effects on derivatives, since 
all non-derivative relations already satisfy the correct requirements. 

We start by considering a scalar field <f>(x) and form its vector derivative 
V</>(x). Suppose now that from <p{x) wc define the new field (j)'(x) by 

cp'(x) = 4>{x'), (3.3) 

where 

x' = f(x) (3.4) 

and f(x) is an arbitrary (differentiable) map between spacetime position vectors. 
The map f(x) should not be thought of as a map between manifolds, or as 
moving points around; rather, the function fix) is merely a rule for relating 
one position vector to another within a single vector space. Note that the new 
function 4>'(x) is given by the old function (j> evaluated at x' . We could have 
defined things the other way round, so that <j>'(x') is given by 4>(x), but the form 
adopted here turns out to be easier to implement in practice. 
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If we now act on the new scalar field (j)' with V we form the quantity V</>(x'). 
To evaluate this we first look at the directional derivatives of <j>', 

o-V#aO = limi(^(/(x + ea)) -*(/(*))) 

= lim - (</>(x' + ef (a)) - <£(x')) 

= f(a)-V^(x'), (3.5) 

where 

f(a) = a-V/(x) (3.6) 

and the subscript on V z ' records that the derivative is now with respect to the 
new vector position variable x'. The function f(a) is a linear function of a and 
an arbitrary function of x. If we wish to make the position-dependence explicit 
we write this as f(a; x). In general, any position-dependent linear function with 
its position-dependence suppressed is to be taken as a function of x. Also — as 
stated previously — the argument of a linear function should be assumed to be 
constant unless explicitly stated otherwise. 
From (3.5) we see that 

V a =f(V x /) (3.7) 

and it follows that 

V0'(x) = f(V x ^(x'))- (3.8) 

The bracketed term on the right-hand side, V x '4>(x'), is the old gradient vector 
V<^ evaluated at x' instead of x. This tells us how to modify the derivative 
operator V: we must introduce a new linear function that assembles with V in 
such a way that the f field is removed when the full object is displaced. The 
resulting object will then have the desired property of just changing its position 
dependence under arbitrary local displacements. We therefore introduce the 
position-gauge field h(a;x), which is a linear function of a and an arbitrary 
function of position x. As usual this is abbreviated to h(a) when the position- 
dependence is taken as a function of x. Under the displacement x i— ► x' = fix), 
h(a) is defined to transform to the new field h'(a;x), where 

h'(a;x) = h(f- 1 (a);/(x)) = h(r 1 (a);;s'). (3.9) 

This ensures that h(V) transforms as 

h(y x ; x) i * h(r 1 (V a ,);ar') = h(V^;x'). (3.10) 

This transformation law ensures that, if we define a vector A{x) by 

A(x) = h(V0(x)), (3.11) 

then A(x) transforms simply as A[x) i— » A'(x) = A(x') under arbitrary dis- 
placements. This is the type of behaviour we seek. The vector A(x) can now be 
equated with other (possibly non-differentiated) fields and the resulting equa- 
tions are unchanged in form under arbitrary repositioning of the fields in space- 
time. 

Henceforth, we refer to any quantity that transforms under arbitrary dis- 
placements as 

M{x) h-> M'(x) = M(x') (3.12) 
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as behaving covariantly under displacements. The position gauge field h en- 
ables us to form derivatives of covariant objects which are also covariant under 
displacements. When we come to calculate with this theory, we will often fix 
a gauge by choosing a labelling of spacetime points with vectors. In this way 
we remain free to exploit all the advantages of representing points with vectors. 
Of course, all physical predictions of the theory will remain independent of the 
actual gauge choice. 

The h-ficld is not a connection in the conventional Yang-Mills sense. The 
coupling to derivatives is different, as is the transformation law (3.9). This is 
unsurprising, since the group of arbitrary translations is infinite-dimensional (if 
we were considering maps between manifolds then this would form the group 
of diffcomorphisms) . Nevertheless the h-field embodies the idea of replacing 
directional derivatives with covariant derivatives, so clearly deserves to be called 
a gauge field. 

A remaining question is to find the conditions under which the h-field can 
be transformed to the identity. This would be the case if h was derived entirely 
from a displacement, which would imply that 

h(a)=f _1 (a), (3.13) 

and hence 

h _1 (a)=f(a). (3.14) 
But, from the definition of f (a), it follows that 

f(a) = d b {ab-Vf(x) = V{f(x)a) (3.15) 

and hence that 

VAf(o) = VAV(/0)a) = 0. (3.16) 
So, if the h(a) field can be transformed to the identity, it must satisfy 

VAh" 1 (a) = 0. (3.17) 

An arbitrary h-field will not satisfy this equation, so in general there is no way 
to assign position vectors so that the effects of the h-field vanish. In the light 
of equations (3.16) and (3.17) it might seem more natural to introduce the 
gauge field as h _1 (V), instead of h(V). There is little to choose between these 
conventions, though our choice is partially justified by our later implementation 
of the variational principle. 

3.2 The rotation-gauge field 

We now examine how the derivative must be modified to allow rotational free- 
dom from point to point, as described in point (ii) at the start of this section. 
Here we give an analysis based on the properties of classical fields. An analysis 
based on spinor fields is given in the following section. We have already seen 
that the gradient of a scalar field is modified to k(V<f>) to achieve covariance 
under displacements. But objects such as temperature gradients are certainly 
physical, and can be equated with other physical quantities. Consequently vec- 
tors such as h(V0) must transform under rotations in the same manner as all 
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other physical fields. It follows that, under local spacetime rotations, the h-ficld 
must transform as 

h(o) i — ► Rh{a)R. (3.18) 

Now consider an equation such as Maxwell's equation, which we saw in 
Section 2.2 takes the simple spacetime algebra form VF = J. Once the position- 
gauge field is introduced, this equation becomes 

h(V)F = J, (3.19) 

where 

F=h(F) and J = det(h) h _1 (J). (3.20) 

(The reasons behind these definitions will be explained in Section 7. The use 
of a calligraphic letter for certain covariant fields is a convention we have found 
very useful.) The definitions of T and J ensure that under local rotations they 
transform as 

T i ^ RTR and J h-> RJR. (3.21) 

Any (multi)vector that transforms in this manner under rotations and is covari- 
ant under displacements is referred to as a covariant (multi)vector. 

Equation (3.19) is covariant under arbitrary displacements, and we now need 
to make it covariant under local rotations as well. To achieve this we replace 
h(V) by h^)^ and focus attention on the term d^T. Under a position- 
dependent rotation we find that 

d^RTR) = Rd^TR + d^RTR + RTd^R. (3.22) 

Since the rotor R satisfies RR = 1 we find that 

d^RR + Rd^R = (3.23) 
=S> df.RR = -Rd^R = -{d,,RR)~. (3.24) 

Hence d^RR is equal to minus its reverse and so must be a bivector. (In a 
geometric algebra the bivectors form a representation of the Lie algebra of the 
rotation group [30].) We can therefore write 

d^RTR) = Rd^TR + 2(d^RR) x {RTR). (3.25) 

To construct a covariant derivative we must therefore add a connection term to 
dfj, to construct the operator 

£V = df, + Sip x . (3.26) 

Here f2 M is a bivector-valued linear function. To make the linear argument 
explicit we write 

fi(e M ) = r> M (3.27) 

so that f2(a) = 0(a; x) is a linear function of a with an arbitrary x-dcpendcncc. 
The commutator product of a multivector with a bivector is grade-preserving 
so, even though it contains non-scalar terms, preserves the grade of the 
multivector on which it acts. 
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Under local rotations the term in T>^ cannot change, and we also expect 
that the T>^ operator be unchanged in form (this is the essence of 'minimal 
coupling'). We should therefore have 

K = d » + % x • ( 3 - 28 ) 

But the property that the covariant derivative must satisfy is 

V'^RTR) = RDuJFR (3.29) 

and, substituting (3.28) into this equation, we find that M transforms as 

M i ► n\ L = RQ^R - 2d^RR. (3.30) 

In coordinate-free form we can write this transformation law as 

0(a) i-> O'(a) = Rtt(a)R — 2a-VRR. (3.31) 

Of course, since 0(a) is an arbitrary function of position, it cannot in general 
be transformed away by the application of a rotor. Wc finally reassemble the 
derivative (3.26) with the h(9 a ) term to form the equation 

h(e")2VF = J. (3.32) 

The transformation properties of h(a), T, J and 0(a) ensure that this equation 
is now covariant under rotations as well as displacements. 

To complete the set of transformation laws, we note that under displacements 
0(a) must transform in the same way as a-Vi?i?, so that 

0(a; x) i-> 0(f (a); f(x)) = 0(f (a); x). (3.33) 

It follows that 

b(d a )Sl(a)xF(x) ^ h(f- x (9 ); x') 0(f(a); x') x T{x') 

= h(d a ;x')n(a;x')xT(x'), (3.34) 

as required for covariance under local translations. 

General considerations have led us to the introduction of two new gauge 
fields: vector- valued linear function h(a; x) and the bivector- valued linear func- 
tion 0(a;x), both of which are arbitrary functions of position x. This gives a 
total of 4x4 + 4x6 = 40 scalar degrees of freedom. The h(a) and 0(a) fields 
are incorporated into the vector derivative to form the operator h(e M )I?^, which 
acts covariantly on multivector fields. Thus we can begin to construct equations 
whose intrinsic content is free of the manner in which wc represent spacctime 
positions with vectors. We next see how these fields arise in the setting of the 
Dirac theory. This enables us to derive the properties of the operator from 
more primitive considerations of the properties of spinors and the means by 
which observables are constructed from them. First, though, let us compare 
the fields that we have defined with the fields used conventionally in general 
relativity . To do this, we first construct the vectors 

ff M = h" 1 (e M ), = h(e"). (3.35) 
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From these, the metric is defined by 



(3.36) 



Note that the inner product here is the standard spacetime algebra inner prod- 
uct and is not related to a curved space metric. Similarly, the h-field can be 
used to construct a vierbein, as discussed in Appendix C. A vierbein in general 
relativity relates a coordinate frame to an orthonormal frame. But, while the 
h-ficld can be used to construct such a vierbein, it should be clear that it serves 
a different purpose in GTG — it ensures covariance under arbitrary displace- 
ments. This was the motivation for the introduction of a vierbein in Kibble's 
work [2], although only infinitesimal transformations could be considered there. 
Furthermore, the h-ficld is essential to enable a clean separation between field 
rotations and displacements, which again is not achieved in other approaches. 

3.3 Gauge fields for the Dirac action 

We now rederive the gravitational gauge fields from symmetries of the Dirac ac- 
tion. The point here is that, once the h-field is introduced, spacetime rotations 
and phase rotations couple to the Dirac field in essentially the same way. To see 
this, we start with the Dirac equation and Dirac action in a slightly unconven- 
tional form [20, 21, 31]. We saw in Section 2 that rotation of a multivector is 
performed by the double-sided application of a rotor. The elements of a linear 
space that is closed under single-sided action of a representation of the rotor 
group arc called spinors. In conventional developments a matrix representation 
for the Clifford algebra of spacetime is introduced, and the space of column vec- 
tors on which these matrices act defines the spin-space. But there is no need to 
adopt such a construction. For example, the even subalgcbra of the spacetime 
algebra forms a linear space that is closed under single-sided application of the 
rotor group. The even subalgcbra is also an eight-dimensional linear space, the 
same number of real dimensions as a Dirac spinor, and so it is not surprising 
that a one-to-one map between Dirac spinors and the even subalgebra can be 
constructed. Such a map is given in Appendix A. The essential details are that 
the result of multiplying the column spinor by the Dirac matrix 7^ is repre- 
sented in the spacetime algebra as ip <— ► 7 At '07O; and that multiplication by the 
scalar unit imaginary is represented as ip > tylcTz- It is easily seen that these 
two operations commute and that they map even multivectors to even multi- 
vectors. By replacing Dirac matrices and column spinors by their spacetime 
algebra equivalents the Dirac equation can be written in the form 



which is now representation-free and coordinate-free. Using the same replace- 
ments, the free-particle Dirac action becomes 



and, with the techniques of Appendix B, it is simple to verify that variation of 
this action with respect to ip yields equation (3.37) with A=Q. 

It is important to appreciate that the fixed 70 and 73 vectors in (3.37) 
and (3.38) do not pick out preferred directions in space. These vectors can 
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(3.37) 




(3.38) 
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be rotated to new vectors RqJqRo and R0J3R0, and replacing the spinor by 
tJjRq recovers the same equation (3.37). This point will be returned to when we 
discuss forming observables from the spinor tp. 

Our aim now is to introduce gauge fields into the action (3.38) to ensure 
invariance under arbitrary rotations and displacements. The first step is to 
introduce the h-ficld. Under a displacement tJj transforms covariantly, 

i/>(x) 1 — ► ip'(x) = i/)(x'), (3.39) 

where x' = f(x). We must therefore replace the V operator by h(V) so that 
h(V)-0 is also covariant under translations. But this on its own does not achieve 
complete invariance of the action integral (3.38) under displacements. The ac- 
tion consists of the integral of a scalar over some region. If the scalar is re- 
placed by a displaced quantity, then we must also transform the measure and 
the boundary of the region if the resultant integral is to have the same value. 
Transforming the boundary is easily done, but the measure does require a little 
work. Suppose that we introduce a set of coordinates x M . The measure |d 4 a;| is 
then written 

\d 4 x\ = -Ie /\ei/\e 2 Ae 3 dx° dx 1 dx 2 dx 3 , (3.40) 

where 

- £■ < 3 - 41 > 

By definition, \d A x\ is already independent of the choice of coordinates, but it 
must be modified to make it position-gauge invariant. To sec how, we note that 
under the displacement x t— > f(x), the {e M } frame transforms to 

It follows that to ensure invariance of the action integral we must replace each 
of the e M by h^ 1 (e M ). Thus the invariant scalar measure is 

-Jh _1 (eo)A- • ■Ah _1 (e 3 ) dx" ■ ■ ■ dx 3 = det(h)-VH (3-43) 

These results lead us to the action 

S = J |d 4 x|dct(h)- 1 (h(V)^/7 3 Vi-TOVV ; ), (3.44) 

which is unchanged in value if the dynamical fields are replaced by 

i//(x)=i>(x') (3.45) 
h'(a;x) = h(f- 1 (a);x'), (3.46) 

and the boundary is also transformed. 

Having arrived at the action in the form of (3.44) we can now consider the 
effect of rotations applied at a point. The representation of spinors by even 
elements is now particularly powerful because it enables both internal phase 
rotations and rotations in space to be handled in the same unified framework. 
Taking the electromagnetic coupling first, we see that the action (3.44) is in- 
variant under the global phase rotation 

%l/ =4^ ICT \ (3.47) 



23 



(Recall that multiplication of \tp) by the unit imaginary is represented by right- 
sided multiplication of ip by /<x 3 .) The transformation (3.47) is a special case 
of the more general transformation 

ip i — * tpR, (3.48) 

where R is a constant rotor. Similarly invariance of the action (3.44) under 
spacetime rotations is described by 

ip h-> Rip (3.49) 
h(o) i — * i?h(a)fl. (3.50) 

In both cases, ip just picks up a single rotor. From the previous section we 
know that, when the rotor R is position-dependent, the quantity d^RR is a 
bivector- valued linear function of e M . Since 

d M (i#) = Rd^j + (d^RRjRiP, (3.51) 

with a similar result holding when the rotor acts from the right, we need the 
following covariant derivatives for local internal and external rotations: 

Internal: D^ip = d^ip + (3.52) 
External: D^ip = d^ip + ^%rp. (3.53) 

For the case of (internal) phase rotations, the rotations are constrained to take 
place entirely in the ier 3 plane. It follows that the internal connection takes the 
restricted form 

n T (a) = 2ea-AI<r 3 (3.54) 

where A is the conventional electromagnetic vector potential and e is the cou- 
pling constant (the charge). The full covariant derivative therefore has the form 

h(e»)(d^ + + eipiase^-A) (3.55) 

and the full invariant action integral is now (in coordinate-free form) 

S = J |d 4 x|(deth)- 1 ^h(a Q )(a-V + if7(a))^7 3 ^-eh(A)V'7o^-TOV'V')- (3.56) 

The action (3.56) is invariant under the symmetry transformations listed in 
Table 1. 

3.4 The coupled Dirac equation 

Having arrived at the action (3.56) we now derive the coupled Dirac equation 
by extremising with respect to ip, treating all other fields as external. When 
applying the Euler-Lagrange equations to the action (3.56) the ip and ip fields 
are not treated as independent, as they often are in quantum theory Instead, 
we just apply the rules for the multivector derivative discussed in Section 2.2 
and Appendix B. The Eulcr-Lagrange equations can be written in the form 

8^ = 8^^) (3.57) 
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Local 
symmetry 


if/(x) 


Transformed Fields 
h'(a;x) fl'(a;x) 


eA'(x) 


Displacements 

Spacetime 
rotations 

Phase 
rotations 


yj(x') 


h(f- 1 (a);s') 
Rh(a)R 

h(a) 


n(f(a);x / ) 
Rn(a)R-2a-X7RR 

fi(a) 


ef(A(x')) 

eA 
eA - V0 



Table 1: TTie symmetries of the action integral (3.56). 



where for simplicity we assume that the a; M are a set of fixed orthonormal co- 
ordinates given by x M = 7 A '-x (see Appendix B). Applied to the action (3.56), 
equation (3.57) yields 

(h(V)V>/73)~ + !/73#(d a )fi(a) + i(h(5 a )f2(a)V/73)~ - e 7o #(A) 

- (eh(A)V>7o)~ - 2mi> = 9 A1 (det(h)- 1 / 73 ^h(7^)) det(h). (3.58) 

Reversing this equation and simplifying gives 

h(0„)(a-V + iQ(a))V/73 - eh(A)V>7o - m^> 

= -idet(h)P M (h(7^)det(h)- 1 )^73, (3.59) 

where we have employed the T>^ derivative defined in equation (3.26). If we now 
introduce the notation 

Dip = h(5 )(a-V + fo(a))fl> = + \%)\l> (3.60) 
A = h(A), (3.61) 

we can write equation (3.59) in the form 

Diplas - eAtp = m^ 7o - | det(h)D (tl (h(7^) dct(h)- 1 )^/cr3. (3.62) 

This equation is covariant under the symmetries listed in Table 1, as must be 
the case since the equation was derived from an invariant action integral. But 
equation (3.62) is not what we would have expected had we applied the gauging 
arguments at the level of the Dirac equation, rather than the Dirac action. 
Instead, we would have been led to the simpler equation 

Diplcr 3 — eAip = rnipj . (3.63) 

Clearly, equation (3.62) reduces to equation (3.63) only if the h(a) and 0(a) 
fields satisfy the identity 

dct(h)P Al (h( 7 ' I )det(h)- 1 ) =0. (3.64) 

The restriction to an orthonormal coordinate frame can be removed by writing 
this equation as 

det(h)P AI (h(e^)dct(h)- 1 ) =0. (3.65) 
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(It is not hard to show that the left-hand side of equation (3.65) is a covari- 
ant vector; later it will be identified as a contraction of the 'torsion' tensor). 
There are good reasons for expecting equation (3.65) to hold. Otherwise, the 
minimally-coupled Dirac action would not yield the minimally-coupled equation, 
which would pose problems for our use of action principles to derive the gauged 
field equations. We will see shortly that the demand that equation (3.65) holds 
places a strong restriction on the form that the gravitational action can take. 

Some further comments about the derivation of (3.62) are now in order. The 
derivation employed only the rules of vector and multivector calculus applied to 
a 'flat-space' action integral. The derivation is therefore a rigorous application of 
the variational principle. This same level of rigour is not always applied when 
deriving field equations from action integrals involving spinors. Instead, the 
derivations are often heuristic — \ip) and {ip \ are treated as independent variables 
and the is just 'knocked off' the Lagrangian density to leave the desired 
equation. Furthermore, the action integral given by many authors for the Dirac 
equation in a gravitational background has an imaginary component [44, 45] , in 
which case the status of the variational principle is unclear. To our knowledge, 
only Hchl & Datta [9] have produced a derivation that in any way matches 
the derivation produced here. Hchl & Datta also found an equation similar 
to (3.65), but they were not working within a gauge theory setup and so did not 
comment on the consistency (or otherwise) of the minimal-coupling procedure. 

3.5 Observables and covariant derivatives 

As well as keeping everything within the real spacetime algebra , representing 
Dirac spinors by elements of the even subalgebra offers many advantages when 
forming observables. As described in Appendix A, observables are formed by the 
double-sided application of a Dirac spinor ip to some combination of the fixed 
{7^} frame vectors. So, for example, the charge current is given by J = ipjoip 
and the spin current by s = '073V'- I n general, an observable is of the form 

M = ipYip, (3.66) 

where T is a constant multivector formed from the {7 M }- AH observables are 
invariant under phase rotations, so T must be invariant under rotations in the 
la 3 plane. Hence T can consist only of combinations of 70, 73, Ier 3 and their 
duals (formed by multiplying by /). An important point is that, in forming the 
observable M, the T multivector is completely 'shielded' from rotations. This is 
why the appearance of the 70 and 73 vectors on the right-hand side of the spinor 
ip in the Dirac action (3.56) does not compromise Lorentz invariance, and does 
not pick out a preferred direction in space [21]. All observables are unchanged 
by rotating the {7^} frame vectors to Ro^^Ro and transforming ip to ipRo- (In 
the matrix theory this corresponds to a change of representation.) 

Under translations and rotations the observables formed in the above man- 
ner (3.66) inherit the transformation properties of the spinor ip. Under transla- 
tions the observable M = ipTip therefore transforms from M(x) to M(x'), and 
under rotations M transforms to RipTipR = RMR. The observable M is there- 
fore covariant. These Dirac observables are the first examples of quantities that 
transform covariantly under rotations, but do not inherit this transformation 
law from the h-ficld. In contrast, all covariant forms of classical fields, such as 
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T or the co-variant velocity along a worldlinc h _1 (i), transform under rotations 
in a manner that that is dictated by their coupling to the h-field. Classical 
general relativity in fact removes any reference to the rotation gauge from most 
aspects of the theory. Quantum theory, however, demands that the rotation 
gauge be kept in explicitly and, as we shall show in Section 8, Dirac fields probe 
the structure of the gravitational fields at a deeper level than classical fields. 
Furthermore, it is only through consideration of the quantum theory that one 
really discovers the need for the rotation-gauge field. 

One might wonder why the observables are invariant under phase rotations, 
but only covariant under spatial rotations. In fact, the h-ficld enables us to form 
quantities like h(JVf), which are invariant under spatial rotations. This gives an 
alternative insight into the role of the h-field. We will find that both covariant 
observables (M) and their rotationally-invariant forms (h(M) and h _1 (M)) play 
important roles in the theory constructed here. 

If we next consider the directional derivative of M, we find that it can be 
written as 

d»M = {d^)T^ + ^r(5 M V)~- (3-67) 

This immediately tells us how to turn the directional derivative d^M into a co- 
variant derivative: simply replace the spinor directional derivatives by covariant 
derivatives. Hence we form 

= M (^r$+ty 1 x(V>r$. (3.68) 

We therefore recover the covariant derivative for observables: 

£> M M = d^M + % x M. (3.69) 

This derivation shows that many features of the 'classical' derivation of gravi- 
tational gauge fields can be viewed as arising from more basic quantum trans- 
formation laws. 

Throughout this section we have introduced a number of distinct gravita- 
tional covariant derivatives. We finish this section by discussing some of their 
main features and summarising our conventions. The operator T> ^ acts on any 
covariant multivector and has the important property of being a derivation, that 
is it acts as a scalar differential operator, 

V^AB) = {V, L A)B + AfppB). (3.70) 

This follows from Leibniz' rule and the identity 

f2 M x(AB) = Sl fl xAB + An^xB. (3.71) 

Neither or T>^ are fully covariant, however, since they both contain the fl^ 
field, which picks up a term in f under displacements (3.33). It is important in 
the applications to follow that we work with objects that are covariant under 
displacements, and to this end we define 

Lu{a) = Q(h(a)) = a-g^ M . (3.72) 
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Gauge fields 



Displacements: h(a) 

Rotations: 0(a), u)(a) = fi(h(a)) 



Spinor derivatives 



a-DV = a-h(V)V' + ^u{a)ip 



Obscrvables derivatives a-VM — a-h(V)M + w(a)xM 

VM = daa-VM = V-M + DAM 



Vector derivative 



V = 



d_ 



Multivector derivative d x = V] e* A • • ■ A e° ' (e 3 -A ■ ■ ■ A ) * d x 



K-<j 



Table 2: Definitions and conventions 



We also define the full covariant directional derivatives a-D and a- 1? by 

a-Lh/' = a-h(V)^ + ^(a)V> (3.73) 
a-VM = a-h(V)M + lj(o)xM. (3.74) 

Note that these conventions imply that V^ = g^-V, with a similar result holding 
for D^. 

For the a-V operator we can further define the covariant vector derivative 

VM = d a a-VM = g^V^M. (3.75) 

The covariant vector derivative contains a grade-raising and a grade-lowering 
component, so that 

VA = V-A + VAA, (3.76) 



where 



V-A = d a -{a-VA)=g»-{V ll A) 
VAA = d a A{a-VA) = g^AiV^A). 



(3.77) 
(3.78) 



As with the vector derivative, V inherits the algebraic properties of a vector. 

A summary of our definitions and conventions is contained in Table 2. We 
have endeavoured to keep these conventions as simple and natural as possible, 
but a word is in order on our choices. It should be clear that it is a good 
idea to use a separate symbol for the spacetime vector derivative (V), as op- 
posed to writing it as d x . This maintains a clear distinction between spacetime 
derivatives, and operations on linear functions such as 'contraction' (d a -) and 
'protraction' (d a A). It is also useful to distinguish between spinor and vector 
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covariant derivatives, which is why we have introduced separate D and T> sym- 
bols. We have avoided use of the d symbol, which already has a very specific 
meaning in the language of differential forms. It is necessary to distinguish be- 
tween rotation-gauge derivatives (£> M ) and the full covariant derivative with the 
h-ficld included (a-T>). Using 2? M and a- T> for these achieves this separation in 
the simplest possible manner. Throughout this paper we will switch between 
using coordinate frames (as in 2? M ) and the frame- free notation (such as u>(a)). 
Seeing which notation is appropriate for a given problem is something of an art, 
though we hope to convey some of the basic rules in what follows. 

4 The field equations 

Having introduced the h and il-ficlds, we now look to construct an invariant 
action integral that will provide a set of gravitational field equations. We start 
by defining the field-strength via 

[D M D v ]i(i=^R ltv i(i, (4.1) 

so that 

V = a„fi„ - <9„fV + VL^ x n„. (4.2) 
It follows that we also have 

[Dp, V V \M = R M „ x M. (4.3) 
A frame-free notation is introduced by writing 

FV = R(e M Ae„). ( 4 - 4 ) 

The field R(a A b) is a bivector- valued linear function of its bivector argument 
a A b. Its action on bivectors extends by linearity to the function R(B), where 
B is an arbitrary bivector and therefore, in four dimensions, not necessarily a 
pure 'blade' a A b. Where required, the position dependence is made explicit by 
writing R(B;x). 

Under an arbitrary rotation, the definition (4.2) ensures that R(B) trans- 
forms as 

R(B) i ^ R'(-B) = RR(B)R. (4.5) 
Under local displacements we find that 

R\e^e u ) = d^\e y ) - $,fi'(e„) + Q'(e u ) x n'(e^) 

= f(e ll )-V x >Sl(f(e v );x') - f(e„)-V*n(f(e M );x / ) + fi'( e/1 ) xft'(e„) 

+ fl(d^(e u )-dj(e fl );x') 
= R(f(e ^1 Ae ! ,);a: / )+0(9 ^1 f(e I/ )-^f(e ^1 );x , ). (4.6) 

But we know that 

d^(e v ) - ft,f(e M ) - d^f(x) - ft,fl^/(x) = 0, (4.7) 

so the field strength has the simple displacement transformation law 

R(B) h-> R'(B) = R(f(B);x'). (4.8) 
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A covariant quantity can therefore be constructed by denning 



11(B) = R(h(B)). (4.9) 

Under arbitrary displacements and local rotations, TZ-(B) has the following trans- 
formation laws: 

Translations: K'(B; x) = 1Z(B; x') 

Rotations: 1Z'(B) = R1Z(RBR)R. ^ ' 

We refer to any linear function with transformation laws of this type as a co- 
variant tensor. 1Z(B) is our gauge theory analogue of the Riemann tensor. We 
have started to employ a notation that is very helpful for the theory developed 
here. Certain covariant quantities, such as ~R-(B) and 2?, are written with cal- 
ligraphic symbols. This helps keep track of the covariant quantities and often 
enables a simple check that a given equation is gauge covariant. It is not nec- 
essary to write all covariant objects with calligraphic symbols, but it is helpful 
for common objects such as 71(B) and V. 

From TZ(B) we define the following contractions: 

Ricci Tensor: 11(b) = d a -TZ(aAb) (4.11) 
Ricci Scalar: K = d a -11(a) (4.12) 
Einstein Tensor: Q(a) = 11(a) - \aK. (4.13) 

The argument of 1Z determines whether it represents the Riemann or Ricci 
tensors or the Ricci scalar. Furthermore, we never apply the extension notation 
of equation (2.50) to 1Z(a), so 1Z(a A b) unambiguously denotes the Riemann 
tensor. Both 1Z(a) and Q(a) are also covariant tensors, since they inherit the 
transformation properties of 1Z(B). 

The Ricci scalar is invariant under rotations, making it our first candidate 
for a Lagrangian for the gravitational gauge fields. We therefore suppose that 
the overall action integral is of the form 

S = J \d A x\det(h)- 1 (\K- KC rn ), (4.14) 

where C m describes the matter content and k = 8t:G. The independent dynam- 
ical variables are h(a) and £l(a), and in terms of these 

U = (h(d b Ad a )(a-Vn(b) - 6-Vf2(a) + 0(a)xf2(6))>. (4.15) 

We also assume that C m contains no second-order derivatives, so that h(a) and 
Q(a) appear undifferentiated in the matter Lagrangian. 

4.1 The h-equation 

The h-ficld is undifferentiated in the entire action, so its Euler-Lagrange equa- 
tion is simply 

9 F(a) (dct(h)- 1 (^/2 - nC m )) = 0. (4.16) 
Employing some results from Appendix B, we find that 

^ (Q) dct(h)- 1 = -dct(h)- 1 h- 1 (a) (4.17) 
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and 



d- h{a) ii = d- h(a) (h(d c Ad b )R(bAc)) 

= 2h(d b )-R(bAa), (4.18) 

so that 

a F(a) (^dct(h)- 1 ) = 2^(h- 1 (a))det(h)- 1 . (4.19) 
If we now define the covariant matter stress-energy tensor T(a) by 

det(h)9 fl(Q) (/: m dct(h)- 1 ) =T(h- 1 (a)), (4.20) 

we arrive at the equation 

G(a) = kT{o). (4.21) 

This is the gauge theory statement of Einstein's equation. Note, however, that 
as yet nothing should be assumed about the symmetry of G{a) or T(a). In 
this derivation only the gauge fields have been varied, and not the properties 
of spacetime. Therefore, despite the formal similarity with the Einstein equa- 
tions of general relativity, there is no doubt that we are still working in a flat 
spacetime. 

4.2 The fi-equation 

The Euler-Lagrange field equation from f2(a) is, after multiplying through by 
det(h), 

d n{a) Tl - det(h)9 6 -V(5 n(a) ^dctth)- 1 ) = 2 K d Q(a) C m , (4.22) 

where we have made use of the assumption that Q(a) does not contain any 
coupling to matter through its derivatives. The derivatives <9fi(a) and <9n( a ) b are 
defined in Appendix B. The only properties required for this derivation are the 
following: 

d nia) {Q(b)M) =a-b(M} 2 (4.23) 
dn(b), a (c-Vfi(d)M) =a-cb-d(M} 2 . (4.24) 

From these we derive 

3n(a)(h(<9 d A<9 c Mc)xft(d)) = tt(d)xh(d d Aa) + h(aAd c )xn(c) 

= 2n(b)xh(d b Aa) (4.25) 

and 

d n{a) b (h(d d Ad c )(c-\7n(d) - d-Vn(c))) = h(aAb) - h(bAa) 

= 2h(aA6). (4.26) 

The right-hand side of (4.22) defines the 'spin' of the matter, 

S(a) = d Q{a) C m , (4.27) 

where S(a) is a bivcctor-valucd linear function of a. Combining (4.22), (4.25) 
and (4.26) yields 

h(V)Ah(a) +det(h)cVV(h(fo) dct(h)" 1 ) Ah(a) 

+ n(d b )xh(bAa) = KS(a). (4.28) 
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Recall that in an expression such as this, the vectors a and b are viewed as being 
independent of position. 

To make further progress we contract equation (4.28) with h _1 (<9 a ). To 
simplify this we require the result 

(&-Vh(a)h- 1 (9 a )) = det(h)- 1 (&-Vh(S a )5h (o) det(h)) 

= det(h) _1 &-Vdet(h). (4.29) 

We then obtain 

-2det(h)&.V(h(&)det(h) _1 ) - 2fi(«9 6 )-n(&) = k^ 1 (d a )- S(a). (4.30) 

To convert this into its manifestly covariant form we first define the covariant 
spin tensor S(a) by 

5(a) =S , (h" 1 (a)). (4.31) 

We can now write 

dct(h)I) AI (h(e' i )dct(h)- 1 ) = -±Kd a -S(a). (4.32) 

In Section 3.3 we found that the minimally-coupled Dirac action gave rise to 
the minimally-coupled Dirac equation only when the term on the left of this 
equation was zero. We now see that this requirement amounts to the condition 
that the spin tensor has zero contraction. But, if we assume that the f2(a) field 
only couples to a Dirac fermion field, then the coupled Dirac action (3.56) gives 

5(a) =S-a, (4.33) 

where 5 is the spin trivector 

S = ±V^73^- (4.34) 
In this case the contraction of the spin tensor does vanish: 

d a -(S-a) = (d a Aa)-S = 0. (4.35) 

There is a remarkable consistency loop at work here. The Dirac action gives 
rise to a spin tensor of just the right type to ensure that the minimally-coupled 
action produces the minimally-coupled equation. But this is only true if the 
gravitational action is given by the Ricci scalar\ No higher-order gravitational 
action is consistent in this way. So, if we demand that the minimally-coupled 
field equations should be derivable from an action principle, we are led to a 
highly constrained theory. This rules out, for example, the type of 'R + R 2 ' 
Lagrangian often considered in the context of Poincare gauge theory [46, 47, 
48]. In addition, the spin sector is also tightly constrained. Satisfyingly, these 
constraints force us to a theory that is first-order in the derivatives of the fields, 
keeping the theory on a similar footing to the Dirac and Maxwell theories. 

The only freedom in the action for the gravitational fields is the possible in- 
clusion of a cosmological constant A. This just enters the action integral (4.14) 
as the term — Adet(h) -1 . The presence of such a term does not alter equa- 
tion (4.28), but changes (4.21) to 

G(a) — Aa = kT(o). (4.36) 
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The presence of a cosmological constant cannot be ruled out on theoretical 
grounds alone, and this constant will be included when we consider applications 
to cosmology. 

Given that the spin is entirely of Dirac type, equation (4.28) now takes the 
form 

X>Ah(a) = d b A(b-Vh(a)) = «5-h(a). (4.37) 

This is the second of our gravitational field equations. Between them, equa- 
tions (4.21) and (4.37) define a set of 40 scalar equations for the 40 unknowns 
in h(et) and Q(a). Both equations are manifestly covariant. In the spin-torsion 
extension of general relativity (the Einstein-Cartan-Sciama-Kibble theory), 
V A h(a) would be identified as the gravitational torsion, and equation (4.37) 
would be viewed as identifying the torsion with the matter spin density. Of 
course, in GTG torsion is not a property of the underlying spacctimc, it simply 
represents a feature of the gravitational gauge fields. Equation (4.37) generalises 
to the case of an arbitrary vector a = a(x) as follows: 

VAh(a) = h(VAa) + nSMa). (4.38) 

Of particular use is the case where a is the coordinate frame vector e^, so that 
we can write 

PA/ = nS-g 11 . (4.39) 



4.3 Covariant forms of the field equations 

For all the applications considered in this paper the gravitational fields are 
generated by matter fields with vanishing spin. So, to simplify matters, we 
henceforth set iS to zero and work with the second of the field equations in the 
form 

X>Ah(a) = h(VAo). (4.40) 

It is not hard to make the necessary generalisations in the presence of spin. 
Indeed, even if the spin-torsion sector is significant, one can introduce a new 
field [49] 

u}'(a)=w(a)-±Ka-S (4.41) 

and then the modified covariant derivative with uj{a) replaced by a/ (a) still 
satisfies equation (4.40). 

The approach we adopt in this paper is to concentrate on the quantities 
that are covariant under displacements. Since both h(V) and tu(a) satisfy this 
requirement, these are the quantities with which we would like to express the 
field equations. To this end we define the operator 

L a =a-h(V) (4.42) 

and, for the remainder of this section, the vectors a, b etc. are assumed to be 
arbitrary functions of position. From equation (4.40) wc write 

h(V)Ah(c) = -d d A(u)(d)-h(c)) 
=> (6Aah(V)Ah(c)) = -(bAad d A(w(d)-h(c))) 

(L a h(b)-L b h(a))-c= (a-u(b)-b-u(a))-h(c) (4.43) 
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where, as usual, the overdots determine the scope of a differential operator, ft 
follows that the commutator of L a and L b is 

[L ai L b ] = (L a h(b) - L b h(a)) -V 

= (L a h(b)-L b h(a))-V+(L a b-L b a)-h(V) 

= (a-uj{b) - b-uj(a) + L a b - L b a) -h(V). (4.44) 

We can therefore write 

[L a ,L b ]=L c , (4.45) 

where 

c = a-cu(b) - b-uj(a) + L a b - L b a = a-Vb - b-Va. (4.46) 

This 'bracket' structure summarises the intrinsic content of (4.37). 

The general technique we use for studying the field equations is to let uj{a) 
contain a set of arbitrary functions, and then use (4.46) to find relations between 
them. Fundamental to this approach is the construction of the Riemann tensor 
1Z(B), which contains a great deal of covariant information. From the definition 
of the Riemann tensor (4.2) we find that 

K(aAb) = L a n(h(b)) - L b Q(h(a)) +uj(a)xu(b) 

= L a uj(b) - L b Lu{a) +u(a)xuj(b) - Q(L a h{b) - L b h(a)), (4.47) 

hence 

K(aAb) = L a Lu(b) - L b Lu(a) + uj(a) x uj(b) - u(c), (4.48) 

where c is given by equation (4.46). Equation (4.48) now enables 7Z(B) to be 
calculated in terms of position-gauge covariant variables. 

Solution of the 'wedge' equation 

Equation (4.40) can be solved to obtain uj(a) as a function of h and its deriva- 
tives. We define 

H{a) = h(VAfi _1 (o)) = -h(V)Ah h-^a), (4.49) 
so that equation (4.40) becomes 

d b A{co(b)-a)=H(a). (4.50) 
We solve this by first 'protracting' with d a to give 

d a hd b h{uj{b)-a) = 2d b Auj(b) = d b AH(b). (4.51) 
Now, taking the inner product with a again, we obtain 

Lu(a) - d b A(a-Lu{b)) = \a-(d b AH(b)). (4.52) 
Hence, using equation (4.50) again, we find that 

uj(a) = -H(a) + \a- (d b AH(b)). (4.53) 
In the presence of spin the term ^na-S is added to the right-hand side. 
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4.4 Point-particle trajectories 

The dynamics of a fcrmion in a gravitational background are described by the 
Dirac equation (3.63) together with the quantum-mechanical rules for construct- 
ing observables. For many applications, however, it is useful to work with clas- 
sical and semi-classical approximations to the full quantum theory. The full 
derivation of the semi-classical approximation is described in [34], but the es- 
sential idea is to specialise to motion along a single streamline defined by the 
Dirac current ^lo^- Thus the particle is described by a trajectory x(X), to- 
gether with a spinor tp(X) which contains information about the velocity and 
spin of the particle. The covariant velocity is h -1 (:c') where, for this and the 
following subsection, dashes are used to denote the derivative with respect to A. 
The covariant velocity is identified with ipjoip and the Lagrange multiplier p is 
included in the action integral to enforce this identification. Finally, an cinbein 
e is introduced to ensure reparameterisation invariance. The resultant action is 

S = J dX (ip'IeT 3 '4> + ^Q(x')ipI(T 3 '4> + p(v - me?/>7o'0) + m 2 e) , (4.54) 

where 

v = h- 1 (x'). (4.55) 

The equations of motion arising from (4.54) are discussed elsewhere [34]. (An 
effect worth noting is that, due to the spin of the particle, the velocity v and 
momentum p are not collinear.) 

We can make a full classical approximation by neglecting the spin (dropping 
all the terms containing ip) and replacing ip^/oip by p/m. This process leads to 
the action 

S = J dX(p-b- 1 {x')- \e(p 2 -to 2 )). (4.56) 

The equations of motion derived from (4.56) are 

v = ep (4-57) 

p 2 = m 2 (4.58) 

dx^ip) = Vp-h^V)- (4-59) 

The final equation yields 

a A p=h(v)p'h- 1 (.T')-i''vh(h- 1 (j ) )) 
^((VAfrV))'^)) 

= H(p)-v (4.60) 

where H(a) is defined by equation (4.49) and we have used equation (2.52). 
From equation (4.53) we see that a-Lu(a) = —a-H(a), hence 

ed\(v/e) = —uj(v)-v. (4.61) 

This is the classical equation for a point-particle trajectory. It takes its simplest 
form when A is the proper time r along the trajectory. In this case e = 1/m 
and the equation becomes 

^ = -w(t>)-t>, (4.62) 
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or, in manifestly covariant form, 



v-Vv = Q. (4.63) 

This equation applies equally for massive particles (v 2 = 1) and photons (v 2 = 
0). Since equation (4.63) incorporates only gravitational effects, any deviation 
oiv-T>v from zero can be viewed as the particle's acceleration and must result 
from additional external forces. 

Equation (4.63) is usually derived from the action 

dx^ dx u \ 1/2 



d\ v/h-i^') 2 = mj dX U^^-J , (4.64) 



which is obtained from (4.56) by eliminating p and e with their respective equa- 
tions of motion. A Hamiltonian form such as (4.56) is rarely seen in conventional 
general relativity, since its analogue would require the introduction of a vier- 
bein. Despite this, the action (4.56) has many useful features, especially when it 
comes to extracting conservation laws. For example, contracting equation (4.59) 
with a constant vector a yields 

d x (a-h- 1 {p)) ^o.Vp-r 1 ^'). (4.65) 

It follows that, if the h-field is invariant under translations in the direction a, 
then the quantity a - h -1 (p) is conserved. In Section 5.5 we show that this result 
extends to the case where h _1 (a) is a Killing vector. 



4.5 The equivalence principle and the Newtonian limit 

In the preceeding section we derived the equation v-T>v = from the classical 
limit of the Dirac action. This equation is the GTG analogue of the geodesic 
equation (see Appendix C). Arriving at such an equation shows that GTG 
embodies the weak equivalence principle — the motion of a test particle in a 
gravitational field is independent of its mass. The derivation also shows the 
limitations of this result, which only applies in the classical, spinlcss approxi- 
mation to quantum theory. The strong equivalence principle, that the laws of 
physics in a freely-falling frame are (locally) the same as those of special relativ- 
ity, is also embodied in GTG through the application of the minimal coupling 
procedure. Indeed, it is clear that both of these 'principles' are the result of 
the minimal coupling procedure. Minimal coupling ensures, through the Dirac 
equation, that point-particle trajectories are independent of particle mass in the 
absence of other forces. It also tells us how the gravitational fields couple to any 
matter field. As we have seen, this principle, coupled with the requirement of 
consistency with an action principle, is sufficient to specify the theory uniquely 
(up to an unspecified cosmological constant). 

The relationship between the minimal coupling procedure (the gauge princi- 
ple) and the equivalence principle has been widely discussed previously. Feyn- 
man [50], for example, argues that a more exact version of the equivalence 
between linear acceleration and a gravitational field requires an equation of the 
form 

gravity' = gravity + acceleration, (4.66) 



36 



which clearly resembles a gauge transformation. What is not often stressed is 
the viewpoint presented here, which is that if gravity is constructed entirely 
as a gauge theory, there is no need to invoke the equivalence principle; the 
physical effects embodied in the principle are simply consequences of the gauge 
theory approach. This further illustrates the different conceptual foundations 
of GTG and general relativity. Similarly, there is no need for the principle 
of general covariance in GTG, which is replaced by the requirement that all 
physical predictions be gauge- invariant. It is often argued that the principle 
of general covariance is empty, because any physical theory can be written in 
a covariant form. This objection cannot be levelled at the statement that all 
physical predictions must be gauge invariant, which has clear mathematical and 
physical content. 

The simplest, classical measurements in GTG are modelled by assuming 
that observers can be treated as frames attached to a single worldlinc. If this 
worldline is written as x(A), then the covariant velocity is v = h~ 1 (x / ), and 
the affine parameter for the trajectory is that which ensures that v 2 = 1. The 
affine parameter models the clock time for an observer on this trajectory. Of 
course, there are many hidden assumptions in adopting this as a realistic model 
— quantum effects are ignored, as is the physical extent of the observer — but 
is is certainly a good model in weak fields. In strong fields a more satisfactory 
model would involve solving the Dirac equation to find the energy levels of an 
atom in the gravitational background and use this to model an atomic clock. 

Equation (4.63) enables us to make classical predictions for freely-falling tra- 
jectories in GTG, and the photon case (v 2 = 0) can be used to model signalling 
between observers. As an example, consider the formula for the redshift induced 
by the gravitational fields. Suppose that a source of radiation follows a world- 
line x\(ti), with covariant velocity V\ = h~ 1 (x' 1 ). The radiation emitted follows 
a null trajectory with covariant velocity u (u 2 = 0). This radiation is received 
by an observer with a worldlinc £2(72) and covariant velocity v 2 = h _1 (x 2 ). 
The spectral shift z is determined by the ratio of the frequency observed at the 
source, u(x\)-v\, and the frequency observed at the receiver, u{x2)-V2, by 

l + z = 4^1. (4.67) 

U{X 2 j-V2 

This quantity is physically observable since the right-hand side is a gauge- 
invariant quantity. This is because each of the four vectors appearing in (4.67) 
is covariant, which eliminates any dependence on position gauge, and taking the 
dot product between pairs of vectors eliminates any dependence on the rotation 
gauge. 

A final point to address regarding the foundations of GTG is the recovery 
of the Newtonian limit. Derivations of this are easily produced by adapting the 
standard works in general relativity . Furthermore, in Section 6.3 we show that 
the description of the gravitational fields outside a static, spherically symmetric 
star is precisely the same as in general relativity. The trajectories defined by 
equation (4.63) are those predicted by general relativity, so all of the predictions 
for planetary orbits (including those for binary pulsars) are unchanged. Simi- 
larly, the results for the bending of light are the same as in general relativity. As 
we show in the applications, any differences between GTG and general relativity 
emerge through the relationship with quantum theory, and through the global 
nature of the gauge fields in GTG. These differences have no consequences for 
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classical tests of general relativity, though they are potentially testable through 
the interaction with quantum spin, and are certainly significant for discussing 
more fundamental aspects of gravitational physics. 

5 Symmetries, invariants and conservation laws 

Having determined the gravitational gauge fields and their field equations, we 
now establish some general results which are applied in the sections that follow. 
Again, we restrict to the case of vanishing torsion. The approach we adopt in 
solving the field equations is to let 10(a) be an arbitrary function, and then work 
with a set of abstract first-order equations for the terms that comprise u>(a). 
However, in letting w(a) be an arbitrary function, we lose some of the informa- 
tion contained in the 'wedge' equation (4.40). This information is recovered by 
enforcing various properties that the fields must satisfy, including the symmetry 
properties of 7t(B) and the Bianchi identities. In addition, it is often necessary 
to enforce some gauge-fixing conditions. For the rotation gauge these conditions 
are applied by studying 1Z(B), so it is important to analyse its general structure. 
We start with the result that, for an arbitrary multivector A(x), 

VA{VAA) = 2?A(h(VAh" 1 (yl))) = h(V A V Ah~ 1 (A)) = 0, (5.1) 

where we have made use of equation (4.40). It then follows from the result 

= §h(e")Ah(e , ')A(R(e, 1 Ae I ,)xA), (5.2) 

that, for any multivector A, 

d a Ad b A(K(aAb)xA) = 0. (5.3) 

This derivation illustrates a useful point. Many derivations can be performed 
most efficiently by working with the 2?^, since these contain commuting partial 
derivatives. However, the final expressions take their most transparent form 
when the h-field is included so that only fully covariant quantities are employed. 

If we now set A in equation (5.3) equal to a vector c, and protract with d Cl 
we find that 

d c Ad a Ad b A{K{aAb)xc) = -2d a Ad b AK(aAb) = 0. (5.4) 

Taking the inner product of the term on the right-hand side with c we obtain 

c-(d a Ad b All{aAb)) = d b ATZ{cAb) - d a ATZ{aAc) - d a Ad b A{TZ(aAb) x c), (5.5) 

in which both the left-hand side and the final term on the right-hand side vanish. 
We arc therefore left with the simple expression 

d a AK{aAb) = 0, (5.6) 

which summarises all the symmetries of 1Z(B). This equation says that the 
trivector d a A lZ(a A b) vanishes for all values of the vector 6, so gives a set 
of 4 x 4 = 16 equations. These reduce the number of independent degrees of 
freedom in TZ-(B) from 36 to the expected 20. It should be clear from the ease 
with which the degrees of freedom arc calculated that the present geometric 
algebra formulation has many advantages over traditional tensor calculus. 
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5.1 The Weyl tensor 

A good example of the power of the present approach is provided by an analysis 
of the Ricmann and Weyl tensors. To illustrate this point a number of examples 
of 72(7?) for physical systems are included in this section. (These are stated here 
without derivation.) The first application of geometric algebra to the analysis of 
the Riemann tensor in classical differential geometry was given by Hestenes & 
Sobczyk [11, 15]. Here the formalism is developed and extended for applications 
relevant to our gauge theory of gravity. 

Six of the degrees of freedom in 72(7?) can be removed by arbitrary gauge 
rotations. It follows that 72(7?) can contain only 14 physical degrees of freedom. 
To see how these are encoded in 72(7?) we decompose it into Weyl and 'matter' 
terms. Since the contraction of 72(a A 6) results in the Ricci tensor 72(a), we 
expect that 72(aA6) will contain a term in 72(a)A6. This must be matched with 
a term in a A 72. (6), since it is only the sum of these that is a function of a A b. 
Contracting this sum we obtain 

d a - (72(a) A6 + aA72(6)) = bK - 72(6) + 472(6) - 72(6) 

= 272(6) + 672, (5.7) 

and it follows that 

d a -{\ (72(a)A6 + aA72(6)) - ±aA672) = 72(6). (5.8) 
We can therefore write 

72(aA6) = W(aA6) + ±(72(a)A6 + aA72(6)) - ±aA672, (5.9) 
where W(7?) is the Weyl tensor, and must satisfy 

d a -W{aAb) = 0. (5.10) 
Returning to equation (5.6) and contracting, we obtain 

<9 6 -(<9 a A72(aA6)) = <9 a A72(a) = 0, (5.11) 
which shows that the Ricci tensor 72(a) is symmetric. It follows that 

<9 Q A(i(72(a)A6 + aA72(6)) - ±aA672) = (5.12) 

and hence that 

<9 a AW(aA6) = 0. (5.13) 
Equations (5.10) and (5.13) combine to give the single equation 

<9 a W(aA6) = 0. (5.14) 

Since the d a - operation is called the 'contraction', and d a A the 'protraction', 
Hestenes & Sobczyk [11] have suggested that the sum of these be termed the 
'traction'. Equation (5.9) thus decomposes 72(73) into a 'tractionless' term 
W(7?) and a term specified solely by the matter stress-energy tensor (which 
determines 72(a) through the Einstein tensor Q{a)). There is no generally ac- 
cepted name for the part of 72(7?) that is not given by the Weyl tensor so, as it 
is entirely determined by the matter stress-energy tensor, we refer to it as the 
matter term. 
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5.2 Duality of the Weyl tensor 

To study the consequences of equation (5.14) it is useful to employ the fixed 
{7^} frame, so that equation (5.14) produces the four equations 

ffiW(ffi) + <x 2 W(er 2 ) + o- 3 W(<t 3 ) = (5.15) 

<7lW(0-l) - I<T 2 W(I(T 2 ) - I(T 3 W{I(T 3 ) = (5.16) 

-IaxWilai) + cr 2 W{a 2 ) - I<t 3 W{I(T3) = (5.17) 

-/<Tl>V(/0-l) - I(T 2 W(I(T 2 ) + (T 3 W((T 3 ) = 0. (5.18) 

Summing the final three equations, and using the first, produces 

Ia k W(Icr k ) = (5.19) 

and substituting this into each of the final three equations produces 

W(Jo- fc ) = JW(o-fc). (5.20) 

It follows that the Weyl tensor satisfies 

W(IB) = IW(B), (5.21) 

and so is 'self-dual'. This use of the term 'self-dual' differs slightly from its use 
in the 2-spinor formalism of Penrose & Rindler [51]. However, the spacetime 
algebra pseudoscalar / a similar role to the Hodge star operation (the duality 
transformation) in differential form theory, so 'self-duality' is clearly an appro- 
priate name for the relation expressed by equation (5.21). 

The fact that tractionless linear functions mapping bivectors to bivectors 
in spacetime satisfy equation (5.21) was first noted in [11]. Equation (5.21) 
means that W(B) can be analysed as a linear function on a three-dimensional 
complex space rather as a function on a real six-dimensional space. This is why 
complex formalisms, such as the Newman-Penrose formalism, are so successful 
for studying vacuum solutions. The unit imaginary employed in the Newman- 
Penrose formalism is a disguised version of the spacetime pseudoscalar [33]. 
Geometric algebra reveals the geometric origin of this 'imaginary' unit, and 
enables us to employ results from complex analysis without the need for formal 
complcxification. Furthermore, the complex structure only arises in situations 
where it is geometrically significant, instead of being formally present in all 
calculations. 

Given the self-duality of the Weyl tensor, the remaining content of equations 
(5.15)-(5.18) is summarised by the relation 

<T k W{<T k ) = 0. (5.22) 

This equation says that, viewed as a three-dimensional complex linear function, 
W(B) is symmetric and traceless. This gives W(B) five complex, or ten real 
degrees of freedom. (Since we frequently encounter combinations of the form 
scalar + pseudoscalar, we refer to these loosely as 'complex' scalars.) The gauge- 
invariant information in W(B) is held in its complex eigenvalues and, since the 
sum of these is zero, only two are independent. This leaves a set of four real 
intrinsic scalar quantities. 
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Overall, H{B) has 20 degrees of freedom, 6 of which are contained in the 
freedom to perform arbitrary local rotations. Of the remaining 14 physical 
degrees of freedom, four are contained in the two complex eigenvalues of W(B), 
and a further four in the real eigenvalues of the matter stress-energy tensor. 
The six remaining physical degrees of freedom determine the rotation between 
the frame that diagonalises Q{a) and the frame that diagonalises W{B). This 
identification of the physical degrees of freedom contained in 1Z(B) appears to 
be new, and is potentially very significant. 

5.3 The Petrov classification 

The algebraic properties of the Weyl tensor are traditionally encoded in its 
Petrov type. Here we present geometric algebra expressions for the main Petrov 
types (following the conventions of Kramer et al. [52] ) . The Petrov classification 
is based on the solutions of the eigenvalue equation 

W(B) = aB, (5.23) 

in which B is a bivector (the 'eigenbivector') and a is a complex scalar. There 
are five Petrov types: I, II, III, D and N. Type I are the most general, with two 
independent eigenvalues and three linearly-independent orthogonal eigenbivec- 
tors. Such tensors have the general form 

W(B) = \a x {B + ZF l BF 1 ) + \u 2 (B + 3F 2 BF 2 ) (5.24) 

where ot\ and a 2 are complex scalars and F\ and F 2 are orthogonal unit bivec- 
tors (F± = F 2 = 1). The eigenbivectors are Fi, F 2 and F 3 = FiF 2 , and the 
corresponding eigenvalues are (2ai — a 2 ), (2a 2 — a\) and — (ai + a 2 ). 

Type D (degenerate) are a special case of type I tensors in which two of the 
eigenvalues are the same. Physical examples are provided by the Schwarzschild 
and Kerr solutions. The region outside a spherically-symmetric source of mass 
M has 

M 

K(B)=W(B) = ^(B + 3a- r B(T r ) (5.25) 

2r 6 

where cr r is the unit radial bivector. The eigenbivectors of this function are er r , 
with eigenvalue —2M/r 3 , and any two bivectors perpendicular to cr r , with eigen- 
value M/r 3 . Similarly, TZ(B) for a stationary axisymmetric source described by 
the Kerr solution is [34, 53] 

M 

K(B)=W(B) = -- — (B + 3* r B* r ). (5.26) 

2(r — lLcosO) 6 

This differs from the radially-symmetric case (5.25) only in that its eigenvalues 
contain an imaginary term governed by the angular momentum L. Verifying 
that (5.25) and (5.26) are tractionlcss is simple, requiring only the result that, 
for an arbitrary bivector B, 

d a Fahb = d a F{ab- a-b) = -bF, (5.27) 

which employs equation (2.40) from Section 2.2. 

The fact that the Riemann tensor for the Kerr solution is obtained from 
that for the Schwarzschild solution by replacing r by r — IL cosf? is reminiscent 
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of a 'trick' used to derive the Kerr solution in the null tetrad formalism [54]. 
This is particularly suggestive given that the unit imaginary employed in the 
Newman-Penrose formalism is a disguised version of the spacetime pseudoscalar 
/. The significance of these observations is discussed further in [55, 34]. 

For tensors of Petrov type other than I, null bivectors play a significant 
role. Type II tensors have eigenvalues a\, — ai and and two independent 
cigcnbivcctors, one timclikc and one null. Type III and type JV have all three 
eigenvalues zero, and satisfy 

type III: W 3 (B) = 0, W 2 {B) ^ (5.28) 
type JV: W 2 (B) = 0. (5.29) 

An example of a type JV tensor is provided by gravitational radiation. For a 
plane-polarised gravitational wave travelling in the 73 direction 7Z(B) is given 

by 

TZ+{B) = W+(B) = i/(t-z)7 + ( 7l B 7l -7 2 S7 2 )7 + (5.30) 
for waves polarised in the direction of the 71 and 72 axes, and 

K X (B) = W X (B) = ±/(t-z)7 + ( 7l B72 + 72^71)7+ (5-31) 

for waves polarised at 45° to the axes. In both cases f(t — z) is a scalar function, 
and 7+ is the null vector 

7+ =7o +73- (5.32) 

The direct appearance of the null vector 7+ in W(B) shows that W 2 (B) = 0, 
and is physically very suggestive. Expressions of the type 7+i?7 + project the 
bivector B down the null vector 7+, and such a structure is exhibited in the 
radiation field generated by an accelerating point charge [12]. 

These examples illustrate the uniquely compact forms for the Riemann ten- 
sor afforded by geometric algebra. In terms of both clarity and physical in- 
sight these expressions are far superior to any afforded by tensor algebra, the 
Newman-Penrose formalism or differential forms. Only Wahlquist and Es- 
tabrook's (3+1) dyadic notation [56, 57] achieves expressions of comparable 
compactness, although their formalism is of limited applicability. 

5.4 The Bianchi identities 

Further information from the wedge equation (4.40) is contained in the Bianchi 
identity. One form of this follows from a simple application of the Jacobi iden- 
tity: 

[D a , [Dp, D 7 ]]A + cyclic permutations = (5.33) 
which implies that 

D a R(ep/\e~[) + cyclic permutations = 0. (5.34) 

Again, use of the D^ derivatives makes this identity straightforward, but more 
work is required to achieve a fully covariant relation. 

We start by forming the adjoint relation to (5.34), which we can write as 

d a Ad b Ad c ((a- VR(fcAc) + Q(a) x R(6Ac))S) = 0, (5.35) 
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where B is a constant bivector. We next need to establish the result that 

B v n(B 2 ) = B 2 -U(B 1 ). (5.36) 

This follows by contracting equation (5.6) with an arbitrary vector and an ar- 
bitrary bivector to obtain the equations 

d c A(a-K{cAb)) = K(aAb) (5.37) 

(B-d a )-n(aAb) = -d a B-K(aAb). (5.38) 

Protracting the second of these equations with d b and using the first, we obtain 

d b A((B-d a )-K(aAb)) = -211(B) = -d b Ad a B -K(aAb). (5.39) 

Taking the scalar product with a second bivector now gives equation (5.36). 
Using this result in equation (5.35), we now obtain 

VA(h~ 1 (^(B))) - d a Ah' 1 {R,(Sl{a)xB)) = 0. (5.40) 

Finally, acting on this equation with h and using equation (4.40), we establish 
the covariant result 

VAK(B)-d a AK(L}(a)xB) = Q. (5.41) 

This result takes a more natural form when B becomes an arbitrary function of 
position, and we write the Bianchi identity as 

d a A(a-VK(B) -K(a-VB)) =0. (5.42) 

We can extend the overdot notation of Section (2.2) in the obvious manner to 
write equation (5.42) as 

±>Ali(B) = 0, (5.43) 

which is very compact, but somewhat symbolic and hard to apply without 
unwrapping into the form of equation (5.42). 

The self-duality of the Wcyl tensor implies that 

VAW(I B) = -IV-W(B), (5.44) 

so, in situations where the matter vanishes and W(B) is the only contribution 
to 1Z(B), the Bianchi identities reduce to 

t>W{B) = 0. (5.45) 

The properties of a first-order equation such as this are discussed in more detail 
in Section 7.1. 

The contracted Bianchi identities are obtained from 

(d a Ad b )-(i>AK(aAb)) = d a ■ (K(aAV) + VK(a)) 

= 2K(V) - VK, (5.46) 

from which we can write 

Q{V) = 0. (5.47) 

An alternative form of this equation is obtained by taking the scalar product 
with an arbitrary vector, and using the symmetry of Q(a) to write 

V-g(a) = 0. (5.48) 

Written out in full, this equation takes the form 

d a ■ {L a G(b) - G(L a b) + w(a) x Q(b) - Q(w{a) x b)) = 0. (5.49) 
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5.5 Symmetries and conservation laws 

We end this section with some comments on symmetries and conservation laws. 
These comments are not all directly relevant to the applications discussed in 
this paper, but concern the general structure of GTG. 

The first significant point is that the theory is founded on an action principle 
in a 'flat' vector space. It follows that all the familiar equations relating sym- 
metries of the action to conserved quantities hold without modification. (The 
geometric algebra approach to Lagrangian field theory is developed in [31, 34].) 
Any symmetry transformation of the total action integral that is parametcriscd 
by a continuous scalar will result in a vector that is conserved with respect to 
the vector derivative V. To every such vector there corresponds a covariant 
equivalent, as is seen from the simple rearrangement 

V-J = IVA(IJ) 

= dct(h)/VA(h- 1 (i\7)) 

= det(h)V-(h( l 7)det(h)" 1 ). (5.50) 

So, if J satisfies V- J = 0, then the covariant equivalent 

J = h _1 (J)dct(h) (5.51) 

satisfies the covariant equation T> ■ J = 0. (This explains the definition of J 
in Section 3.2.) Note that if we attempt to form the canonical stress-energy 
tensor conjugate to translations we obtain the quantity Q(a) — nT(a), which the 
field equations set to zero. The overall stress-energy tensor is therefore clearly 
conserved, but this does not yield any new information. 

A second feature of our use of flat background spacetime is that all differ- 
ential equations can be recast in integral form. The integral equation form is 
not always useful, since it often forces one to deal with non-covariant quantities. 
But integral equations are particularly well-suited to handling singularities in 
the gravitational fields. Just as Gauss' theorem in elcctromagnctism can be used 
to determine the structure of an electric field source, so integral equations can 
be used to uncover the structure of the matter sources of gravitational fields. 
An example of this is provided in Section 6.4 where the Schwarzschild solution is 
shown to arise from a matter stress-energy tensor containing a single (5-function 
source of strength M. A less obvious example is contained in [55, 58], where it 
is shown that the matter generating the Kerr solution takes the form of a ring 
rotating at the speed of light, supported by a disk of tension. Such notions are 
quite different from classical general relativity. 

Killing vectors have played a significant role in the analysis of symmetries 
and conserved quantities in general relativity, and their properties are largely 
unchanged in GTG. The simplest covariant form of Killing's equation for a 
Killing vector K is that 

a-(b-VK) + b-(a-VK) = (5.52) 

for any two vector fields a and b. Contracting with d a ■ db immediately yields 
the result that K is divergenceless: 

V-K = 0. (5.53) 
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Killing vectors are frequently obtained when, in some coordinate system, the 
metric is independent of one of the coordinates. Suppose that g^ v = g^ ■ g v is 
independent of the x° coordinate. It follows that 

d 

q-q9hu = 9fi ■ (9o-T>9u) + 9u ■ {go-Vg^) = 0. (5.54) 
But, for a coordinate frame, 

g^Vg u - g v -Vg^ = h _1 (5 M e„ - <9„e M ) = (5.55) 

and using this in equation (5.54) we find that 

g fl -(g v -VK)+g I/ -(g fl -VK) = 0, (5.56) 

where K = go. Equation (5.56) is entirely equivalent to the frame-free equa- 
tion (5.52). 

A further consequence of equation (5.52) is that, for any vector a, 

a-{a-VK) = 0. (5.57) 

So, for a particle satisfying the geodesic equation v-Vv = (4.63), we see that 

d T (v-K) = v-V{vK) = K -{v-Vv) + v(v VK) = 0. (5.58) 

It follows that the quantity v ■ K is conserved along the worldline of a freely- 
falling particle. 
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Part II - - Applications 
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6 Spherically-symmetric systems 



The first full application of our formalism is to time-dependent spherically- 
symmetric fields. For simplicity, we consider only the case where the matter 
is described by a perfect fluid. The equations derived here are applicable to 
static and collapsing stars, radially-symmetric black holes and many aspects of 
cosmology, including inflation. Furthermore, in a suitable gauge the relevant 
equations are essentially Newtonian in form, making their physical interpreta- 
tion quite transparent. Applications discussed here include an analytic solution 
to the equations governing collapsing dust, and the new understanding of hori- 
zons forced by our gauge theory. This section includes an extended version of 
the work presented in [59] . 

6.1 The 'intrinsic' method 

The traditional approach to solving the gravitational field equations is to start 
with the metric g^ v , which is usually encoded as a line element 

ds 2 = g M „ dx^ dx v . (6-1) 

The analogous quantity in GTG is derived from the h-function via 

g ia , = h- 1 (e ll )-b- 1 (e v ), (6.2) 

where the {e^} comprise a coordinate frame (see also Appendix C). For a given 
matter stress-energy tensor, the field equations then yield a set of non-linear, 
second-order differential equations for the terms in g^ v . These equations are 
notoriously hard to solve. On the other hand, any metric is potentially a solu- 
tion to the Einstein equations — one where the matter stress-energy tensor is 
determined by the corresponding Einstein tensor. This approach, in which the 
tail wags the dog, has recently probably been more popular! Here we present a 
new approach to solving the gravitational field equations. The method is closely 
tied to our gauge-theoretic understanding of gravity, but can always be used to 
generate a metric which solves the Einstein equations. So, even if one rejects 
our gauge-theory description of gravity, the techniques developed below can still 
be viewed as providing a new method for studying the Einstein equations. 
Under a local rotation the vector h~ 1 (a) transforms as 

h _1 (a) i?h" 1 (a) J R. (6.3) 

It follows that the metric (6.2) is invariant under rotation-gauge transforma- 
tions. This is in keeping with our earlier observation that, at the classical level, 
it is possible to work with a set of equations that are invariant under rotation- 
gauge transformations, and this is precisely what general relativity docs. This 
approach has the advantage of removing a number of degrees of freedom from 
the theory, but one pays a heavy price for this: the equations become second- 
order and one has to deal with complicated non-linear terms. The approach we 
develop here is quite different. We keep the rotation-gauge field explicit, and 
work entirely with quantities that transform covariantly under position-gauge 
transformations. Such quantities include h(V), u>(a) and 71(B). We therefore 
work with directional derivatives of the form L a = a-h(V), and treat u)(a) as 



47 



an arbitrary field. The relationship between h(a) and ui(a) is then encoded in 
the commutation relations of the L a . This setup is achieved by initially making 
a suitably general ansatz for the h-function. This trial form is then substi- 
tuted into equation (4.53) to find the general form of Lu(a). An arbitrary u)(a) 
field consistent with this general form is then introduced, resulting in a set of 
equations relating commutators of the L a derivatives to the variables in u>(a). 

Next, the Riemann tensor 1Z(B) is constructed in terms of abstract first- 
order derivatives of the oj(a) and additional quadratic terms. The rotation- 
gauge freedom is then removed by specifying the precise form that 1Z(B) takes. 
For example, one can arrange that W(B) is diagonal in a suitable frame. This 
gauge fixing is crucial in order to arrive at a set of equations that are not 
under-constrained. With TZ(B) suitably fixed, one arrives at a set of relations 
between first-order abstract derivatives of the u>(a), quadratic terms in io(a), 
and matter terms. The final step is to impose the Bianchi identities, which 
ensure overall consistency of the equations with the bracket structure. Once all 
this is achieved, one arrives at a fully 'intrinsic' set of equations. Solving these 
equations usually involves searching for natural 'integrating factors'. These in- 
tegrating factors provide 'intrinsic' coordinates, and many of the fields can be 
expressed as functions of these coordinates alone. The final step is to 'coordi- 
natise' the solution by making an explicit (gauge) choice of the h-function. The 
natural way to do this is to ensure that the coordinates used in parameterising 
h(a) match the intrinsic coordinates defined by the integrating factors. 

The method outlined above is quite general and can be applied to a wide 
range of problems [34, 60, 61]. Here we employ it in the analysis of time- 
dependent spherically-symmetric systems. 

6.2 The intrinsic field equations 

We start by introducing a set of spherical polar coordinates 
fixed {7^} frame we define: 

t = x-'fo cos0 = x-j 3 /r 

r = y (.TA70) 2 tan0 = (.t-7 2 )/(x-7 1 ). 
The associated coordinate frame is 
e* = 7o 

e r = a;A7o7o/r 

eg = rcos#(cos07! + sin</>7 2 ) — rsin0 7 3 
= r sin(9(— sin07i +cos</>72) 

and the dual-frame vectors are denoted by {e*, e r , e , e^}. We will also frequently 
employ the unit vectors 6 and 4> defined by 

9 = eg/r, (f) = e^/(r sin#). (6.6) 

Associated with these unit vectors are the unit timelike bivectors 

cr r = e r e u crg=9e t , cr^ = 4>e t , (6.7) 

which satisfy 

(TrfTgcr^ = ete r 9cj) = I- (6-8) 



In terms of the 
(6.4) 
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The dual spatial bivectors are given by 



Icr r = —9(f), lag = e r 4>, Icr^ = —e r 9. (6-9) 
Throughout we use the abbreviations 

*=Jr (6 - io) 

The h-function 

Our first step towards a solution is to decide on a general form of the h-function 
that is consistent with spherical symmetry. Suppose that B is a constant spatial 
bivector (e t -B = 0), and define 

R = e B ' 2 (6.11) 
x' = RxR. (6.12) 

Then, in analogy with electromagnetism, the gravitational fields will be spheri- 
cally symmetric if rotating h(et) to Rh(a)R and displacing it to the back-rotated 
position x' leaves h(a) unchanged. Hence rotational symmetry is enforced 
through the requirement that 

Rh x ,(RaR)R = h(a). (6.13) 

This symmetry immediately implies that the {e r , e'} and {e , e^} pairs decouple 
from each other, and the action of h(a) on the 9 and <fi vectors is further restricted 
to the form 

h(0) = u9 + B6 

. (6.14) 
b(<j>) = a<j>-/39. 

However, the skew-symmetric term parameterised by (3 can always be removed 
by a rotation in the Icr r plane, so we can assume that h(a) is diagonal on 
{e e ,e^}. No such assumption can be made for the {e r ,e*} vectors, so we take 
h(a) as having the general form 

h(e 4 ) = he} + f 2 e r 
h(e r ) - 9l e r + g 2 e t 

He«)=ae° ^ 
h(e' l> ) = ae*, 

where /i, f 2 , <?i, <?2 and a are all functions of t and r only. We retain the gauge 
freedom to perform a boost in the cr r direction, and this freedom is employed 
later to simplify the equations. Our remaining position-gauge freedom lies in 
the freedom to reparameterise t and r, which does not affect the general form 
of (6.15). A natural parametcrisation will emerge once the 'intrinsic' variables 
have been identified. 

The u;-function 

To find a general form ui(a) consistent with (6.15) we substitute equation (6.15) 
into equation (4.53) for uj(a) as a function of h(a). Where the coefficients contain 
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<T8 


<T0 


e t -V 







GI*t 


-GIa e 


e r -T> 







FI<T4, 


-Flag 


9-V 


T(T0 


- S ICTj, 


— T(T r 


S la r 


4>-v 


Ter^ 


+ SI(Tg 


— SIcr r 


-T<T r 



Tabic 3: Covariant derivatives of the polar-frame unit timelike bivectors. 



derivatives of terms from h(a) new symbols are introduced. Undifferentiated 
terms from h(a) appearing in w(a) are left in explicitly. These arise from frame 
derivatives and the algebra is usually simpler if they are included directly This 
procedure results in the following form for uj(a): 

uj(e t ) = Ge r e t 
oj(e r ) = Fe r e t 

u){6) = S6e t + (T-a/r)e r 

uo{(t>) = S<f>e t + (T - a/r)e r 4>, 

where G, F, S and T are functions of t and r only. The important feature of 
these functions is that they are position-gauge covariant. 

Substituting this definition for u(a) into equations (4.45) and (4.46) we find 
that the bracket relations are as follows: 

[Lt, L r ] = GL t — FL r [L r , Lg] = —TLq 

[Lt, L S ] = -SL e [L r , L $ ) = -TL $ ( 6 . 17 ) 

[L t ,L i ] = -SL i [L e ,L^}=0, 



where 



L t = e r h(V) L § = 0-h(V) 

L r = e r -h(V) L,=fh(V). 



(6.18) 



The use of unit vectors in these derivatives eliminates the need to calculate 
irrelevant coordinate derivatives. A set of bracket relations such as (6.17) is 
precisely what one aims to achieve — all reference to the h-function is removed 
and one deals entirely with position-gauge covariant quantities. 



The Riemann tensor 

Having found a suitable form for ui{a) we next use equation (4.48) to calculate 
1Z(B). This derivation is simplified by judicious use of the results listed in 
Table 3. The only subtlety in the derivation is the removal of terms involving 
derivatives of a/r using the bracket relations (6.17). Since a/r = Lg8 we have 

L t (a/r) = L t L § 6 = [Lt, L § }9 = -Sa/r (6.19) 

and 

L r (a/r) = L r L § = [L r , L § }9 = -Ta/r. (6.20) 
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Application of equation (4.48) is now straightforward, and leads to the Riemann 
tensor 



K(a r ) 


= (L r G - 


L t F 


+- G 2 - F 2 ' 


)<T r 






n{* 6 ) 


= (-L t S 


+ GT 


- S 2 )a e + 


(L t T H 


-ST~ 


SG)Itr, 




= (-L t S 


+ GT 


- S 2 )<r^ - 


(L t T- 


h5T- 


SG)I<r 


11(1(7$) 


= (L r T + 


T 2 - 


FS)I<rt - 


{L r S^ 


-ST- 


FT)a e 


n(Icrg) 


= (L r T + 


T 2 - 


FS)I<t 9 + 


(L r S-i 


-ST- 


FT)a^ 


K(I(Tr) 


= (S 2 + 




[a/rf)Ia 


r • 







The matter field and gauge fixing 

Now that we have found 1Z(B) in terms of 'intrinsic' functions and their first 
derivatives, we must next decide on the form of matter stress-energy tensor that 
the gravitational fields couple to. We assume that the matter is modelled by an 
ideal fluid so we can write 

T(a) = (p + p)a-vv — pa, (6.22) 

where p is the energy density, p is the pressure and v is the covariant fluid 
velocity (v 2 = 1). Radial symmetry means that v can only lie in the et and e r 
directions, so v must take the form 

v = coshit et + sinhu e r . (6.23) 

But, in restricting the h-function to the form of equation (6.15), we retained 
the gauge freedom to perform arbitrary radial boosts. This freedom can now 
be employed to set v = e t , so that the matter stress-energy tensor becomes 

T(a) = (p + p)a-e t e t - pa. (6.24) 

There is no physical content in the choice v = et as all physical relations must 
be independent of gauge choices. In particular, setting v = e t does not mean 
that the fluid is 'at rest', or that we are 'comoving with the fluid'. An observer 
comoving with the fluid will have a covariant velocity et, but this implies no 
special relationship with the time coordinate t, since the observer's trajectory 
would have d\x = h(et) and nothing has yet been said about the specific form 
of h(a). 

In setting v = e t all rotation-gauge freedom has finally been removed. This 
is an essential step since all non-physical degrees of freedom must be removed 
before one can achieve a complete set of physical equations. Note that the 
rotation gauge has been fixed by imposing a suitable form for 1Z(B), rather 
than restricting the form of h(a). The reason for working in this manner is 
obvious — TZ-{B) deals directly with physically-measurable quantities, whereas 
the algebraic structure of the h-function is of little direct physical relevance. 

From equation (5.9) the source term in 1Z(B) is given by 

K(aAb) - W(aAb) = An(aAT(b) + T(a)Ab - ^TaAb) (6.25) 

where T = d a -T(a) is the trace of the matter stress-energy tensor. With T(a) 
given by equation (6.24), IZ(B) is restricted to the form 

K(B)=W(B)+^((p + p)B-e t et-lpB). (6.26) 
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Comparing this with equation (6.21) we find that W(B) has the general form 
W(ov) = <X\(T r W(I(T r ) = a±Icr r 

W(<r fl ) = a 2 o- e + W(/er e ) = a 3 /<x + /? 2 e>> (6.27) 
W(a>) = a 2 o> - /3i/cr e >V(/(T0) = a 3 Itr^ - /3 2 <x e . 

But W(-B) must be self-dual, so a.\ = a±, a 2 = a 3 and /3i = — /3 2 . It must also 
be symmetric, which implies that /3i = /3 2 . It follows that /3i = [3 2 = 0. Finally, 
W(i?) must be traceless, which requires that a>i + 2a 2 = 0. Taken together, 
these conditions reduce W(B) to the form 

W(B) = -j(B + 3<r r Bcr r ), (6.28) 

which is of Petrov type D. It follows from the form of TZ(I<r r ) that if we set 

A= \(-S 2 + T 2 -(a/r) 2 ) (6.29) 

then the full Riemann tensor must take the form 

K{B) = (A + lnp){B + 3(T r B(T r ) + An((p + p)B-e t e t - \pB). (6.30) 

Comparing this with equation (6.21) yields the following set of equations: 

L t S = 2A + GT - S 2 - 4ttp (6.31) 

L t T = S(G — T) (6.32) 

L r S = T(F-S) (6.33) 

L r T = —2A + FS — T 2 - Airp (6.34) 

L r G - L t F = F 2 -G 2 + AA + Att( p + p). (6.35) 

The Bianchi identity 

We are now close to our goal of a complete set of intrinsic equations. The 
remaining step is to enforce the Bianchi identities. The contracted Bianchi 
identity (5.47) for a perfect fluid results in the pair of equations 

V-(pv) +pV-v = (6.36) 
(p + p)(v-Vv)Av- (Dp) Av = 0. (6.37) 

Since (v-T> v)/\v is the acceleration bivector, the second of these equations relates 
the acceleration to the pressure gradient. For the case of radially-symmetric 
fields, equations (6.36) and (6.37) reduce to 

L tP = -(F + 2S)(p + P ) (6.38) 
L rP = -G(p + p), (6.39) 

the latter of which identifies G as the radial acceleration. The full Bianchi iden- 
tities now turn out to be satisfied as a consequence of the contracted identities 
and the bracket relation 

[L t , L r ] = GLt — FL r . (6.40) 
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Equations (6.19), (6.20), (6.31)-(6.35), the contracted identities (6.38) and 
(6.39), and the bracket condition (6.40) now form the complete set of intrinsic 
equations. The structure is closed, in that it is easily verified that the bracket 
relation (6.40) is consistent with the known derivatives. The derivation of such 
a set of equations is the basic aim of our 'intrinsic method'. The equations deal 
solely with objects that transform covariantly under displacements, and many 
of these quantities have direct physical significance. 

Integrating factors 

To simplify our equations we start by forming the derivatives of A. From equa- 
tions (6.19), (6.20) and (6.31)-(6.35) it follows that 

L t A + 3SA = 2irSp (6.41) 
L r A + 3TA = -2nTp. (6.42) 

These results, and equations (6.32) and (6.33), suggest that we should look for 
an integrating factor for the L t + S and L r + T operators. Such a function, X 
say, should have the properties that 

L t X = SX, L r X = TX. (6.43) 

A function with these properties can exist only if the derivatives are consistent 
with the bracket relation (6.40). This is checked by forming 

[L u L r ]X = L t (TX) - L r (SX) 
= X(L t T - L r S) 
= X{SG-FT) 

= GL t X - FL r X, (6.44) 

which confirms that the properties of X are consistent with (6.40). Estab- 
lishing the existence of integrating factors in this manner is a key step in our 
method, because the integrating factors play the role of intrinsically defined co- 
ordinates. If the h-function is parameterised directly in terms of these functions, 
the physical status of the quantities in it becomes clearer. In the present case, 
equations (6.19) and (6.20) show that r/a already has the properties required 
of X, so it is r/a that emerges as the intrinsic distance scale. It is therefore 
sensible that the position-gauge freedom in the choice of r should be absorbed 
by setting a = 1. This then sets the intrinsic distance scale equal to r, lifting 
r from the status of an arbitrary coordinate to that of a physically-measurable 
quantity. 

Having fixed the radial scale with the position-gauge choice 

r = X, a = l, (6.45) 

we can make some further simplifications. From the form of h(a) (6.15) and 
equations (6.43) and (6.45) we see that 

gi = L r r = Tr (6.46) 
g 2 = L t r = Sr, (6.47) 
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which gives two of the functions in h(et). We also define 

M = -2r 3 A = \r{g 2 2 - g x 2 + 1), (6.48) 

which satisfies 

L t M = -^Tir 2 g 2 p (6.49) 

and 

L r M = 47rr 2 5 ip. (6.50) 
The latter shows that M plays the role of an intrinsic mass. 

The 'Newtonian' gauge 

So far a natural distance scale has been identified, but no natural time coordinate 
has emerged. To complete the solution it is necessary to make a choice for the 
t coordinate, so we now look for additional criteria to motivate this choice. We 
are currently free to perform an arbitrary r and i-dependent displacement along 
the et direction. This gives us complete freedom in the choice of f 2 function. An 
indication of how this choice should be made is obtained from equations (6.49) 
and (6.50) for the derivatives of M (6.49), which invert to yield 

dM _ ~^ 9l g 2 r 2 {p + p) 



dt figi - hgi 

dM = inr 2 (f igi p + f 2 g 2 p) 
dr fxgi ~ /2ff2 



(6.52) 



The second equation reduces to a simple classical relation if we choose f 2 = 0, 
as we then obtain 

d r M = Airr 2 p, (6.53) 

which says that M(r, t) is determined by the amount of mass-energy in a sphere 
of radius r. There are other reasons for choosing the time variable such that 
f 2 = 0. For example, we can then use the bracket structure to solve for f\. 
With f 2 = we have 

L t = fid t + g 2 d r (6.54) 
L r = gid r , (6.55) 

and the bracket relation (6.40) implies that 

L r fi = -Gfi. (6.56) 

It follows that 

h = <t) cxp (- J G/ 9l dr^j . (6.57) 

The function e{t) can be absorbed by a further f-dependent rescaling along et 
(which does not change f 2 ), so with f 2 = we can reduce to a system in which 



fi = cxp I - J G/ 9l dr) . (6.58) 

Another reason why f 2 — is a natural gauge choice is seen when the pressure 
is zero. In this case equation (6.39) forces G to be zero, and equation (6.58) 
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then forces fi = 
fluid) then has 



1. A free-falling particle with v = et (i.e. comoving with the 
ie t + fe r — e t + g 2 e r , (6.59) 



where henceforth the dots denote differentiation with respect to the proper 
time. Since 1 = 1 the time coordinate t matches the proper time of all observers 
comoving with the fluid. So, in the absence of pressure, we are able to recover a 
global 'Newtonian' time on which all observers can agree (provided all clocks are 
correlated initially). Furthermore it is also clear from (6.59) that 172 represents 
the velocity of the particle. Hence equation (6.51), which reduces to 



in the absence of pressure, has a simple Newtonian interpretation — it equates 
the work with the rate of flow of energy density. Equation (6.48), written in the 
form 



is also now familiar from Newtonian physics — it is a Bernoulli equation for 
zero pressure and total (non-relativistic) energy (<?i 2 — l)/2. 

For these reasons we refer to f 2 = as defining the 'Newtonian' gauge. The 
applications discussed in the following sections vindicate our claim that this 
is the natural gauge for radially-symmetric systems. The full set of equations 
in the Newtonian gauge are summarised in Table 4. They underlie a wide 
range of phenomena in relativistic astrophysics and cosmology. The closest 
analogue of the Newtonian gauge description of a spherically-symmetric system 
is provided by Gautrcau's 'curvature coordinates' [62] (see also [63, 64]). This 
description employs a set of geodesic clocks in radial freefall, comoving with 
the fluid. However, such a description can only be applied if the pressure is 
independent of radius, whereas the Newtonian gauge description is quite general. 

One aspect of the equations in Table 4 is immediately apparent. Given an 
equation of state p = p(p), and initial data in the form of the density p(r,to) 
and the velocity g 2 (r, to), the future evolution of the system is fully determined. 
This is because p determines p and M on a time slice, and the definition of M 
then determines g\. The equations for L r p 1 L r g\ and L r g 2 then determine the 
remaining information on the time slice. Finally, the L t M and L t g 2 equations 
can be used to update the information to the next time slice, and the process 
can then start again. The equations can thus be implemented numerically as 
a simple set of first-order update equations. This fact considerably simplifies 
the study of collapsing matter, and should be particularly significant in cur- 
rent studies of the critical phenomena associated with horizon and singularity 
formation [65, 66]. 

6.3 Static matter distributions 

As a simple first application we consider a static, radially-symmetric matter 
distribution. In this case p and p are functions of r only. Since M(r, t) is now 



d t M = -4nr 2 g 2 p 



(6.60) 



given by 




(6.62) 
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The h-function 



h(e*) 


= /ie* 


h(e r ) 




He 6 ) 


= e e 


h(e*) 





The w function 



w(et) 


= Ge r et 




w(e r ) 


= -Fe r e t 






= 92/rOet 4 


- (51 - l)/re r 6 




= 92/r^et - 


\- (91 - l)Ae r 



Directional derivatives 



it = ji<9 t + ,g 2 9 r 
i r = g\d r 



G, F and h 



L t gi = G52, = Fgi 

h = cxp(/ r -G/gidr) 



Definition of M 



M = \r{g 2 2 - 9l 2 + 1) 



Remaining derivatives 



L t g2 = Ggi - 


- M/r 2 


— Airrp 


L r g\ = Fg 2 - 


f- M/r 2 


— Airrp 



Matter derivatives 



L t M = -4:irr 2 g 2 p, L r M = \-nr 2 gxp 
L t p=-(2g 2 /r + F){p + p) 
L r p = -G{p + p) 



Ricmann tensor 



11(B) =4n((p + p)B-e t et-2p/3B) 
-±(M/r 3 - 4tt / 9/3)(B + 3o- r Bcr r 



Stress-energy tensor T{a) = (p + p)a-e t e-t — pa 



Table 4: Equations governing a radially-symmetric perfect fluid. An equation 
of state and initial data p(r, to) and g 2 (r, to) determine the future evolution of 
the system. 



it follows that 



L t M = 4nr g 2 p = —Anr g 2 p. 



(6.63) 
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For any physical matter distribution p and p must both be positive, in which 
case equation (6.63) can be satisfied only if 

g 2 = F = 0. (6.64) 

Since g 2 = we see that g\ is given simply in terms of M(r) by 

.9i 2 = 1 - 2M{r)/r, (6.65) 

which recovers contact with the standard line element for a static, radially- 
symmetric field. It is immediately clear that a solution exists only if 2M(r) < r 
for all r. This is equivalent to the condition that a horizon has not formed. 
The remaining equation of use is that for L t g 2 , which now gives 

Ggi = M(r) /r 2 + Anrp. (6.66) 

Equations (6.65) and (6.66) combine with that for L r p to give the famous 
Oppenhcimcr-Volkov equation 

dp _ {p + p){M{r)+A^p) 

dr r(r-2M(r)) ' ( j 

At this point we have successfully recovered all the usual equations governing 
a non-rotating star, and the description is therefore unchanged from that of 
standard general relativity. The work involved in recovering these equations 
from the full time-dependent case is minimal, and the final form of h(a) is very 
simple (it is a diagonal function). Furthermore, the meaning of the t and r 
coordinates is clear, since they have been defined operationally. 

The solution extends straightforwardly to the region outside the star. We 
now have M constant, and 

fi = 1/gi = (1 - 2Af/r)" 1 / 2 , (6.68) 

which recovers the Schwarzschild line element. It follows that all predictions for 
the behaviour of matter in the star's gravitational field, including those for the 
bending of light and the perihelion precession of Mercury, are unchanged. 

6.4 Point source solutions — black holes 

The next solution of interest is obtained when the matter is concentrated at a 
single point (r = 0). For such a solution, p — p — everywhere away from the 
source, and the matter equations reduce to 



L t M = 
L r M = 



M = constant. (6.69) 



Retaining the symbol M for this constant we find that the equations reduce to 

L t gi = Gg 2 (6.70) 

L r g 2 = Fg x (6.71) 

and 

91 -g-2 2 = l-2M/r. (6.72) 
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No further equations yield new information, so we have an under-determined 
system of equations and some additional gauge fixing is needed to choose an 
explicit form of h(a). The reason for this is that in the vacuum region the 
Riemann tensor reduces to 

M 

K(B) = -jj (B + 3a r B(T r ). (6.73) 

This tensor is now invariant under boosts in the cr r plane, whereas previously the 
presence of the fluid velocity in the Riemann tensor vector broke this symmetry. 
The appearance of this new symmetry in the matter-free case manifests itself 
as a new freedom in the choice of h-function. 

Given this new freedom, we should look for a choice of g\ and 52 that sim- 
plifies the equations. If we attempt to reproduce the Schwarzschild solution we 
have to set gi = 0, but then we immediately run into difficulties with g\, which 
is not defined for r < 2M. We must therefore look for an alternative gauge 
choice. We show in the following section that, when p = 0, g\ controls the en- 
ergy of infalling matter, with particles starting at rest at r = 00 corresponding 
to gi = 1. A sensible gauge choice is therefore to set 



5i = 1 (6.74) 



so that 



.92 



V2M/r (6.75) 
G = (6.76) 
F = -M/(g 2 r 2 ) (6.77) 



and 



fi = I- (6.78) 
In this gauge the h-function takes the remarkably simple form 



h(a) = a- sp2Mjra-e r e u (6.79) 

which only differs from the identity through a single term. From the results of 
Section 2.3 the extension to the action of h on an arbitrary multivector A is 

straightforward: 

h(A) = A - x /2M/r(A-e r )Ae t . (6.80) 

It follows that det(h) = 1 and the inverse of the adjoint function, as defined 
by (2.53), is given by 



h-\A) = A + y/2M/r(A-e t )Ae r . (6.81) 



Point-particle trajectories 

To study the properties of the solution (6.79) wc consider the equation of motion 
for infalling matter. For a particle following the trajectory x(r), with r the 
proper time, we have 

v = ie t + (i^2M/r + r)e r + 0e g + e</> , (6.82) 
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where dots denote the derivative with respect to r. Since the h-function is 
independent of t we have, from equation (4.65), 

h -1 (et)-y = (1 - 2M/r)i-r^2M/r = constant, (6.83) 

and, for particles moving forwards in time (t > for r — ► oo), we can write 

(1 - 2M/r)i = a + fy/2M/r, (6.84) 

where the constant a satisfies a > 0. The f equation is found from the constraint 
that v 2 = 1, which gives 

r 2 = a 2 - (1 -2M/r)(l + r 2 (0 2 +sin 2 6></> 2 )). (6.85) 

The horizon lies at r = 2M since, for r < 2M, the velocity f must be negative. 
It might appear that an attempt to integrate equation (6.84) will run into diffi- 
culties with the pole at horizon, but this not the case. At r = 2M we find that 
f = —a and this cancels the pole. All particles therefore cross the horizon and 
reach the singularity in a finite coordinate time. 

Specialising to the case of radial infall, we see from equation (6.85) that the 
constant a 2 — 1 can be identified with twice the particle's initial energy (for 
a unit mass particle). Furthermore, equation (6.85) shows immediately that f 
satisfies 

r=-M/r 2 , (6.86) 

a feature of motion in a spherically-symmetric gravitational field that is rarely 
emphasised in general relativity texts. Some possible matter and photon tra- 
jectories are illustrated in Figure 1. In the case where the particle is dropped 
from rest at r = oo equations (6.84) and (6.85) reduce to 

f = -y/2M/r, t=l, (6.87) 

and we recover an entirely Newtonian description of the motion. The properties 
of a black hole are so simple in the gauge defined by (6.79) that it is astonish- 
ing that this gauge is almost never seen in the literature (see [67, 64] for some 
exceptions). Presumably, this is because the line element associated with (6.79) 
does not look as natural as the h-function itself and hides the underlying sim- 
plicity of the system. Part of the reason for this is that the line element is not 
diagonal, and relativists usually prefer to find a coordinate system which diag- 
onaliscs g^ v . Even when the frecfall time coordinate t is employed, a different 
radial coordinate is usually found to keep the metric diagonal [68, 69]. 

Since the gauge defined by g\ = 1 and g 2 = —y/2M/r extends our aim of 
keeping the equations in a simple Newtonian form, we refer to this solution 
as defining the 'Newtonian gauge' vacuum solution. We show in Section 6.5 
that this gauge arises naturally from the description of collapsing dust. In 
the Newtonian gauge one hardly needs to modify classical reasoning at all to 
understand the processes involved — all particles just cross the horizon and 
fall into the singularity in a finite coordinate time. And the horizon is located 
at r = 2M precisely because we can apply Newtonian arguments! The only 
departures from Newtonian physics lie in relativistic corrections to the proper 
time taken for infall, and in modifications to the equations for angular motion 
which lead to the familiar results for orbital precession. 
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Figure 1: Matter and photon trajectories for radial motion in the in the New- 
tonian gauge. The solid lines are photon trajectories, and the horizon lies at 
r=2. The broken lines represent possible trajectories for infalling matter. Tra- 
jectory I is for a particle released from rest at r = 4. Trajectory II is for a 
particle released from rest at r = oo. 

When extracting physical predictions in the Newtonian, or any other, gauge, 
it is important to ensure that the problem is posed in a gauge-invariant manner. 
For example, one can envisage a simple experiment with two observers initially 
at rest outside a black hole at a distance ro, where this distance is defined 
in terms of the magnitude of 7Z(B). One observer can then start free-falling, 
and agree to emit photons of a chosen frequency at regular intervals. If one 
then computes what the remaining, stationary observer sees as a function of 
their proper time, this is clearly something physically meaningful. It is not 
hard to show that the predictions for this are gauge invariant, as they must 
be. Furthermore, if everything takes place outside the horizon, one can work 
in the 'Schwarzschild' gauge with g 2 = 0. But, to see what happens as the 
free-falling observer crosses the horizon, a global solution such as (6.79) must 
be used. One still finds that the signal from the free-falling observer becomes 
successively more red-shifted and less intense, as predicted when working with 
the Schwarzschild metric, but the free- falling observer crosses the horizon in a 
finite coordinate time. 

Horizons and time-reversal asymmetry 

Our picture of the gravitational fields due to a radially-symmetric point source 
is rather different from that of general relativity. These differences arc seen most 
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clearly in the effects of time reversal. Time reversal is achieved by combining 
the displacement 

f(x) = —e t xe t = x (6.88) 

with the reflection 

h» = -e t h(a)e t , (6.89) 
resulting in the the time-reversed solution 

h*(a) = e t b x ,(e t ae t )e t . (6.90) 

As an example, the identity function h(a) = a is time- reverse symmetric — as it 
should be. The displacement (6.88) is a gauge transformation and cannot have 
any physical consequences. The reflection (6.89) is not a gauge transformation, 
however, and can be used to transform between physically distinct gauge sectors. 
The reflection (6.89) is lost when the metric is formed, so general relativity 
cannot handle time-reversal in the same manner as GTG. 

With the h-function described by equation (6.15), and with the {fi} and 
{gi} functions of r only, the effect of (6.90) is simply to change the sign of the 
off-diagonal elements f 2 and g 2 . For example, applied to the solution (6.79), 
the transformation (6.90) produces the time-reversed solution 

h*(a) = a + y/2M/ra-e r e t . (6.91) 

The result is a solution in which particles inside the horizon are swept out. Once 
outside, the force on a particle is still attractive but particles cannot re-enter 
through the horizon. 

This lack of time-reversal symmetry is not a feature of the various gauge 
choices made in arriving at (6.79); it is an inevitable result of the presence of 
a horizon. To see why, we return to the equations in the form prior to the 
restriction to the Newtonian gauge. The h-field is as defined by equation (6.15) 
with a = 1 and the {fi} and {gi} functions of r only. The general set of 
time-independent vacuum equations still have M constant and 

gi - .92 2 = 1 - 2M/r (6.92) 



with 

d r gi = G, 

The bracket relation (6.40) now gives 

from which it follows that 

9 r {fm - hdl) 



d r g 2 = F. (6.93) 

(6.94) 

= <9 r det(h) = 0. (6.95) 



gid r h - gxd r h = Gfi - Ff 2 



Hence det(h) is a constant, with its value dependent on the choice of position 
gauge. Since TZ{B) tends to zero at large r, we can always choose the gauge 
such that h(a) tends to the identity as r — > oo. In this case det(h) must be one, 
so we can write 

fxQl - / 2 .92 = 1. (6.96) 
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If we form the line element derived from our general h-function we obtain 



ds 2 = (1 - 2M/r) dt 2 + 2(/ifla - f 2 gi)dtdr - (/i 2 - f 2 2 ) dr 2 
~r 2 {d8 2 + sin 2 6> <# 2 ). 



(6.97) 



The off-diagonal term here is the one that breaks time-reversal symmetry. But 
we must have g± = ±g 2 at the horizon, and we know that figi — f 2 g 2 = 1 
globally. It follows that 



so the line element (6.97) must break time reversal symmetry at the horizon [43]. 
In fact, the asymmetry is even more pronounced. Once inside the horizon, 
equation (6.92) forces a non-zero g 2 , so the h-function cannot be time-reverse 
symmetric anywhere inside the horizon. This link between the existence of a 
horizon and the onset of time-reversal asymmetry is one of the most satisfying 
aspects of GTG. Furthermore, the requirement that a sign be chosen for fig 2 — 
f 2 gi at the horizon shows that a black hole has more memory about its formation 
than simply its mass M — it also remembers that it was formed in a particular 
time direction. We will see an example of this in the following section. 

At the level of the metric the discussion of time-reversal is much less clear. 
For example, inside the horizon a valid h-function is obtained by setting /i and 
<7i to zero. Since f 2 and g 2 are non-zero, this h-function is manifestly not time- 
reverse symmetric. However, the line element generated by this h-function is 
just the Schwarzschild line clement 



ds 2 = (1 - 2M/r) dt 2 - (1 - 2A//?-)" 1 dr 2 - r 2 (d6 2 + sin 2 6d<j) 2 ), (6.99) 



which is usually thought of as being time-reverse symmetric. Clearly, our gauge 
theory probes questions related to time-reversal symmetry at a deeper level than 
general relativity. The consequences of our new understanding of time-reversal 
will be met again in Section 8, where we study the Dirac equation in a black-hole 
background. 

In general relativity one of the most important results for studying radially- 
symmetric fields is Birkhoff's theorem. This can be stated in various ways, 
though the most usual statement is that the line element outside a radially- 
symmetric body can always be brought to the form of (6.99). As we have 
seen, this statement of Birkhoff's theorem is correct in GTG only if no hori- 
zon is present. However, the more general statement, that the fields outside a 
spherically-symmetric source can always be made stationary, does remain valid. 

The Kruskal extension and geodesic completeness 

In modern general relativity, the line element (6.99) is not viewed as representing 
the final form of the metric for a radially-symmetric black hole. The full 'maxi- 
mal' solution was obtained by Kruskal [70] , who employed a series of coordinate 
transformations that mixed advanced and retarded Eddington-Finkelstein co- 
ordinates. The Kruskal extension describes a spacetime that contains horizons 
and is time-reverse symmetric, so can have no counterpart in GTG. Further- 
more, the Kruskal solution has two distinct regions for each value of r [16] and 
so is, topologically, quite distinct from the solutions admitted in GTG. This is 



fig 2 - hgi = ±1 at r = 2M, 



(6.98) 
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because any solution of our equations must consist of h(a) expressed as a func- 
tion of position x. The form of this position dependence is arbitrary, but it must 
be present. So, when the coordinate r is employed in defining the h-function, 
this always represents a particular function of the vector x. The point r = is, 
by definition, a single point in space (an unbroken line in spacetime). No fields 
can alter this fact. The Kruskal solution contains two separate regions with the 
label r = 0, so immediately fails in GTG. Instead of the full Kruskal extension 
with 4 sectors (usually denoted I, II, I' and II' [16]), GTG admits two distinct 
solutions, one containing the sectors I and II, and the other containing I' and 
II . These solutions are related by the discrete operation of time reversal, which 
is not a gauge transformation. This splitting of a single time-reverse symmet- 
ric solution into two asymmetric solutions is typical of the transition from a 
second-order to a first-order theory. Similar comments apply to the maximal 
extensions of the Reissner-Nordstrom and Kerr solutions. The infinite chain 
of 'universes' general relativity admits as solutions have no counterpart in our 
theory. Instead, we use integral equations to determine the nature of the matter 
singularities, precisely as one would do in electromagnctism [55] 

In general relativity the Kruskal solution is the unique maximal continuation 
of the Schwarzschild metric. The fact that it has no analogue in GTG means 
that our allowed solutions are not 'maximal' and forces us to address the issue 
of geodesic incompleteness. For the solution (6.79) geodesies exist that cannot 
be extended into the past for all values of their affine parameter. But we have 
already seen that the presence of a horizon commits us to a choice of time 
direction, and in the following section we show how this choice is fixed by the 
collapse process. So, if we adopt the view that black holes arise solely as the 
endpoint of a collapse process, then there must have been a time before which 
the horizon did not exist. All geodesies from the past must therefore have 
come from a period before the horizon formed, so there is no question of the 
geodesies being incomplete. We therefore arrive at a consistent picture in which 
black holes represent the endpoint of a collapse process and the formation of 
the horizon captures information about the direction of time in which collapse 
occurred. This picture is in stark contrast with general relativity, which admits 
eternal, time-reverse symmetric black-hole solutions. 

Coordinate transformations and displacements 

The coordinate transformations employed in general relativity have two distinct 
counterparts in GTG: as passive re-labcllings of the coordinates employed in a 
solution, such as changes of variables used for solving differential equations; and 
as disguised forms of position-gauge transformations. An example of the lat- 
ter is the transformation between the Schwarzschild and advanced Eddington- 
Finkelstein forms of the spherically-symmetric line element. This is achieved 
with the coordinate transformations 



t'-r' = t-(r + 2Mln(r - 2M)) 



(6.100) 
(6.101) 



r = r, 



which can be viewed as the result of the displacement defined by 



f(x) =x' = x- 2Mln(r - 2M)e t . 



(6.102) 
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This displacement is to be applied to the solution 



h(o) = A~ 1 / 2 a-e t e* + A^ 2 a-e r e r + aAa r o> (6.103) 

where 

A = l-2M/r. (6.104) 

Clearly the gravitational fields are defined only outside the horizon, and the aim 
is to achieve a form of h(a) that is globally valid. 
Differentiating the definition (6.102) we find that 

2M . 
f(a) = a + - - 2M a ' e ' r &t (6.105) 

and hence 

2M 

f _1 (a) = a — — a-e t e r . (6.106) 

r — 2M 

Now, the function (6.103) is independent of t, so h(a;x') = h(a;x). It follows 
that the transformed function h'(a) is given by 

2 M 

h'(a) = hf-^a) = h(o) A-^ 2 a-e t e r . (6.107) 

r 

This new solution is not yet well-defined for all r, but if we now apply the boost 
defined by the rotor 

fl = exp(cr r x/2), (6.108) 

where 

sinhx = ^(A~ 1/2 - A 1 / 2 ), (6.109) 

we obtain the solution 

h"(o) = a+ — a-e-e-, (6.110) 
r 

where 

e_ = e t — e r . (6.111) 

The solution (6.110) is now globally defined. It is the GTG equivalent of the 
Kerr-Schild form of the Schwarzschild solution, and has the property that in- 
falling null geodesies are represented by straight lines on a t—r plot. It is not 
hard to find a transformation between (6.110) and the Newtonian gauge solu- 
tion (6.79). This transformation consists of a displacement and a rotation, both 
of which are globally well-defined. On the other hand, if one starts with the 
solution (6.110) and tries to recover a version of the Schwarzschild solution by 
working in reverse, it is clear that the process fails. The boost needed is infi- 
nite at the horizon and ill-defined for r < 2M, as is the required displacement. 
Such transformations fail to meet our requirement that gauge transformations 
be well-defined over the whole region of physical interest. 

Integral equations and the singularity 

The Riemann tensor 7Z{B) contains derivatives of terms from uj{a) which fall 
off as 1/r 2 . When differentiating such terms, one must take account of the fact 
that 

V-{x/r 3 ) = 4ir6{x), (6.112) 
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where x = xAet- This fact will not affect the fields away from the origin, 
but will show up in the results of integrals enclosing the origin. To see how, 
we again return to the setup before the Newtonian gauge was chosen. From 
equation (6.34) we see that 

e t -g(e t ) = 8np = 2{—L r T — T 2 + FS) - AA (6.113) 

and, using equations (6.46) and (6.93), this gives 

4np = -{gid r gx - g2d r g 2 )/r + M/r 3 
= d r {M/r)/r + M/r 3 
= (o- r /r)-V(M/r) + (M/r)V-(a r /r) 
= V-{Mx/r 3 ). (6.114) 

It follows that p = M6(x), so the singularity generating the radially-symmetric 
fields is a simple <5-function, of precisely the same kind as the source of the 
Coulomb field in electrostatics. 

The presence of the 5-function source at the origin is most easily seen when 
the solution is analysed in the gauge defined by (6.110). Solutions of this type 
are analysed in [55], and we restrict ourselves here to a few basic observations. 
For the solution (6.110), Tl{B) is given by 

K(a + lb) = M(a- V(x/r 3 ) + IV-{bhx/r 3 )) (6.115) 

and it is simple to see that, away from the origin, (6.115) reduces to (6.73). In 
this section we adopt the convention that, when all symbols are in bold face, 
the inner and outer products drop down to their three-dimensional definitions. 
The significance of (6.115) is that it allows us to compute the integral of the 
Riemann tensor over a region enclosing the origin simply by converting the 
volume integral to a surface integral. Taking the region of integration to be a 
sphere of radius ro centred on the origin, we find that 

/ d 3 x1l(a) = M / d<f> dO sin6> a cr r a r = -a, (6.116) 

Jr<r JO JO 3 

and 

/ d 3 xTZ(Ib)= d<p dd sm6I<T r -(bA<T r ) = — ^— lb. (6.117) 

Jr<r a Jo Jo 3 

These results arc independent of the radius of the spherical shell, reflecting the 
spherical symmetry of the solution. The above results combine to give 

/ d 3 xTZ(B) = ^!-B-4:nMB/\e t e t = -^!-(B + 3etBet), (6.118) 

Jr<r " " 

which makes it clear what has happened. The angular integral of the Weyl 
component of 71(B) has vanished, because 

/ d(j> / dO sm6(B + 3crBcr r ) = 0, (6.119) 
Jo Jo 
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and what remains is the contribution from the stress-energy tensor, which is 
entirely concentrated at the origin. On contracting we find that 



and, since TZ(a) = everywhere except for the origin, the integrals in (6.120) 
can be taken over any region of space enclosing the origin. It is now apparent 
that the solution represents a point source of matter, and we can therefore write 



for the matter stress-energy tensor. This is consistent with the definition of M 
as the integral of the density for a static system (6.62). 

Analysing singularities in the gravitational fields by means of integral equa- 
tions turns out to be very powerful in GTG. While the above application does 
not contain any major surprises, we show in [55] that the same techniques ap- 
plied to axisymmetric fields reveal that the Kerr solution describes a ring of 
rotating matter held together by a disk of isotropic tension — a quite different 
picture to that arrived at in general relativity. This clearly has implications for 
the ultimate fate of matter falling onto the singularity, and could yield testable 
differences between GTG and general relativity. 

6.5 Collapsing dust 

The equations in Table 4 can be used to determine the future evolution of a 
system given an equation of state and the initial p and 172 distributions. They 
are therefore well-suited to the description of radial collapse and the forma- 
tion of horizons and singularities. The simplest model, in which the pres- 
sure is set to zero, describes collapsing dust. This situation was first studied 
by Oppenheimer & Snyder [71] and has been considered since by many au- 
thors [63, 72, 73, 74]. A feature of these studies is the appearance of formulae 
which have a suggestively Newtonian form. This is usually dismissed as a 'coin- 
cidence' [73, Section 32.4]. Here we study the collapse process in the Newtonian 
gauge and show that, far from being coincidental, the Newtonian form of the 
results is a natural consequence of the equations. The distinguishing feature of 
the Newtonian gauge approach is that the associated line element is not diag- 
onal. This manifestly breaks time-reversal symmetry, as is appropriate for the 
description of collapsing matter. Working in this gauge enables us to keep all 
fields globally defined, so the horizon is easily dealt with and the matching onto 
an exterior vacuum region is automatically incorporated. This is quite different 
from previous work [71, 72, 73], which usually employs two distinct diagonal 
metrics, one for the matter region and one for the vacuum. Finding the cor- 
rect matching conditions between these metrics is awkward, and difficulties are 
encountered once the horizon has formed. 

If p = it follows immediately that G — and hence /1 = 1 . This ensures 
that the global time coordinate t agrees with the time measured by observers 




(6.120) 



T(a) = M8(x)a-et et 



(6.121) 
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comoving with the fluid. Since v-Vv = in the absence of pressure, such 
observers are also freely falling. The function gi defines a velocity since, for a 
particle comoving with the fluid, 172 is the rate of change of r (which is defined 
by the Weyl tensor) with proper time t. To emphasise its role as a velocity we 
replace 92 with the symbol u for this section. The equations of Table 4 now 
reduce to 

F = d r u (6.122) 
M(r,t)= [ 4irr' 2 p{r',t)dr\ (6.123) 



which define F and M on a time slice, together with the update equations 

d t u + ud r u = -M/r 2 (6.124) 
d t M + ud r M = 0. (6.125) 

Equations (6.124) and (6.125) afford an entirely Newtonian description of the 
fluid. Equation (6.124) is the Euler equation with an inverse-square gravi- 
tational force, and (6.125) is the equation for conservation of mass. The L t 
derivative plays the role of the 'matter' or 'comoving' derivative for the fluid 
since, when acting on a scalar, v-T> — L t . 

The fact that L t M = in the absence of pressure (6.125) is a special case 
of a more general result. Consider the integral 

I(r,t)= [ 4irs 2 p(s,t)f(s,t)ds (6.126) 
Jo 

where f(r, t) is some arbitrary function that is conserved along fluid streamlines, 
that is, it obeys 

L t f = 0. (6.127) 
If we now construct L t I we find that 

L t I = u4irr 2 pf + f Airs 2 (pdj + fd t p) ds. (6.128) 
Jo 

But, from 

L t p = -(2u/r + F)p (6.129) 
and equation (6.122), we have 

d t {r 2 p) = -d r (ur 2 p). (6.130) 

Similarly, from equation (6.127), we see that 

d t f = -ud r f, (6.131) 

so it follows that 

L t I = uAirr 2 pf - [ 4tt(us 2 pdj + fd s (us 2 p)) ds = 0. (6.132) 
Jo 

Any integral of the type defined by I leads to a quantity that is conserved 
along the fluid streamlines. The integral for M(r, t) (6.123) is one such example, 
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with / set to 1. It is clear from its appearance in the Ricmann tensor that M 
represents the 'gravitating energy' of the region enclosed inside r. 
Since G = 0, we have 

Ltgi = (6.133) 
and an alternative conserved quantity is therefore defined by 

r ds 
(j,(r,t)= 4Trs 2 p(s,t) —. (6.134) 
Jo 3i 

This is the covariant integral of the density, so is also a covariant scalar quantity; 
it is simply the total rest-mass energy within r (sec Box 23.1 of 'Gravitation' [73] 
for a discussion of this point in the static case). The relationship between the 
rest-mass energy /i and the gravitating energy M can be seen more clearly by 
recalling that 

9l 2 = l- 2M/r + M 2 . (6.135) 

Since 

r ds 
M(r,t) - ii(r,t) = 4ns 2 p(s,t)( gi -l)—, (6.136) 
Jo 9i 

the difference between the rest energy p and the total energy M is governed by 
— This is then multiplied by the term Airr 2 pdr / g\, which is the rest mass 
of a shell of width dr. For \2M/r — u 2 \ <C 1 wc can approximate (6.135) to give 

3i - 1 ps -M/r + \u 2 (6.137) 

which explicitly shows the decomposition of the energy difference into the sum of 
the Newtonian gravitational potential energy (always negative) and the energy 
due to the bulk kinetic motion (always positive). It is clear that for a shell of 
material to escape it must have 31 — 1 > so, with no approximation necessary, 
we recover the Newtonian escape velocity u 2 = 2M/r. 

As a further example of the insight provided by the Newtonian gauge, con- 
sider the case where the interior of the shell is empty. In this case M = 0, 
so 

31 = (l + u 2 ) 1/2 , (6.138) 

which shows that 31 can be interpreted as a relativistic 7-factor associated u. 
This identification is justified if we put u = sinha, which is reasonable since 
we know that u can be greater than 1. It is the presence of this additional 
boost factor in the formula for M compared to p that, in this case, makes the 
total gravitating energy greater than the rest mass energy. These results should 
demonstrate that the physical picture of gravitational collapse in the absence 
of pressure is really no different from that afforded by Newtonian physics and 
special relativity. It is therefore no surprise that many of the results agree with 
those of Newtonian physics. Furthermore, abandoning a description in terms of 
distorted volume elements and spacetime geometry has allowed us to recover a 
clear physical picture of the processes involved. 

Analytical solutions 

A useful property of the system of equations obtained when p = is that it 
is easy to construct analytical solutions [75, 76]. To see this in the Newtonian 
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gauge we write equation (6.125) in the form 



Since M is a function of r and t only we can employ the reciprocity relation 

-1 (6.140) 



to deduce that 



fdM\ fdt\ f dr\ 
\dt) r \d~r) M \dM) t 

(-) 



(6.141) 



dr). 



■ M 

But we know that u is determined by equation (6.135), and we also know that 
both M and g\ are conserved along the fluid streamlines. We can therefore 
write gi = gi(M), and equation (6.141) can be integrated straightforwardly to 
give t as a function of r and M. 

To perform the integration it is necessary to make a choice for the sign of 
u. For collapsing matter we clearly require u < 0, while for cosmology it turns 
out that u > is the appropriate choice [77] . For this section we can therefore 
write 

= -(gi(M) 2 - 1 + 2M/r)~ 1/2 . (6.142) 

' M 

Finally, we need to choose a form for g\. This amounts to making an initial 
choice of u, since u and g\ are related via equation (6.135). For this section we 
simplify to the case in which the matter is initially at rest. This might provide 
a reasonable model for a star at the onset of a supernova, in which there is a 
catastrophic loss of pressure due to vast amounts of neutrino production, and 
the central core is suddenly left with no supporting pressure. With u(r, 0) = 
we can write 

5l 2 = l-2A/(r )/r (6.143) 

where ro labels the initial r coordinate. We can view the value of ro as carried 
along the streamline defined by it at t = 0, so can write ro = ro(t, r) and treat 
M and g\ as functions of ro only. Equation (6.142) now becomes 



dt\ /2M 2M\ 



-1/2 



dr) V r 



ro J 



(6.144) 



which integrates to give 
3 \ V2 



t = 



= (Sf) ( 7r / 2 - sin ' 1 (^Ao) 1/2 + (r/ro) 1/2 (l-r/r ) 1/2 ) (6.145) 



where we have chosen the initial conditions to correspond to t = 0. 

Equation (6.145) determines a streamline for each initial value ro, and can 
therefore be treated as implicitly determining the function ro(r, t). Since M(tq) 
and .gi (ro) are known, the future evolution of the system is completely deter- 
mined. Furthermore, quantities such as p or d r u can be found directly once 
ro(r, t) is known. The above approach is easily extended to deal with initial 
conditions other than particles starting from rest since, once M(tq) and <?i(ro) 
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are known, all one has to do is integrate equation (6.142). The ability to give a 
global description of the physics in a single gauge allows for simple simulations 
of a wide range of phenomena [77] . 

An important restriction on the solution (6.145) is that the streamlines 
should not cross. Crossed streamlines would imply the formation of shock fronts, 
and in such situations our physical assumption that p = will fail. Streamline 
crossing is avoided if the initial density distribution p(r ) is chosen to be either 
constant or a monotonic-decreasing function of ro- This is physically reasonable 
and leads to sensible simulations for a collapsing star. 



Singularity formation 

An immediate consequence of equation (6.145) is that the time taken before a 
given streamline reaches the origin is given by 

h =7r(r3/8M) 1/2 . (6.146) 

Since the global time t agrees with the proper time for observers comoving with 
the fluid, equation (6.146) is also the lapse of proper time from the onset of 
collapse to termination at the singularity as measured by any particle moving 
with the dust. As pointed out in Section 32.4 of [73], the formula (6.146) agrees 
with the Newtonian result. For the reasons given above, this should no longer 
be a surprise. 

Since <?i and M are conserved along a streamline, equation (6.135) shows 
that u = <72 must become singular as r — * 0. Thus the central singularity forms 
when the first streamline reaches the origin. Near = the initial density 
distribution must behave as 

p(r ) « po - O(r 2 ) (6.147) 
so the mass function M{r$) is 

47T 

M(r ) = yr 3 - O(r 5 ). (6.148) 



It follows that 

rl/8M) 1/2 =^'^ 1/2 
so the central singularity forms after a time 



lim 7r(rg/8Af) ' = (37r/32p ) (6.149) 



"-(^r <6 - i5o) 

where po is the initial density at the origin and the gravitational constant G has 
been included. 

A simulation of this process is shown in Figure 2, which plots the fluid 
streamlines for the initial density function 

9=1-, (6.151) 

(l + (r/a) 2 ) 2 v ' 

The streamlines all arrive at the singularity after a finite time, and the bunching 
at tf can be seen clearly. Solutions can be extended beyond the time at which 
the singularity first forms by including an appropriate 5-function at the origin. 
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Figure 2: Fluid streamlines for dust collapsing from rest. The initial density is 
given by equation (6.151) with po = 0.22 and a = 1. The singularity Grst forms 
at t f = 1.16. 

A further point revealed by such simulations is that it is possible to have 
quite large differences between the total rest mass /i^ = /i(r = oo) and the final 
mass of the singularity, M . For example, for the case plotted in Figure 2 we 
find that fioo = 6.26, whereas M = 2.17. In this case nearly 3 times as many 
baryons have gone into forming the black hole than is apparent from its mass 
M. The possibility of large differences between M and fi is usually ignored in 
discussions of the thermodynamics of black holes (see footnote 14 of [78]). 

Horizon formation 

Any particle on a radial path has a covariant velocity of the form 

v = coshw e t + sinhu e r . (6.152) 

The underlying trajectory has x = h(t>), so the radial motion is determined by 

f = coshu .92 + sinhu gi . (6.153) 

Since gi is negative for collapsing matter, the particle can only achieve an out- 
ward velocity if g\ 2 > 52 2 ■ A horizon therefore forms at the point where 

2M(r,t)/r=l. (6.154) 

To illustrate the formation of a horizon, we again consider the initial density 
profile of equation (6.151). By inverting (6.145) at fixed t, ro is found as a 
function of r. From this, (1— 2A/(r, t) /r) is calculated straightforwardly, and this 
quantity is plotted on Figure 3 at different time slices. The plots show clearly 
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Figure 3: Simulation of collapsing dust in the Newtonian gauge. Successive time 
slices for the horizon function (1 — 2M(r, t)/r) versus x are shown, with the top 
curve corresponding to t = and lower curves to successively later times. The 
initial velocity is zero, and the initial density prohlc is given by equation (6.151) 
with po = 0.22 and a = 1. There is no horizon present initially, but a trapped 
region quickly forms, since in regions where 1 — 2M/r < photons can only 
move inwards. 

that the horizon forms at a finite distance from the origin. It is conventional to 
extend the horizon back in time along the past light-cone to the origin (r = 0), 
since any particle inside this surface could not have reached the point at which 
1 — 2M/r first drops to zero, and hence is also trapped [79]. The ease with 
which horizon formation is treated again illustrates the advantages of working 
in a non-diagonal gauge. Such considerations will clearly be important when 
performing numerical studies of more realistic collapse scenarios. 

A final point is that, since it is negative, it follows that fig% — f%g\ = u must 
also be negative. This tells us that the fields that remain after the collapse 
process has finished are in the class defined by /152 ~ ./2.91 = — 1 at the horizon. 
This time direction is then frozen into the fields, as discussed in Section 6.4. 

6.6 Cosmology 

The equations of Table 4 are sufficiently general to deal with cosmology as well 
as astrophysics. In recent years, however, it has once more become fashionable 
to include a cosmological constant in the field equations. The derivation of 
Section 6.2 is largely unaffected by the inclusion of the cosmological term, and 
only a few modifications to Table 4 arc required. The full set of equations with 
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The h-function 



h(e*) 


= /ie* 


h(e r ) 




He 6 ) 


= e e 


h(e*) 





The w function 



w(et) 


= Ge r et 




w(e r ) 


= -Fe r e t 






= 92/rOet -\ 


- (51 - l)/re r 6 




= 92/r^et - 


\- (91 - l)Ae r 



Directional derivatives 



it = ji<9 t + ,g 2 9 r 
i r = g\d r 



G, F and fx 



Definition of M 



L t gi = Gg 2 , 
fi = exp(/ r 


L r g2 = 
-G/gi dr) 


= F 9l 




M = \r{g 2 2 - 


9i 2 + l- 


Ar 2 /3) 



Remaining derivatives 


Lt 32 = 


G91- 


M/r 2 + 


rA/3- 


- Anrp 






L r gi = 


Fg 2 + 


M/r 2 - 


rA/3- 


- Aitrp 


















L t M = 


-47T.9; 


r 2 p, 


L r M 


= Airgir 2 


P 


Matter derivatives 


L t p = 


"(2.92/ 


r + F){p 


+ P) 






LrP = 


-G(p- 


hp) 







Ricmann tensor 



n{B)An{p+p)B-e t e t - ±(87175 + A)5 
-±(M/r 3 - 47rp/3)(B + 3o- r Bcr r 



Stress-energy tensor T{a) = (p + p)a-e t e-t — pa 



Table 5: Equations governing a radially-symmetric perfect fluid with a non-zero 
cosmological constant A. The shaded equations differ from those of Table 4- 



a cosmological constant incorporated are summarised in Table 5. 

In cosmology we are interested in homogeneous solutions to the equations of 
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Table 5. Such solutions are found by setting p and p to functions of t only, and 
it follows immediately from the L r p equation that 

G = 0, and /i = 1. (6.155) 

For homogeneous fields the Weyl component of the Riemann tensor must vanish, 
since this contains directional information through the e r vector. The vanishing 
of this term requires that 

4 , 

M(r,t) = -Trr 3 p, (6.156) 
o 

which is consistent with the L r M equation. The L t M and L t p equations now 
reduce to 

F = g 2 /r (6.157) 

and 

p = -3g 2 (p + p)/r. (6.158) 
But we know that L r g 2 = Fg±, which can only be consistent with (6.157) if 

F = H(t), g 2 (r,t)=rH(t). (6.159) 

The L t g2 equation now reduces to a simple equation for if, 

H + H — A/3 = — -(p + 3p). (6.160) 
Finally, we are left with the following pair of equations for g\\ 

L t gi = (6.161) 
Lrfh = (gi 2 - l)/r. (6.162) 

The latter equation yields gi 2 = 1 + r 2 cf>(t) and the former reduces to 

4> = -2H{t)<p. (6.163) 

Hence g\ is given by 

gi 2 = l-kr 2 exp (-2 J H(t') dt'^j , (6.164) 

where k is an arbitrary constant of integration. It is straightforward to check 
that (6.164) is consistent with the equations for H and p. The full set of equa- 
tions describing a homogeneous perfect fluid arc summarised in Table 6. 

At first sight, the equations of Table 6 do not resemble the usual Friedmann 
equations. The Friedmann equations are recovered straightforwardly, however, 
by setting 

H(t) = §| (6.165) 
With this substitution we find that 

9l 2 = 1 - kr 2 /S 2 (6.166) 
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The h-function 



h(a) = a + a-e r (( gi - l)e r + H^re 1 
<?i 2 = 1 - kr 2 exp(-2 /* H(t') dt') 



The u) function u(a) = H(t)aAe t — (gi — \)/ra/\(e r et)et 



The density 

Dynamical 
equations 



87rp/3 = H{t) 2 - A/3 + fccxp(-2 /' H(t') dt') 



H + H 2 - A/3 = -4tt/3 (p + 3p) 
p = -3H(t)(p + p) 



Table 6: Equations governing a homogeneous perfect fluid. 



and that the H and density equations become 

S A 47T, 

s" 3 =-Y {p + 3p) 



and 



g 2 + fc 

^ 2 



A 

3 



Stt 



-P- 



(6.167) 



(6.168) 



recovering the Fricdmann equations in their standard form [80]. The intrinsic 
treatment has therefore led us to work directly with the 'Hubble velocity' H(t), 
rather than the 'distance' scale S(t). There is a good reason for this. Once the 
Wcyl tensor is set to zero, the Riemann tensor reduces to 



11(B) = 4ir(p + p)B-e t e t — |(87rp + A)B, 



(6.169) 



and we have now lost contact with an intrinsically-defined distance scale. We 
can therefore rescale the radius variable r with an arbitrary function of t (or r) 
without altering the Riemann tensor. The Hubble velocity, on the other hand, 
is intrinsic and it is therefore not surprising that our treatment has led directly 
to equations for this. 

Among the class of radial rescalings a particularly useful one is to rescale r 
to r' = S(t)r. This is achieved with the transformation 

f(x) = x-etet + SxAetet, (6.170) 

so that, on applying equation (3.9), the transformed h-function is 

h'(o) = a-e t e t + -^((1- kr 2 ) 1/2 a-e r e r + aAa r oy) . (6.171) 



The function (6.171) reproduces the standard line element used in cosmology. 
We can therefore use the transformation (6.170) to move between the 'Newto- 
nian' gauge developed here and the gauge of (6.171). This is useful for later 
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sections, where the Maxwell and Dirac equations are solved in a cosmological 
background described by (6.171). The differences between these gauges can 
be understood by considering geodesic motion. A particle at rest with respect 
to the cosmological frame (defined by the cosmic microwave background) has 
v = et- In the gauge of (6.171) such a particle is not moving in the flatspace 
background (the distance variable r is equated with the comoving coordinate). 
In the Newtonian gauge, on the other hand, comoving particles are moving out- 
wards radially at a velocity r = H(t)r, though this expansion centre is not an 
intrinsic feature. Of course, attempting to distinguish these pictures is a point- 
less exercise, since all observables must be gauge-invariant. All that is of physical 
relevance is that, if two particles are at rest with respect to the cosmological 
frame (defined by the cosmic microwave background) , then the light-travel time 
between these particles is an increasing function of time and light is redshiftcd 
as it travels between them. 

Dust models 

The utility of the Newtonian gauge in cosmology has been independently dis- 
covered by other authors [62, 81]. Here we illustrate its advantages for dust 
models (p = 0). Setting p to zero implies that 

H(t) = -p/3p, (6.172) 

so 

gi 2 = 1 - kr 2 p 2/3 (6.173) 

and 

1/2 

H(t)=(^p-k^ + ^j ■ (6.174) 

We are therefore left with a single first-order differential equation for p. Explicit 
solutions of this equations are often not needed, as we can usually parameterise 
time by the density p(t). 

If we now look at the trajectories defined by the fluid, these have v = et 
which implies that 

x = e t +rH(t)e r . (6.175) 

It follows that 

r/r = H{t) = -p/'Sp (6.176) 

and hence that 

r/r = (p/p r 1/3 - (6.177) 

The fluid streamlines form a family of spacetime curves spreading out from the 
origin at the initial singularity (when p = oo). The Newtonian gauge therefore 
describes an expanding universe in a very simple, almost naive manner. Since all 
points in a homogeneous cosmology are equivalent, we can consider ourselves to 
be located at r = 0. The Newtonian gauge then pictures our observable universe 
as a ball of dust expanding outwards radially from us. 

While the picture provided by the Newtonian gauge has no physical reality 
of its own, it does have some heuristic merit and can provide a useful aid to 
one's physical intuition. For example, consider the familiar relationship between 
angular size and redshift. An initially surprising feature of this relationship is 
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Figure 4: Null geodesies for a dust-filled k = A = universe in the Newtonian 
gauge. All geodesies start from the origin at the initial singularity (here set at 
t = 0). The dashed line gives the critical distance where photons 'turn round' 
in this gauge. 

that, beyond a certain redshift, angular sizes stop decreasing and start increas- 
ing. This result is easily understood in the Newtonian gauge. Consider the 
photon paths shown in Figure 4. (These paths are for a k = and A = 
universe, though the comments are applicable more generally.) Suppose that 
at some finite time to we receive a photon from the distant past. This photon 
must have followed part of a path which begins on the origin at t = 0. Before 
a certain time in the past, therefore, photons received by us must have initially 
travelled outwards before turning round and reaching us. The angular size of an 
observed object is then that appropriate to the actual r coordinate of the object 
when it emitted the photons. Since the value of r decreases before a certain 
time, angular sizes appear larger for objects that emitted photons before this 
time. 

A radial null geodesic has a trajectory 



t{r)e t + r(r)e r . 



(6.178) 



For these trajectories 




(6.179) 



and the condition that v 2 = reduces to 



f/ 5 i = i(rH(t)/ 9l ± 1) 



(6.180) 




dr 



(6.181) 
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Since H (t) is positive in a cooling universe, the distance at which photons 'turn 
round' is defined by 



9i _ (Sirp AN 1/2 
H(t) V 3 3 



(6.182) 



(In [59] we mistakenly referred to this distance as a particle horizon.) For 
k = A = it is simple to show that this distance corresponds to a redshift of 
1.25. This result relates physically measurable quantities, so is gauge invariant. 



Closed universe models 



The Newtonian gauge presents a particularly simple picture for a closed universe 
(k > 0). For k > the requirement that 1 — fcr 2 p 2 / 3 is positive means that 

r 3 p > fc~ 3/2 . (6.183) 

This places a limit on the speed with which the dust can expand, so the solution 
describes a finite ball of dust expanding into a vacuum. This ball expands out 
to some fixed radius before turning round and contracting back to the origin. 
The turning point is achieved where H(t) = so, for A = 0, this occurs when 

p=(3fc/87r) 3 . (6.184) 

The maximum radius is therefore 

87T 

rmax = ^572 . (6.185) 

The time taken to reach the future singularity is given by (6.150), since this 
cosmological model is a special case of spherical collapse in which the density 
is uniform. This picture of a model for a closed universe is both simple and 
appealing. It allows us to apply Newtonian reasoning while ensuring consistency 
with the full relativistic theory. 

The finite ball model for a k > cosmology is clearly useful when considering 
experiments with particles carried out near the origin, but globally one must 
consider the boundary properties of the ball. A crucial question is whether the 
particle horizon (the largest region of the universe with which an observer at 
the origin is in causal contact) extends past the edge of the ball or not. It can 
be shown that this horizon always lies inside the radius at which g\ becomes 
imaginary, except at the turnaround point (the point at which the ball reaches 
its maximum radius), where the two radii coincide. A suitable choice of cutoff 
radius is therefore available in either the expanding or contracting phase sep- 
arately, but what happens at the turnaround point is potentially ambiguous. 
When discussing field theory, however, the finite ball model is inadequate. One 
must instead use a global gauge so, in Section 7.4, we introduce the 'stereo- 
graphic projection gauge'. This provides a global solution which can be shown 
to be spatially closed. This solution is used in the study of electromagnctism 
(Section 7.4) and the Dirac equation (Section 8.3) in a cosmological background. 
It is possible to treat the stereographic projection gauge solution in a form of 
the Newtonian gauge, though this possibility is not explored here. 
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7 Electromagnetism in a gravitational 
background 

In Section 4 we derived field equations for the gravitational and Dirac fields. 
We now turn to the derivation of the Maxwell equations in a gravitational 
background. A number of applications of these equations are then discussed, 
including a simple derivation of the characteristic surfaces for both the Maxwell 
and Dirac equations. 

The basic dynamical variable is the electromagnetic vector potential A, for 
which the coupling to spinor fields was derived in Section 3.3. Under phase 
rotations of the spinor field, A transforms as 

eA eA - V<£. (7.1) 

It follows that, under a displacement, A must transform in the same way as V0, 
that is, 

A[x) h- > l(A(x')). (7.2) 
The covariant form of the vector potential is therefore 

A=h(A), (7.3) 

which is the term that appeared in the Dirac equation (3.63). 
From A, the Faraday bivector F is defined by 

F = VAA (7.4) 

This definition implies that, under displacements, F (x) is transformed to F'(x), 
where 

F'(x) = WMA(x') 

= ~f(V x ,AA(x')) 

= fF(x'). (7.5) 

It follows that the covariant analogue of F is defined by 

T=h(F), (7.6) 

which is covariant under position and rotation-gauge transformations, and is 
invariant under phase changes. 

The same covariant quantity T is obtained if one follows the route used for 
the construction of R(a A b) at (4.2). In particular, the contracted commutator 
of two covariant derivatives gives (in the absence of torsion) 

h^Ah^)^, A,]V> - ih(d b )h(d a )R(aA6)V- 

= \RaP. (7.7) 

The analogous construction for the 'internal' covariant derivative 

D J a ip = a-Vip - ea-AtpIa 3 (7.8) 
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gives 

h (e" ) A h {e» ) [Dl , D* v ] i> = eh (d b ) A h (0„ ) (a A 6) • ^W<r 3 

= 2eT^I(T 3 , (7.9) 

which clearly identifies J as a covariant quantity. Unlike 1Z, however, T is a 
bivector and equation (7.9) exhibits a curious interaction between this bivector 
on the left of ip, an d the fixed bivector icr 3 on the right. 

Having defined the covariant bivector J 7 , it is clear that the appropriate 
generalisation of the electromagnetic action to include gravitational interactions 
is 

S= J |d 4 x|det(h)- 1 {\T-T-A-J), (7.10) 

where here J is the covariant charge current. Unlike the Dirac action, the 
rotation-gauge field f2(a) does not appear in this action. It follows that the 
electromagnetic field does not act as a source of spin. The action (7.10) is 
varied with respect to A, with h(a) and J treated as external fields. The result 
of this is the equation 

VG = J, (7.11) 

where 

G= hh(F)det(h)- 1 (7.12) 

and 

J=h(J)det(b)- 1 . (7.13) 

Equation (7.11) combines with the identity 

VAF = (7.14) 

to form the full set of Maxwell equations in a gravitational background. We 
again see that the classical field equations can be expressed in a form from 
which all reference to the rotation gauge has been removed. 

Some insight into the equations (7.11) and (7.14) is obtained by performing 
a space-time split (see Section 2.1) and writing 

F = E + IB (7.15) 
G = D + IH (7.16) 
J l0 = P + J. (7.17) 

In terms of these variables Maxwell's equations in a gravitational background 
take the familiar form 

V-B = V-D = p 

^ dB „ TT T dD (7-18) 

VxE = — — VxH = J+ — 

dt dt 

where V = 70 A V = <Jid x i is the 3D vector derivative, and the bold cross X is 
the traditional vector cross product: 

axb = -Iaxb. (7.19) 
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Equation (7.18) shows that the h-field defines the dielectric properties of the 
space through which the electromagnetic fields propagate, with det(h) _1 h h de- 
termining the generalized permittivity/permeability tensor. Many phenomena, 
including the bending of light, can be understood easily in terms of the proper- 
ties of the dielectric defined by the h-field. 



Covariant form of the Maxwell equations 

So far we have failed to achieve a manifestly covariant form of the Maxwell 
equations. We have, furthermore, failed to unite the separate equations into a 
single equation. In the absence of gravitational effects the equations V ■ F = J 
and VA-F = combine into the single equation 

VF = J. (7.20) 

The significance of this equation is that the V operator is invertible, whereas 
the separate V- and VA operators are not [11, 26]. Clearly, we expect that 
such a unification should remain possible after the gravitational gauge fields 
have been introduced. To find a covariant equation, we first extend the 'wedge' 
equation (4.37) to include higher-grade terms. To make the derivation general 
we include the spin term, in which case we find that 

X>Ah~(aA&) = (X>Ah(a))Ah(6) - h(a) A(X>Ah(6)) 
= K{b(a)-S)Ab{b) - Kb{a)A(b{b)-S) 
= KSxh(aAb) (7.21) 

and, more generally, we can write 

VAb(A r ) = h(VAA r ) + K(Sb(A r )) r+1 . (7.22) 

We can therefore replace equation (7.14) by 

VAT - kSxT = 0. (7.23) 

We next use the rearrangement 

V-(h(J r )dct(h)- 1 ) = 7VA(/h(J : -)dct(h)- 1 ) 
= 7VA(h- 1 (/JC)) 

= Ib~- 1 (VA{IT) + k{IT)xS) (7.24) 

to write equation (7.11) as 

V-T - nS-T = Ib{JI) = J. (7.25) 

Equations (7.23) and (7.25) now combine into the single equation 

VT - kST = J, (7.26) 

which achieves our objective. Equation (7.26) is manifestly covariant and the 
appearance of the T>T term is precisely what one might expect on 'minimal- 
coupling' grounds. The appearance of the spin term is a surprise, however. 
Gauge arguments alone would not have discovered this term and it is only 



81 



through the construction of a gauge-invariant action integral that the term is 
found. Equation (7.26) should be particularly useful when considering electro- 
magnetic effects in regions of high spin density, such as neutron stars. 

To complete the description of electromagnetism in a gravitational back- 
ground we need a formula for the free-field stress-energy tensor. Applying the 
definition (4.20) we construct 

TUb-'ia)) = idet(h)9 R(a) (h(F)h(F)det(h)- 1 ) 

= h(a-F)-F-±h- 1 {a)J 7 -J r . (7.27) 

Hence, 

T cm (a) = h{U(a)-F)-T - \aT-T 
= (a-F)-F — \aT-T 

= -\TaT, (7.28) 

which is the natural covariant extension of the gravitation-free form —FaF/2. 
The tensor (7.28) is symmetric, as one expects for fields with vanishing spin 
density. 

7.1 Characteristic surfaces 

In their spacetime algebra form the Maxwell and Dirac equations are both first- 
order differential equations involving the vector derivative V. For electromag- 
netism, this first-order form of the equations offers many advantages over the 
equivalent second-order theory [14, 26] . We have now seen that gravitational in- 
teractions modify both these equations in such a way that the vector derivative 
V is replaced by the position-gauge covariant derivative h(V). As an illustration 
of the utility of first-order equations, both without and with gravitational effects, 
we now give a simple derivation of the properties of characteristic surfaces. 
Consider, initially, a generic equation of the type 

VV = /(M (7-29) 

where ip is some arbitrary field, and / is some known function. Suppose that an 
initial set of data is given over some three-dimensional surface in spacetime, and 
we wish to propagate this information off the surface into the adjoining region. 
We pick three vectors, a, b and c, which are tangent to the surface. From our 
initial data we can construct a-V?/>, b-Vip and c- V?/>. Now define 

n = IaAbAc (7.30) 

and use nV = n- V + nAV to decompose (7.29) into 

n-Wtp = -nAWtp + nf(ip). (7.31) 

The right-hand side of (7.31) contains the term 

nAVijj = I(aAbAc)-Vi' 

= I(aAbc-V4> - aAcb-Vip + bAca-Vip), (7.32) 
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which is therefore known. It follows that we know all the terms on the right- 
hand side of equation (7.31) and can therefore construct n-V ip. This gives us the 
derivative required to propagate off the surface. The only situation for which 
propagation is impossible is when n remains in the surface. This occurs when 

nA(aAfeAc) = 
=*> nA(nI) = 

=> n-n = 0. (7.33) 

It follows that the characteristic surfaces for any first-order equation of the type 
defined by (7.29) are null surfaces. These considerations automatically include 
the Maxwell and Dirac equations. It is notable how this result follows from 
purely algebraic considerations. 

The generalisation to a gravitational background is straightforward. Equa- 
tion (7.26) is generalised to 

h(V)V = /WO, (7.34) 

and we assume that a gauge choice has been made so that all the fields (apart 
from ?p) are known functions of x. Again, we assume that the initial data consist 
of values for ip over some three-dimensional surface, so we can still determine 
a-~Vip etc. Since 

a-V = h _1 (a)-h(V) (7.35) 
it follows that the vector of interest is now 

Jh- 1 (o)Ah -1 (6)Ah- 1 (c) = Ih-^nl) = h(n) det(h)- 1 . (7.36) 

This time we multiply equation (7.34) by h(n) and find that h(n)-h(V)'0 can 
be constructed entirely from known quantities. We can therefore propagate in 
the h h(n) direction, so we now require that this vector does not lie in the initial 
surface. The analogue of (7.33) is therefore 

hh(re)-7i = 0, or h(n) 2 = 0, (7.37) 

and the characteristic surfaces are now those for which h(n) is null. This is the 
obvious covariant extension of n being a null vector. 



7.2 Point charge in a black-hole background 

The problem of interest here is to find the fields generated by a point source held 
at rest outside the horizon of a radially-symmetric black hole. The h-function in 
this case can be taken as that of equation (6.79). The solution to this problem 
can be found by adapting the work of Copson [82] and Linct [83] to the present 
gauge choices. Assuming, for simplicity, that the charge is placed on the z-axis 
a distance a from the origin (a > 2M), the vector potential can be written in 
terms of a single scalar potential V(r, 8) as 



A = V(r,6) U + ^-^er . (7.38) 
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It follows that 



<9y 1 <9v 

D = (7.40) 

and 

T = -f * r - ^% (tr g + V2MA^). (7.41) 
The Maxwell equations now reduce to the single partial differential equation 

^ cV V ~dr~) + r(r - 2M) sm9 86 \ ~d6 ) ~ ~ P ' ( ' 

where p = q5(x — a) is a 5-function at z = a. This was the problem originally 
tackled by Copson [82] who obtained a solution that was valid locally in the 
vicinity of the charge, but contained an additional pole at the origin. Linet [83] 
modified Copson's solution by removing the singularity at the origin to produce 
a potential V(r, 6) whose only pole is on the z-axis at z = a. Linet's solution is 

= ± (r-Af)( -M)-MWfl + *M 
ar D ar 

where 

D = (r(r - 2M) + (a - Af) 2 - 2(r - M)(a - M) cos6> + A/ 2 cos 2 6») 1/2 . (7.44) 

An important feature of this solution is that once (7.43) is inserted back 
into (7.41) the resultant T is both finite and continuous at the horizon. Fur- 
thermore, since h is well-defined at the horizon, both F and G must also be 
finite and continuous there. Working in the Newtonian gauge has enabled us 
to construct a global solution, and we can therefore study its global properties. 
One simple way to illustrate the global properties of the solution is to plot the 
streamlines of D which, from equation (7.18), is divergenceless away from the 
source. The streamlines should therefore spread out from the charge and cover 
all space. Since the distance scale r is fixed to the gravitationally-defined dis- 
tance, the streamlines of D convey genuine intrinsic information. Hence the 
plots arc completely unaffected by our choice for the g\ or g 2 functions, or in- 
deed our choice of i-coordinate. Figure 5 shows streamline plots for charges 
held at different distances above the horizon. Similar plots were first obtained 
by Hanni & Ruffini [84] although, as they worked with the Schwarzschild metric, 
they were unable to extend their plots through the horizon. The plots reveal 
an effective contraction in the radial direction. It is not hard to show that 
the contraction is precisely that of a particle moving with the free-fall velocity 
(2M/0) 1 / 2 relative to a fixed observer. 

The description presented here of the fields due to a point charge is very 
different from that advocated by the 'Membrane Paradigm' [85]. The mem- 
brane paradigm was an attempt to develop the theory of black holes in a way 
that, as far as possible, employed only familiar physical concepts. In this way, 
gravitational effects could be incorporated correctly without requiring an un- 
derstanding of the full theory The hope was that astrophysicists would adopt 
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Figure 5: Streamlines of the D field. The horizon is at r = 2 and the charge 
is placed on the z-axis. The charge is at z = 3 and z = 2.01 for the top and 
bottom diagrams respectively. The streamlines are seeded so as to reflect the 
magnitude of D. The streamlines are attracted towards the origin but never 
actually meet it. Note the appearance of a l cardiod of avoidance' as the charge 
gets very close to the horizon. The equation for this cardiod is r = M(l + cosO), 
which is found by setting D = when a — 2M. 

this paradigm when modelling regions where black-hole physics could be signif- 
icant, such as at the heart of a quasar. The paradigm works by drawing a veil 
over the horizon (the membrane) and concentrating on the physics outside the 
horizon as seen by observers remaining at a fixed distance (fiducial observers). 
Our view, however, is that it is the Newtonian gauge which provides the clearest 
understanding of the physics of black holes while requiring minimal modification 
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to Newtonian and special-relativistic ideas. Furthermore, writing the Maxwell 
equations in the form (7.18) removes any difficulties in applying conventional 
reasoning to the study electromagnetism in a gravitational background. 

There are other ways that our approach offers advantages over the mem- 
brane paradigm. When applying the membrane paradigm one has to work with 
quantities which are singular at the horizon, and this is hardly a recipe for ap- 
plying traditional intuition! As we have seen, once formulated in the Newtonian 
gauge (or any other gauge admissible in GTG) all physical quantities arc finite. 
The membrane paradigm also warns physicists against producing plots such as 
Figure 5, because such plots depend on the choice of radial coordinate. But 
our intrinsic approach makes it clear that such plots are meaningful, because 
r is determined uniquely by the Ricmann tensor. Presenting the plots in the 
form of Figure 5 enables direct physical information to be read off. In short, 
the simple physical picture provided by our intrinsic method and Newtonian- 
gauge solution disposes of any need to adopt the artificial ideas advocated by 
the membrane paradigm. 

7.3 Polarisation repulsion 

An interesting feature of the above solution (7.43) is the existence of a repulsive 
'polarisation' force [86], one effect of which is that a smaller force is needed 
to keep a charged particle at rest outside a black hole than an uncharged one. 
In their derivation of this force, Smith & Will [86] employed a complicated 
energy argument that involved renormalising various divergent integrals. Here 
we show that the same force can be derived from a simple argument based on the 
polarisation effects of the dielectric described by a black hole. First, however, 
we must be clear how force is defined. In the presence of an electromagnetic 
field the equation of motion for a point particle (4.63) is modified to 

mv = (qj- — mu>(v))-v. (7-45) 

We therefore expect that any additional force should also be described by a 
covariant bivector which couples to the velocity the same way that T does. So, 
if we denote the externally applied force as W, the equation of motion for a 
neutral test particle becomes 

mv = (W — mw(v))-v. (7.46) 

Now, suppose that W is chosen so that the particle remains at a fixed distance 
a outside the horizon of a black hole. The equation for the trajectory is 

x(t) = t(r)et + ae r (9 , fa), (7.47) 

where the constants 6q and 4>q specify the angular position of the particle. The 
covariant velocity is 

v = ih- 1 (e t ), (7.48) 
and the condition that v 2 = 1 forces 

i 2 (l -2M/a) = 1, (7.49) 
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so that t, and hence v, arc constant. Equation (7.46) now reduces to 



W - mu)(v) = (7.50) 

-1/2 

<T a , (7.51) 



v./ x Mm / 2M\ 

=S> W = mtfi(et) = — =- 1 

or \ a / 

where <x a is the unit outward spatial vector from the source to the particle. The 
magnitude of the force is therefore 

Mm ( 2M X ~ 1/2 



\W\ = — 1 (7.52) 

a z \ a J 

and it is not hard to check that this result is gauge-invariant. Equation (7.52) 
reduces to the Newtonian formula at large distances and becomes singular as 
the horizon is approached, where an infinite force is required to remain at rest. 

We now want to see how this expression for the force is modified if the particle 
is charged and feels a force due to its own polarisation field. From (7.45) the 
extra term in the force is simply T , and only the radial term contributes. From 
equation (7.41) this is just — d r V <r r . Since the charge lies on the z axis, we 
need only look at V along this axis, for which 

V W -« J" + (7.53) 
\z — a\ \z — a\ \a zj az 

The singular terms must be due to the particle's own Coulomb field, and so 
do not generate a polarisation force. The only term that generates a force is 
therefore the final one, which is precisely the term that Linet added to Copson's 
original formula! This term produces an outward-directed force on the charge, 
of magnitude q 2 M/a 3 . The applied force is therefore now 

„. /ifflA 2M\~ 1/2 Mq 2 \ 

w -{-^{ 1 -— ) -^r)*- (7 ' 54) 

so a smaller force is needed to keep the particle at rest outside the horizon. This 
result agrees with that found in [86] , though our derivation avoids the need for 
infinite mass renormalisation and is considerably simpler. This result is a good 
example of the importance of finding global solutions. The polarisation force is 
felt outside the horizon, yet the correction term that led to it was motivated by 
the properties of the field at the origin. 

7.4 Point charge in a k > cosmology 

We saw in Section 6.6 that one form of the h-function for a homogeneous cos- 
mology is defined by (6.171) 

h(a) = a-e t e t + a((l - fcr 2 ) 1 / 2 a-e r e r + af\cr r ay), (7.55) 

where a = l/S. When k is positive, however, the function (7.55) is undefined 
for r > k~ x / 2 and so fails to define a global solution. A globally- valid solution 
is obtained with the displacement 

f{x) =x-e t e t + 1 + l r2/4 er, (7.56) 
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which results in the particularly simple solution 

h'(a) = a-e t e t + a(l + kr 2 /A)a/\e t e t . (7.57) 

This solution is well-defined for all x and generates an 'isotropic' line element. 
(The solution can also be viewed as resulting from a stereographic projection of 
a 3-sphere.) 

We want to find the fields due to a point charge in the background defined 
by (7.57). Since the h-function is diagonal we start with the obvious ansatz 

A = aV(x)e u (7.58) 

so that 

E = -aVV (7.59) 

and 

D = -(l + fer 2 /4)- 1 W. (7.60) 
It follows that the equation we need to solve is simply 

-V- ((1 + kr 2 /A)- 1 W) = qd(x - a), (7.61) 

where the charge q is located at x = a. This equation can be solved using the 
general technique described by Hadamard (see also Copson [82] for a discussion 
of a similar problem.) The solution to (7.61) turns out to be 

y l + fc«74-fc S /2 

(s(l + fca 2 /4-/cs/4)) 1 /2' v ' > 

where 

(x — a) 2 r^r 

a = Va 2 . (7.63) 



1 + kr 2 /A' 

As a simple check, V reduces to the usual Coulomb potential when k = 0. 

Fieldline plots of the D field defined from (7.62) are shown in Figure 6. 
The fieldlines follow null geodesies and clearly reveal the existence of an image 
charge. The reason for this is can be seen in the denominator of V. Not only is 
V singular at s = 0, it is also singular at 

1 + ka 2 /4-jfes/4 = 0. (7.64) 

So, if the charge lies on the z axis, the position of its image is found by solving 

= 1 + "4 (7-65) 
4 + fcz 2 4 v ' 

from which we find that the image charge is located at 

z = -4/(fca). (7.66) 

(This is the stereographic projection of the opposite point to the charge on a 
3-sphere.) However, if we try to remove this image charge by adding a point 
source at its position, we find that the fields vanish everywhere, since the new 
source has its own image which cancels the original charge. The image charge 
is therefore an unavoidable feature of k > 1 cosmologies. We comment on this 
further in Section 9. 
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Figure 6: Ficldlincs of the D-Rcld for a point charge in a k > universe. The 
hcldlincs follow null geodesies, which are circles in this gauge. The existence of 
an image charge is clear. 

8 The Dirac equation in a gravitational back- 
ground 

In this paper we began by introducing the gravitational gauge fields and min- 
imally coupling these to the Dirac action. For our final major application we 
return to the Dirac equation in a gravitational background. In the absence of 
an A field, the minimally-coupled equation (3.63) is 



Here we consider two applications: the case of a black-hole background; and 
cosmological models in which the universe is not at critical density. 

8.1 Black hole background 

We have already demonstrated that the Newtonian gauge solution dramatically 
simplifies the study of black- hole physics, so this is the natural gauge in which 
to study the Dirac equation. We therefore start with an analysis in this gauge, 
and then consider the gauge-invariance of the predictions made. With the grav- 
itational fields as described in Section 6.4, the Dirac equation (8.1) becomes 



Driers = rmpjo. 



(8.1) 



- (2M/r) 1 / 2 7o (a i .V + 3/(4r)^)Jo- 3 - rrnpjo- 



(8.2) 
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If we pre-multiply by 70 and introduce the symbol i to represent right-sided 
multiplication by ier 3 , so that itp = iplcr^, then equation (8.2) becomes 

idtj) = -iVtJj + i{2M/r) 1/2 r- 3/4 d r (r 3/4 0) + mip, (8.3) 

where ip = 70 "070 • One feature that emerges immediately is that the Newtonian 
gauge has recovered a Hamiltonian form of the Dirac equation (see [27] for a 
discussion of operators, Hamiltonians and Hermiticity in the spacetime algebra 
approach to Dirac theory). Since the Newtonian gauge involves the notion of 
a global time, it might have been expected that this gauge would lend itself 
naturally to a Hamiltonian formulation. An open question of some interest is 
whether it is always possible to make such a gauge choice (see [87, 34] for an 
extension to the Kerr solution). 

The Hamiltonian (8.3) contains a subtlety: it is Hcrmitian only away from 
the origin. To see why, consider the interaction term 

Hity) = i(2M/r) 1/2 r- 3/4 d r (r 3/4 ^). (8.4) 

For this we find that 

d 3 x(^H I (^))s = V2M I dQ I r 2 drr- 5/i {^d r {^ /A tj})I(T S ) t 

dr{r 3 ' 4 ^d r {r 3 l 4 ^)IcT 3 )s 





(8.5) 




d 3 x{(H I (<p)^) s +V2M / dn r 3 ' 2 (tfipltr 3 )s 



where ( ) g denotes the projection onto the 'complex' 1 and ier 3 terms and $ = 
7o07o (see Appendix A and [27]). For all normalised states the final term in (8.5) 
tends to zero as r — > 00. But it can be shown from (8.3) that wavefunctions 
tend to the origin as r~ 3 ' 4 , so the lower limit is finite and Hj is therefore not 
a Hermitian operator. This immediately rules out the existence of normalisable 
stationary states with constant real energy [88]. 

Equation (8.3) can be used to propagate a spinor defined over some initial 
spatial surface, and numerical simulations based on this equation give a good 
picture of the scattering induced by a black hole. Here, however, we wish to 
focus on an analytical approach. A result that follows immediately from the 
Hamiltonian form of the Dirac equation is that (8.3) is manifestly separable, so 
that we can write 

ip(x) = ip(x)a(t). (8.6) 
As usual, the solution of the t-equation is 

a(t) = exp(-7<r 3 Et), (8.7) 

where E is the separation constant. The non-Hermiticity of Hi means that E 
cannot be purely real if ip is normalisable. The imaginary part of E is determined 
by equation (8.5) and, for suitably normalised states, we find that 

Jm(E) = - lim 2 7 rv / 2M(V> t ?/')r 3/2 . (8.8) 
1 — >o 
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This equation shows that the imaginary part of E is necessarily negative, so 
the wavefunction decays with time. This is consistent with the fact that the 
streamlines generated by the conserved current ipJo^P are timelike curves and, 
once inside the horizon, must ultimately terminate on the origin. 
With the i-dependence separated out, equation (8.3) reduces to 

Vtp - (2M/r) 1/2 r- 3/4 d r (r 3/ V) = iEip - imtp. (8.9) 

To solve this equation we next separate out the angular dependence. This 
is achieved using the spherical monogenics, which are the spacetime algebra 
equivalent of the Xjim Pauli spinors. These are described in detail in [27] and 
here we quote the necessary formulae. The unnormalised monogenic ipV 1 is 
defined by 

= ((I + ro + l)P, m (cos(9) - Pl n+1 (cos8)I<T ri ,)e m 'f ,I(T 3, (8.10) 

where I > 0, — (Z + l) < m < I, and P™ are the associated Legendre polynomials. 
The two properties of the ip™ relevant here are 

vvr = -(!M<r r i>r ( 8 - n ) 

and 

V(<r r *T) =(i + 2)/r^- (8-12) 

Now, the operator K(ip) = 70(1 — a; A V)?/>7o commutes with the Hamiltonian 
defined by (8.3). Constructing eigenstates of this operator with eigenvalue k, 
leads to solutions of the form 

Mx k) = {^ ru{r) + a ^ v ^ I(Ts * = i + i (813) 
\o- r vrwW<T 3 +vr^w k = -(i+i), 

where k is a non-zero integer and u(r) and v(r) are complex functions of r (i.e. 
sums of a scalar and an ier 3 term). Substituting (8.13) into equation (8.9), and 
using the properties of the spherical monogenics, we arrive at the coupled radial 
equations 

- ( T ,,/2 )©=0 <-> 

where 

/ k/t l (E + m)-(2M/ r y/ 2 (4r)- 1 \ 

\i(E-m)-(2M/r) 1 / 2 (Ar)- 1 -«/r )' K ' 

u\ and U2 are the reduced functions defined by 

mi = ru U2 — irv, (8.16) 

and the primes denote differentiation with respect to r. (We continue to employ 
the abbreviation i for Icr^.) 

To analyse (8.14) we first rewrite it in the equivalent form 



( 1 (2M/r) 1 / 2 
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This makes it clear that the equations have regular singular points at the origin 
and horizon (r = 2M), as well as an irregular singular point at r = oo. To 
our knowledge, the special function theory required to deal with such equations 
has not been developed. Without it we either attempt a numerical solution, or 
look for power series with a limited radius of convergence. Here we consider the 
latter approach, and look for power-series solutions around the horizon. To this 
end we introduce the scries 



■ = V s Y,hv k , (8-18) 

fc=0 fe=0 



Ul=T) a kV > "2 



where 77 = r — 2M. The index s controls the radial dependence of ip at the 
horizon, so represents a physical quantity. To find the values that s can take, 
we substitute equation (8.18) into (8.17) and set 77 = 0. This results in the 
equation 

s fa \ _(l l\ ( k/(2M) i{E + m) - {MI)-^ ( V 

2M \boJ~ V 1 V V( E - m ) - (8A/)- 1 -k/(2M) ) \b 0/ 

(8.19) 

Rewriting this in terms of a single matrix and setting its determinant to zero 
yields the two indicial roots 

s = and s = -\ + MM E. (8.20) 

The s = solution is entirely sensible — the power series is analytic, and 
nothing peculiar happens at the horizon. The existence of this root agrees 
with our earlier observation that one can evolve the time-dependent equations 
without encountering any difficulties at the horizon. The second root is more 
problematic, as it leads to solutions that are ill-defined at the horizon. Before 
assessing the physical content of these roots, however, we must first check that 
they are gauge invariant. 

If, instead of working in the Newtonian gauge, we keep the gauge unspecified 
then, after separating out the angular dependence, the equations reduce to 



L r L t \ fuA _ f K/r-G/2 im-F/2\( Ul 
L t L r ) \u 2 ) ~ \-im - F/2 -«/r - G/2 ) \u 2 



.21) 



We can still assume that the i-dcpcndcncc is of the form cxp(-iEt), so that 
equation (8.21) becomes 

(91 ^f^ =B M (8 . 22) 
\92 9iJ \u' 2 J \u 2 ) 

where 



n _( K/r-G/2 + if 2 E i{m + hE)-F/2\ 

\-i(m- f 1 E)-F/2 -K/r-G/2 + if 2 EJ ■ {0 

Now, since gi 2 — g 2 2 = 1 — 2M/r holds in all gauges, we obtain 

( i - 2 ^)cn = (_ fl ; 2 ;fW:;j< ^ 
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and substituting in the power series (8.18) and setting 77 = produces the 
indicial equation 



det 



9i -92 
-92 9i 



B -I 



= 0, (8.25) 

r=2M 



where / is the identity matrix. Employing the result that 

9iG - g 2 F = \d r ( 9l 2 - g 2 2 ) = M/r 2 (8.26) 

we find that the solutions to the indicial equation are now 

s = Q and s = -±+4iME(g 1 f 2 -g 2 f 1 ). (8.27) 

But in Section 6.4 we established that 31/2 — .92/1 = +1 at the horizon for 
all solutions with a forward time direction. This demonstrates that the indices 
are indeed intrinsic, with the sign of the imaginary term for the singular root 
picking up information about the time direction implicit in the presence of the 
horizon. 

The fact that s = is always a solution of the indicial equation means that 
solutions always exist that arc analytic at the horizon. Determining the split 
between ingoing and outgoing states of these solutions enables one to calculate 
reflection coefficients and scattering amplitudes [34]. The question we wish to 
consider is whether the second, singular, root can be physically significant. To 
address this we look at the current. The covariant current J is given by ipjQijj, 
and satisfies T>-J = 0. The corresponding non-covariant quantity is therefore 

J = h("07oV>) det(h) -1 (8.28) 

which satisfies the fiatspace conservation equation V • J = 0. It is the stream- 
lines of J that are plotted as functions of x and determine the flow of density. 
The crucial terms in J are the time component and radial component, which 
(ignoring the overall exponential decay term) are given by 

1q-J = \(fi{uiui + u 2 u 2 ) + f 2 { Ul u 2 + u lU2 ))^r (8-29) 



and 

e r -J = -j(gi(intj 2 + uiu 2 ) + g 2 {u 1 u 1 + Uzuz))^™^™ ■ (8.30) 

The {fi} and {gi} are finite for all admissible solutions so, for the s = solution, 
the components of J are well-defined at the horizon. Furthermore, it is easily 
shown that for s = the radial flux at the horizon always points inwards. The 
s = root therefore describes the case where the flux crosses the horizon and 
continues onto the singularity. 

For the singular root we must first decide on a branch for the solution so that 
ip is fully specified on both sides of the horizon. We can then assess whether 
the discontinuity in -0, and the discontinuity in the current generated by it, arc 
physically acceptable. To do this, we first write 77 s as 

rf = exp((-| +4dME)1n(r-2M)). (8.31) 



93 



We can now write 



ln(r - 2M) = In |r - 2M\ + i arg(r - 2M ), (8.32) 
and for the choice of argument we take 

arg(r-2M) = (° r >^' (8.33) 
I — 7T r < 2M. 

(The choice of sign will be discussed further below.) If we now take the limit 
r — ► 2AI from above and below we find that the 70 component of J is given by 

l0 .J = A(e,^e-^\r-2Mr^^x\ 1 " >2 f (8.34) 

I exp (8irM E r ) r < 2M 

where A{6, 4>) is a positive-definite, finite term and we have split E into real and 
imaginary parts as 

E = E r - it. (8.35) 

Equation (8.34) is valid in all gauges for which g\ji — .92/1 = +1 at the horizon. 
While the density 70 -J is singular at the horizon, the presence of the positive 
term 8Me ensures that any integral over the horizon is finite and the solution is 
therefore normalisablc. This link between the properties of ip at the horizon and 
at the origin (where e is determined) provides another example of the importance 
of finding global solutions to the field equations. The radial current now turns 
out to be 

^•J=^#e^V-2Mr^xj 1 , , r>2M (8.36) 



4M 1 1 -exp(87rM£L) r < 2M 



and is therefore zero at the horizon, and inward-pointing everywhere inside the 
horizon. It appears that the existence of the imaginary contribution to E does 
ensure that that the singular solutions have sensible physical properties, and the 
singularity in i/> at the horizon is no worse than that encountered in the ground 
state of the hydrogen atom [89]. What is less clear, however, is the extent 
to which the properties of ?/> at the horizon are compatible with the original 
equation (8.1). In particular, since tp is both singular and non-differentiable at 
the horizon, it does not appear that the singular root can be viewed as defining 
a solution of (8.1) over all space. 



8.2 The Hawking temperature 

A number of authors have attempted to give derivations of the Hawking tem- 
perature and particle flux due to a black hole from an analysis of first-quantised 
theory, i.e. from the properties of wave equations alone [90, 91, 92]. This 
work has generated some controversy [93] , so it is interesting to assess how the 
ideas stand up in GTG. These derivations focus on the singular solutions to 
the wave equation (cither Klein-Gordon or Dirac), and study the properties of 
these solutions under the assumption that the energy is real. If one ignores the 
problems that e = introduces for the normalisability of ip and presses ahead, 
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then from (8.36) there is now a non-zero current at the horizon and, further- 
more, there is a net creation of flux there. The ratio of the outward flux to the 
total flux is simply 

e r -J+-t r -J- = e ME + 1 (8 ' 3?) 
which defines a Fermi-Dirac distribution with temperature 

T= — — . (8.38) 

8nMk B 

Remarkably, this is the temperature found by Hawking [94] . The fact that both 
the correct Fermi-Dirac statistics and Hawking temperature can be derived in 
this manner is astonishing, since both arc thought to be the result of quantum 
field theory. But what can we really make of this derivation? The first prob- 
lem is that setting e to zero means that the density is no longer normalisable 
at the horizon — any integral of the density over the horizon region diverges 
logarithmically, which is clearly unphysical. A further problem relates to the 
choice of branch (8.33). Had the opposite branch been chosen we would not have 
obtained (8.37) and, as pointed out in [93], there is no a priori justification for 
the choice made in (8.33). 

For the above reasons, the 'derivation' of (8.37) cannot be viewed as being 
sound. The remarkable thing is that the same techniques can be used to 'de- 
rive' the correct temperatures for the Reissner-Nordstrom and Rindler cases, as 
well as the Schwinger production rate in a constant electric field. This will be 
discussed elsewhere. These further analyses contain another surprise: the tem- 
perature at the interior horizon of a Reissner-Nordstrom black hole is necessarily 
negative! However, while these analyses are both interesting and suggestive, it is 
only through a study of the full quantum field theory in a black hole background 
that one can be sure about particle production rates. 

A final point in this section is that all our analyses have been based on work- 
ing with the correct time-asymmetric solutions admitted in GTG. On attempt- 
ing to force through the analysis in the 'Schwarzschild' gauge (172 = ji = 0), 
one discovers that the indices are now given by 

s = -\±AiME. (8.39) 

In this case no analytic solution is possible, and even the presence of an expo- 
nential damping factor does not produce a normalisable current at the horizon. 
This only serves to reinforce the importance of working with global solutions, 
since there is no doubt that the presence of non-singular, normalisable solutions 
is an intrinsic feature of horizons. 



8.3 The Dirac equation in a cosmological background 

As a second application we consider the Dirac equation in a cosmological back- 
ground. We have a choice of form of h-function to use, of which the simplest is 
that defined by (6.171), 

h'(a) = a-e t e t + -((1 - kr 2 ) l/2 a-e r e r + aAov a r ). (8.40) 
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The Dirac equation in this background takes the form 



e% + - kr 2 ) 1/2 e r d r + e 8 d e + e+d*)) yjla 



S 

+ (^H(t)e t - ^((1 - fcr 2 ) 1 / 2 - l)e r J ^Ia 3 = m^yo, (8.41) 

where the various functions are as defined in Section 6.6. Our question is this: 
can we find solutions to (8.41) such that the observables are homogeneous? 
There is clearly no difficulty if k = since, with p = 0, equation (8.41) is solved 

by 

and the observables are fixed vectors which scale as p(t) in magnitude [95]. 
But what happens when k ^ 0? It turns out that the solution (8.42) must be 
modified to [95] 

V; = P ~^rnt (8 43) 

1 + V 1 — fcr 2 

For the case of k > both h(a) and ip are only defined for r < fc -1 / 2 . This 
problem is overcome by using the displacement (7.56) to transform to the global 
solution of equation (7.57). In this case ip is given by 

V' = |(1 + kr 2 /A)p 1 l 2 e- I<T * mt , (8.44) 

which diverges as r — > oo. 

In both the k > and k < cases, ip contains additional r-depcndence and 
so is not homogeneous. Furthermore, the observables obtained from ip are also 
inhomogeneous. In principle one could therefore determine the origin of this 
space from measurements of the current density. This clearly violates the prin- 
ciple of homogeneity, though it is not necessarily inconsistent with experiment. 
The implications for cosmology of this fact are discussed in the following section. 
(Some consequences for self-consistent solutions of the Einstein-Dirac equations 
are discussed in [95, 96] and, in the context of general relativity, in [97].) 

The fact that quantum fields see this 'preferred' direction in fc ^ models, 
whereas classical phenomena do not, reflects the gauge structure of the theory. 
Dirac spinors are the only fields whose action couples them directly to the w(a)- 
function. All other matter fields couple to the gravitational field through the 
h-ficld only. Dirac spinors therefore probe the structure of the gravitational 
fields directly at level of the w-field, which is inhomogeneous for fc ^ models. 
On the other hand, classical fields only interact via the covariant quantities 
obtained from the gravitational fields, which are homogeneous for all values of 
fc. This conclusion is reinforced by the fact that the Klein-Gordon equation, 
for which the action does not contain the w(a)-field, does have homogenous 
solutions in a fc 7^ universe. 



9 Implications for cosmology 

In Section 6.6 we discussed some aspects of cosmology as examples of the gen- 
eral treatment of time-varying spherically-symmetric systems. There we drew 
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attention to the utility of the Newtonian gauge as a tool for tackling problems 
in cosmology. In addition, in Sections 7.4 and 8.3 we studied the Maxwell and 
Dirac equations in various cosmological backgrounds. In this section we draw 
together some of our conclusions from these sections. Specifically, we discuss 
redshifts, difficulties with k ^ models, and the definitions of mass and energy 
for cosmological models. 

9.1 Cosmological redshifts 

As a final demonstration of the use of the Newtonian gauge, consider a photon 
following a null path in the 9 = n/2 plane. In this case the photon's momentum 
can be written as 

P = <f>R( lo + j^R, (9.1) 

where 

R = e a/2Icr 3 ( 92 ) 

and $ is the frequency measured by observers comoving with the fluid. We 
restrict to the pressureless case, so G = and /i = 1, but will allow p to be 
r-dependent. A simple application of equation (4.63) produces 

d T $ = -$ 2 (— sin 2 X + d r g 2 cos 2 X ) , (9.3) 

where \ = 4> + a - But, since /i = 1, we find that d T t = <I>, so 

^ = -$(^sin 2 x + a i , 92 cos 2 X ), (9.4) 

which holds in any spherically-symmetric pressureless fluid. 

For the case of a cosmological background we have 52 = H(i)r, so the angular 
terms drop out of equation (9.4), and we are left with the simple equation 

^ = -H(t)$=f p $>, (9.5) 

which integrates to give the familiar redshift versus density relation 

l + z= ( Pl /p Q ) 1/3 . (9.6) 

Other standard cosmological relations, such as the luminosity distance and an- 
gular diameter versus rcdshift formulae, can be easily derived in this gauge (see 
also Section 6.6). 

In [95] equation (9.6) was derived in the gauge of equation (6.171), in which 
all particles comoving with the cosmological fluid are at rest in the background 
spacetime. In this gauge the redshift can be attributed to a loss of energy 
to the gravitational field, although this is a gauge-dependent viewpoint — the 
only physical statement that one can make is embodied in equation (9.6). The 
explanation of cosmological redshifts in our theory has nothing to do with 'tired 
light', or spacetime playing a dynamic role by expanding, or even anything to do 
with Doppler shifts. The redshift is simply a consequence of the assumption of 
homogeneity. Ultimately, all physical predictions are independent of the gauge 
in which they are made, although certain gauges may have useful computational 
or heuristic value. 
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9.2 k 7^ cosmologies 

At the level of classical (i.e. non-quantum) physics, there is no doubt that k 7^ 
cosmologies are homogeneous. This is true in both general relativity and GTG. 
No prediction derived for classical systems of point particles or for electromag- 
netic fields can reveal a preferred spatial direction in these models. However, 
we saw in Section 8.3 that it is impossible to find homogeneous solutions of the 
Dirac equation in a k 7^ universe. The consequences of this are, in princi- 
ple, observable, since the local density gradient will reveal a preferred radial 
direction. It has already been pointed out that it is impossible to find a self- 
consistent solution of the combined system of Einstcin-Dirac equations for any 
case other than a spatially flat cosmology [95, 97]. We believe that this is the 
first time that it has been pointed out that even a non-sclf-consistcnt Dirac field 
would be observably inhomogeneous. This is a more damaging result for k =/= 
cosmologies, since it reveals inhomogeneity without assuming that spin-torsion 
effects have anything to do with the dynamics of the universe. 

While the properties of Dirac fields pose theoretical difficulties for k 7^ 
models, there is no contradiction with present observations. Furthermore, one 
could question the validity of inferences drawn in extrapolating the Dirac equa- 
tion to cosmological scales. There is, however, a purely classical effect that does 
lead one to question the validity of k > models. As we have seen, when looking 
at the properties of fields in a k > background it is necessary to work in a 
globally-defined gauge, such as that of equation (7.57). In this case the Maxwell 
equations show that each point charge must have an image charge present in a 
remote region of the universe. This is a consequence of a closed universe that 
we have not seen discussed, although it has doubtless been pointed out before. 
The necessity for this image charge raises many problems in attempting to take 
such a universe seriously. 

9.3 Mass and energy in cosmological models 

Setting aside the problems with k 7^ models, a further issue on which our 
theory sheds some light is discussions of the total matter and energy content 
of the universe. In Section 6.5 we discussed the distinction between the rest- 
mass energy and the total gravitating energy inside a sphere of radius r. Since 
cosmological models are a special case of the general theory outlined in Section 6, 
this same distinction should be significant in cosmology. 

In Section 6.5 we identified the total gravitating energy of a sphere of radius 
r with the function M(r, t). For all cosmological models, this is given by (6.156) 

M(r,t) = p r 3 p(t). (9.7) 

In strictly homogeneous models the Weyl tensor vanishes, and we lose an intrin- 
sically-defined distance scale. But, if we consider cosmological models as the 
limiting case of spherically-symmetric systems, then there seems little doubt 
that (9.7) is still the correct expression for the gravitating energy within a 
sphere of radius r surrounding the origin. Moreover, attempts to discover the 
gravitating content of a region rely on perturbations away from ideal uniformity. 
In these cases an intrinsic distance scale is well defined, since a Weyl tensor is 
again present. (Determining the gravitating content of a region is important 
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in, for example, determinations of — the ratio of the actual density of the 
universe to the critical density.) On the other hand, the total rest mass energy 
within a sphere of radius r centred on the origin must still be given by (6.134) 



p(r, t)= 4irs 2 p(t) 



ds 



(9.8) 



o 



This remains a covariant scalar quantity, and is just p{i) multiplied by the 
covariant volume integral (this is the volume one would measure locally using 
light paths or fixed rods). 

We have now defined two covariant scalar quantities, M(r,t) and p(r,t), 
both of which are conserved along fluid streamlines in the absence of pressure. 
If the identifications made in the spherically-symmetric case remain valid in the 
homogeneous case, then the difference between these should give the additional 
contribution to the total energy beyond the rest-mass energy. For spatially-flat 
universes we have g\ = 1 so there is no difference. (In terms of the Newtonian- 
gauge description of Section 6.5, the gravitational potential energy cancels the 
kinetic energy.) But, for k ^ models, there is a difference because p{r,t) is 
now given by 



An interesting place to study the difference between M and p is in a k > 
universe at its 'turnaround' point, as described in Section 6.6. There one finds 
that, to lowest order in r, the difference is given by 



which is precisely the Newtonian formula for the self-potential of a uniform 
sphere of mass M{r). This is what we would have expected since, at the 
turnaround point, the kinetic energy vanishes. 

The above should only be viewed as suggestive, but one idea that it appears 
to rule out is the popular suggestion that the total energy density of the universe 
should be zero [98, 99]. If the above analysis is correct then there is no possibility 
of the total energy density M(r, t) ever being zero. Furthermore, for spatially- 
flat models — which we consider the most likely — the total energy density 
resides entirely in the rest-mass energies of the particles in the universe and 
cannot be cancelled by a negative gravitational contribution. 

10 Conclusions 

In this paper we developed a theory of gravity consisting of gauge fields defined 
in a flat background spacetime. The theory is conceptually simple, and the 
role of the gauge fields is clearly understood — they ensure invariance under 
arbitrary displacements and rotations. While it is possible to maintain a classical 
picture of the rotation gauge group, a full understanding of its role is only 
achieved once the Dirac action is considered. The result is a theory which offers a 
radically different interpretation of gravitational interactions from that provided 
by general relativity. Despite this, the two theories agree in their predictions 
for a wide range of phenomena. Differences only begin to emerge over issues 




(9.9) 



M(r) - p(r) 



3M(r) 2 
5r 



(9.10) 
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such as the role of topology, our insistence on the use of global solutions, and 
in the interaction with quantum theory. Furthermore, the separation of the 
gauge fields into one for displacements and one for local rotations is suggestive 
of physical effects being separated into an inertia field and a force field. Indeed, 
there is good reason to believe that mass should enter relativistic multiparticle 
wave equations in the manner of the h-field [27]. It is possible that, in the 
development of a multiparticle theory, the h and fl fields will be extended in 
quite distinct ways. Such possibilities do not appear to be open to a theory such 
as general relativity, with its reliance on the metric as the 'foundation of all' [73]. 
Probably the closest approach to the theory developed here is the spin-2 field 
theory discussed by many authors (see Box 17.2 of [73] and [50]). This theory 
is usually viewed as reproducing general relativity exactly, albeit in a somewhat 
ugly form due to the existence of a background spacetime and the reliance on 
infinite series of the field variable. By contrast, we hope to have demonstrated 
that GTG has an internal attractiveness of its own, as well as simplicity due to 
its first-order nature. 

A crucial question to address is whether any experimental tests are likely to 
distinguish between general relativity and GTG in the immediate future. The 
biggest differences between general relativity and GTG to emerge to date lie in 
the treatment of black hole singularities [55, 58], but these arc unlikely to be 
testable for some considerable time! A more promising area is the link between 
gravity and quantum spin. GTG makes a clear prediction for the type and 
magnitude of this interaction, whereas it is not uniquely picked out in general 
spin-torsion theories, or in more general Poincare gauge theory. Any experiment 
measuring this interaction would therefore provide a clear test of GTG. A partial 
exploration of the effects of spin interactions in GTG is contained in [49] . 

The techniques developed here reveal some remarkable properties of spher- 
ically-symmetric systems. It has been known since the 1930's [100] that, in the 
absence of pressure, the dynamical equations of cosmology can be cast in a New- 
tonian form. We have now shown that a single, unified, Newtonian treatment 
can be given for all spherically-symmetric pressurclcss fluids, whether homoge- 
neous or not. Furthermore, effects that have hitherto been viewed as the result 
of spacetime curvature can now be understood in a simple alternative fashion. 
The result is a physical picture in which the background spacetime has no effect 
on either dynamics or kinematics. This, we believe, is both new and poten- 
tially very useful. For example, simulations of black hole formation, and studies 
of the behaviour of the universe as a whole, can be carried out in exactly the 
same framework. All previous studies have relied on cutting and pasting various 
metrics together, with the result that no clear, global view of the underlying 
physics can be achieved. These advantages arc exploited in [77] to model the 
growth of a spherically-symmetric perturbation in a homogeneous background 
cosmology, and to study the effect of the perturbation on the cosmic microwave 
background. 

The intrinsic method described here, and used to study spherically-symmet- 
ric systems, is quite general and can be applied to a wide range of problems. 
In [61] the method is applied to a restricted class of cylindrically-symmetric 
systems, and elsewhere we have presented treatments of axisymmctric sys- 
tems [60, 34]. In all cases studied to date, the intrinsic method has brought 
considerable clarity to what would otherwise be a largely mysterious mess of al- 
gebra. This is achieved by removing the dependence on an arbitrary coordinate 
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system, and instead working directly with physical quantities. The same tech- 
nique also looks well suited to the study of cosmological perturbation theory, 
about which there has been considerable recent debate [101]. 

The interaction between Dirac theory and the gauge theory developed here 
revealed a number of surprises. The first was that consistency of the action 
principle with the minimal-coupling procedure restricted us to a theory that 
is unique up to the possible inclusion of a cosmological constant. The second 
was that spatially flat cosmologies are the only ones that are consistent with 
homogeneity at the level of the single-particle Dirac equation. The final surprise 
was provided by a study of the Dirac equation in a black hole background, 
which revealed a remarkable link with the Hawking temperature and quantum 
field theory. Much work remains to settle the issues raised by this final point, 
however. 

As a final remark, we also hope to have demonstrated the power of geometric 
algebra in analysing many physical problems. Many of the derivations performed 
in this paper would have been far more cumbersome in any other language, and 
none arc capable of the compact expressions provided by geometric algebra 
for, say, the Ricmann tensor. In addition, use of geometric algebra enabled us 
to remove all reference to coordinate frames from the fundamental equations. 
This is a real aid to providing a clear physical understanding of the mathematics 
involved. Wc would encourage anyone interested in studying the consequences 
of our theory to take time to master the techniques of geometric algebra. 
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A The Dirac operator algebra 

In the Dirac Pauli representation the 7-matriccs arc defined as 
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, % = 1 





(A.l) 



where the <Xfc are the standard Pauli matrices [89, 102]. The Dirac 7-matrices 
act on spinors, which are four-component complex column vectors. A spinor 
\tp) is placed in one-to-one correspondence with an 8-component even element 
of the spacetime algebra via [21, 103] 
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-b 3 + ib° 
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The action of the {7^}, i and 75 = - 
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a k Ia k +I(b° + b k I<r k ). (A.2) 



170717273 operators maps to 
<->■ 7^70 

<-> ipI(T 3 (A.3) 
<-> 7/10-3. 
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The Dirac equation. 

r{id ll -eA lt )\1>)=m\il>), (A.4) 
now takes the spacetime algebra form 

^{d^Ia-^ - eA M V)7o = mip. (A.5) 

Recombining to form the vectors V = 7 M 9 M and A = j^A^ , and postmultiplying 
by 7o, we arrive at the Dirac equation in the form 

VipI(Ts — eAip = w07o. (A. 6) 

Under Lorentz transformations the spinor ip transforms single-sidedly to Rij), 
hence the presence of the fixed 70 and 73 vectors on the right-hand side of tp 
does not break Lorentz invariance. 

The role of the Dirac adjoint is played by the geometric operation of re- 
version, and the quantum inner product projects out the {1,/<t 3 } components 
from a general multivcctor. So. for example, the real part of the inner product 
{ipihps) is given in the spacetime algebra by ("01 "02) and the imaginary part by 
— ('0iV , 2^o'3)- The Dirac current J M = (V>|7 M |V') is now replaced by the set of 
components 

(■fa^iho) =7 M -W>7oV0 ■ (A.7) 
These are simply the components of the vector ipjoip, decomposed in the {7^} 
frame. Reference to the frame is removed from the vector by defining the current 
as 

J = ^7o</>- (A.8) 
Similarly, the role of the spin current is played by the vector 

s = ip<y 3 i>, (A.9) 

and the spin trivector is simply Is. The Dirac Lagrangian has the equivalent 
form 

{ipK^iid^ - eA M ) - m)\ip) *-> (Vipl~f 3 ib - eAipwip - mipi>), (A.10) 

which is the form used in the main text. A more detailed discussion of the 
spacetime algebra formulation of Dirac theory is contained in [27]. 

B Some results in multivector calculus 

We begin with a set of results for the derivative with respect to the vector a in 
an n-dimensional space [11]: 



d a a-b = b d a a 2 = 2a 

d a -a = n d a a-A r = rA r 

<9 a Aa = d a ahA r = (n — r)A r 
d a a = n d a A r a = (— l) r (n — 2r)A r . 

The results needed for the multivector derivative in this paper are: 

d x (XA) =P X (A) 
d x (XA)=P x (A), 



(B.l) 



(B.2) 
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where Px{A) is the projection of A onto the grades contained in X. These 
results are combined using Leibniz' rule; for example, 

(Hj) = a (#) + a (v4> =2$. (b.3) 

For the action principle we also require results for the multivector derivative 
with respect to the directional derivatives of a field ip. The aim is again to 
refine the calculus so that it becomes possible to work in a frame-free manner. 
(The derivations presented here supersede those given previously in [31].) Wc 
first introduce a fixed orthonormal frame {e J }, with reciprocal {e^}, so that 
&*-&k — £fe- The partial derivative of tp with respect to the coordinate x J = eJ -x 
is abbreviated to ip t j so that 

tp d =e j -'Vip. (B.4) 
Wc can now define the frame-free derivative 



The operator d^, a is the multivector derivative with respect to the a-derivative 
of ip. The fundamental property of a is that 

fy,a(&-VV>M) = a-bP i ,{M). (B.6) 

Again, more complicated results are built up by applying Leibniz' rule. The 
Euler-Lagrange equations for the Lagrangian density £ = £(-0, a- V^) can now 
be given in the form 

d i ,L = d a -V{d^ a C), (B.7) 

which is the form applied in the main text of this paper. 

We also need a formalism for the derivative with respect to a linear function. 
Given the linear function h(a) and the fixed frame {e^}, we define the scalar 
coefficients 

/',, ',-h;«,;. (B.8) 

The individual partial derivatives . are assembled into a frame-free derivative 
by defining 

#h(a) = "•' .<>>, (B.9) 
The fundamental property of <9h(a) is that 

9h(o)h(6)-c = a-ejeid hi] {hi k b k c l ) 
= a-ej e^b 

= a-bc (B.10) 

which, together with Leibniz' rule, is sufficient to derive all the required prop- 
erties of the (?h( a ) operator. For example, if B is a fixed bivector, 

d h(a) (h(bAc)B) = 9 h (a)(h(6)h(c)S) - d h(a) (h(c)h(b)B) 
= a-bh(c)-B-a-ch(b)-B 

= h(a-(bAc))-B. (B.ll) 
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This result extends immediately to give 



d h{a) (h(A r )B r ) = (h(a-Ar)B r )i. (B.12) 

In particular, 

9 h(a) det(h)=^ (o) (h(/)/- 1 ) 
= h(a-/)/- 1 

= det(h)fr 1 (a)> (B.13) 

where the definition of the inverse (2.53) has been employed. The derivation 
of (B.13) affords a remarkably direct proof of the formula for the derivative of 
the determinant of a linear function. 

The above results hold equally if h is replaced throughout by its adjoint 
h, which is the form of the derivative used throughout the main text. Note, 
however, that 

<9h(a)h(6) = d h{a) {h(c)b)d c 
= a-cbd c 

= ba. (B.14) 

Thus the derivatives of h(6) and h(b) give different results, regardless of any 
symmetry properties of h. This has immediate implications for the symmetry 
(or lack of symmetry) of the functional stress-energy tensors for certain fields. 

We finally need some results for derivatives with respect to the bivector- 
valued linear function 0(a). The extensions are straightforward, and we just 
give the required results: 

d n{a) (n(b)M) =a-b(M) 2 (B.15) 
<9fi(&), (c-Vfi(d)M) =a-cb-d(M) 2 . (B.16) 



C The translation of tensor calculus 

The reformulation of the gauge theory presented in this paper in terms of con- 
ventional tensor calculus proceeds as follows. A choice of gauge is made and 
a set of scalar coordinates is introduced. We define the coordinate frame 

with the reciprocal frame denoted {e M }. In terms of this the vector derivative 
is 

V = e"-^-. (C.2) 
From the frame vectors we construct the vectors 

.o^h- 1 ^), ^ = h(e"). (C.3) 
The metric is then given by the 4x4 matrix 

g l iv = g IJ ,-gv (C.4) 
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If the x-dependence in g^ v is replaced by dependence solely on the coordinates 
{x^} then we recover Riemann-Cartan geometry, where all relations are between 
coordinates and the concept of a point as a vector is lost. 

The connection is defined by (following the conventions of [9, 44]) 

V»g v = T^,g a , (C.5) 

where X> M = g^-V = +ui(g fl )x . We can therefore write 

^ v = 9 X -(Pm9u). (C.6) 

Since 

dn9v\ = (p^g v )-g\+ gv^V^gx), (C.7) 

we find that 

d fl gux = T« v g a x + T« x g av , (C.8) 

which recovers 'metric compatibility' of the connection. This is nothing more 
than the statement that the operator satisfies Leibniz' rule. Equation (C.8) 
can be inverted to show that the connection contains a component given by the 
standard Christoffcl symbol. The connection can then be written 

where K x ^ is the contorsion tensor and is given by 

K\n — —Sx/j, + S\ v ^ — S v a m - (C.10) 

Here Sxp" is the torsion tensor, equal to the antisymmetric part of the connec- 
tion: 

*~>\p. — 2 \ L Aju 1 nX) 

= -±{ 9fl -(Vxgn-9A-(V^)) 
= -\{9^9\>{VAg v ) 

= l(gx/\g»)-S{g v ), (C.ll) 
where S(a) is the torsion bivector. The contorsion is formed from S(a) by 

Kx^ v = -\{gxf\g,i)-S{g lJ ) + l(9nAg v )-S(gx) - ^{g v /\gx)-S(g li ) 
= \{gn^9v)-S{g\) - \gx-{. 9l L-S{g v ) - g v -S{g ll )) 
= {g^gu)-{S{gx) - \gx-{d a hS{a))). (C.12) 

For our gauge theory of gravity, the torsion is exclusively of the type S(a) = a-S, 
where S is the spin trivector, in which case 

Kx^ v = -Sx^v = -\{gxf\9nf\g v )-S. (C.13) 

If we now consider the covariant derivative of a covariant vector A = A a g a = 
A a g a , we find that 

V^A = V^(A a g a ) 

= {d^A a )g a +A a Tl a g p 

= (d^A a +r« p A?)g a , (C.14) 
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so that the components of the vector T>^A are those expected for tensor calculus. 
Obviously the fact that A a g a = A a g a implies that = A a g all , so indices are 
raised and lowered in the expected manner. 

For covariant quantities such as the Ricmann tensor the translation to tensor 
calculus is straightforward: 

R^u P a = (9^g v )-K(g a Ag p ). (C.15) 

The general scheme is that any covariant quantity in GTG can be decomposed 
into tensor components by applying either the or {g^}, or a mixture of 

both, to yield a tensor with the appropriate number of upstairs and downstairs 
indices. So, for example, T can be decomposed to = T ■ {g^^gu), = 
T ■ (g^Ag") or F^ v = T ■ (g^ l\g v ). Tensor calculus is poor at revealing which, 
if any, of the components represent a physical observable. Such issues are much 
clearer in GTG, which focuses attention on the single entity T . 
A vierbein e^ 1 (essentially an orthonormal tetrad) is given by 

e/=^'7 4 (C16) 
e^ = g"- 7l (C.17) 

where {7*} is a fixed orthonormal frame. Any position dependence in the {7*} 
is eliminated with a suitable rotor transformation. When matrix operators {j 1 } 
are required these are replaced by the {7*} frame vectors using the method 
described in Appendix A. In this way frame-free vectors can be assembled. For 
example, the Dirac operator [45, Chapter 11] 

W) = c< f + 1 w J '7 fc ) H) (c.18) 

has the spacetimc algebra equivalent 

ff"-7<7 < + \u{g l i)-( lk A lj ) 1 J 1 k ^ V>7o = h(V)V>7o + Is'MflyOV'To 

= Dt/rfo. (C.19) 

The above relations enable many results from Ricmann-Cartan geometry to 
be carried over into our formalism, though the theory of gravity presented here 
is restricted in the structures from Riemann-Cartan geometry that it admits 
(the torsion is of trivector type, for example). 

A similar translation scheme is easily set up for the language of differential 
forms, which is much closer to the spirit of geometric algebra than tensor cal- 
culus. Differential forms are scalar- valued functions of an antisymmetrised set 
of vectors. They can easily be mapped to an equivalent multivector, and a full 
translation into geometric algebra is quite straightforward [11]. Here we note 
in passing the geometric algebra equivalent of the Hodge dual of a differential 
form, which is 

*a r h-> - dct(h)^ 1 hh(i r )i (C.20) 
where A r is the multivector equivalent of a r . 



106 



References 

[1] R. Utiyama. Invariant theoretical interpretation of interaction. Phys. 
Rev., 101(5):1597, 1956. 

[2] T.W.B. Kibble. Lorentz invariance and the gravitational field. J. Math. 
Phys., 2(3):212, 1961. 

[3] D. Ivanenko and G. Sardanashvily. The gauge treatment of gravity. Phys. 
Rep., 94(1):1, 1983. 

[4] E. Cartan. Sur un generalisation de la notion de courbure de Riemann et 
les espaces a torsion. C. R. Acad. Sci. (Paris), 174:593, 1922. 

[5] H. Weyl. A remark on the coupling of gravitation and electron. Phys. 
Rev., 77(5):699, 1950. 

[6] D. Sciama. The physical structure of general relativity. Rev. Mod. Phys., 
36:463 and 1103, 1964. 

[7] J. Weyssenhoff and A. Raabc. Rclativistic dynamics of spin-fluids and 
spin-particles. Acta Phys. Pol, 9:7, 1947. 

[8] O. Costa de Beauregard. Translational inertial spin effect. Phys. Rev., 
129(1):466, 1963. 

[9] F.W. Hehl, P. von der Heyde, G.D. Kerlick, and J.M. Nester. General 
relativity with spin and torsion: Foundations and prospects. Rev. Mod. 
Phys., 48:393, 1976. 

[10] T. Eguchi, P.B. Gilkey, and A.J. Hanson. Gravitation, gauge theories and 
differential geometry Phys. Rep., 66(6):213, 1980. 

[11] D. Hestenes and G. Sobczyk. Clifford Algebra to Geometric Calculus. 
Rcidcl, Dordrecht, 1984. 

[12] S.F. Gull, A.N. Lasenby, and C.J.L. Doran. Imaginary numbers arc not 
real — the geometric algebra of spacetime. Found. Phys., 23(9):1175, 
1993. 

[13] D. Hestenes. Curvature calculations with spacetime algebra. Int. J. Theor. 
Phys., 25(6):581, 1986. 

[14] D. Hestenes. Space-Time Algebra. Gordon and Breach, New York, 1966. 

[15] G. Sobczyk. Space-time approach to curvature. J. Math. Phys., 22(2):333, 
1981. 

[16] S.W. Hawking and G.F.R. Ellis. The Large Scale Structure of Space- Time. 
Cambridge University Press, 1973. 

[17] R. d'Inverno. Introducing Einstein's Relativity. Oxford University Press, 
1992. 

[18] W.J. Kaufmann. The Cosmic Frontiers of General Relativity. Penguin 
Books, 1979. 



107 



S.W. Hawking. Black Holes and Baby Universes and Other Essays. Ban- 
tam Books, London, 1993. 

D. Hestenes. Observables, operators, and complex numbers in the Dirac 
theory. J. Math. Phys., 16(3) :556, 1975. 

C. J.L. Doran, A.N. Lasenby, and S.F. Gull. States and operators in the 
spacetime algebra. Found. Phys., 23(9):1239, 1993. 

D. Hestenes. New Foundations for Classical Mechanics (Second Edition). 
Kluwcr Academic Publishers, Dordrecht, 1999. 

D. Hestenes. Proper particle mechanics. J. Math. Phys., 15(10):1768, 
1974. 

T.G. Void. An introduction to geometric algebra with an application to 
rigid body mechanics. Am. J. Phys., 61(6):491, 1993. 

D. Hestenes. Proper dynamics of a rigid point particle. J. Math. Phys., 
15(10):1778, 1974. 

S.F. Gull, A.N. Lasenby, and C.J.L. Doran. Electron paths, tunnelling and 
diffraction in the spacetime algebra. Found. Phys., 23(10):1329, 1993. 

C.J.L Doran, A.N. Lasenby, S.F. Gull, S.S. Somaroo, and A.D. Challinor. 
Spacetime algebra and electron physics. Adv. Imag. & Elect. Phys., 95:271, 
1996. 

T.G. Void. An introduction to geometric calculus and its application to 
electrodynamics. Am. J. Phys., 61(6):505, 1993. 

A.N. Lasenby, C.J.L. Doran, and S.F. Gull. Grassmann calculus, pseu- 
doclassical mechanics and geometric algebra. J. Math. Phys., 34(8):3683, 
1993. 

C.J.L. Doran, D. Hestenes, F. Sommen, and N. van Acker. Lie groups as 
spin groups. J. Math. Phys., 34(8):3642, 1993. 

A.N. Lasenby, C.J.L. Doran, and S.F. Gull. A multivector derivative 
approach to Lagrangian field theory. Found. Phys., 23(10):1295, 1993. 

C.J.L. Doran, A.N. Lasenby, and S.F. Gull. Grassmann mechanics, mul- 
tivector derivatives and geometric algebra. In Z. Ozicwicz, B. Jancewicz, 
and A. Borowiec, editors, Spinors, Twistors, Clifford Algebras and Quan- 
tum Deformations, page 215. Kluwer Academic, Dordrecht, 1993. 

A.N. Lasenby, C.J.L. Doran, and S.F. Gull. 2-spinors, twistors and su- 
persymmetry in the spacetime algebra. In Z. Ozicwicz, B. Jancewicz, and 
A. Borowiec, editors, Spinors, Twistors, Clifford Algebras and Quantum 
Deformations, page 233. Kluwer Academic, Dordrecht, 1993. 

C.J.L Doran and A.N. Lasenby. Geometric Algebra for Physicists. Cam- 
bridge University Press, 2003. 



108 



D. Hestenes. A unified language for mathematics and physics. In J.S.R. 
Chisholm and A.K. Common, editors, Clifford Algebras and their Appli- 
cations in Mathematical Physics (1985), page 1. Reidel, Dordrecht, 1986. 

D. Hestenes. The design of linear algebra and geometry. Acta Appl. Math., 
23:65, 1991. 

C.J.L. Doran and A.N. Lasenby. Physical Applications of Geometric Al- 
gebra. Cambridge Univcsity lecture course. Lecture notes available from 
http : //www.mrao . cam. ac . uk/^clif f ord. 

J.S.R. Chisholm and A.K. Common, eds. Clifford Algebras and their Ap- 
plications in Mathematical Physics (1985). Reidel, Dordrecht, 1986. 

A. Micali, R. Boudct and J. Helmstetter, eds. Clifford Algebras and their 
Applications in Mathematical Physics (1989). Kluwcr Academic, Dor- 
drecht, 1991. 

F. Brackx and R. Delanghe and H. Serras, eds. Clifford Algebras and 
their Applications in Mathematical Physics (1993). Kluwer Academic, 
Dordrecht, 1993. 

W.K. Clifford. Applications of Grassmann's extensive algebra. Am. J. 
Math., 1:350, 1878. 

H. Grassmann. Dcr ort der Hamilton'schcn quaternionen in der Aus- 
dehnungslehre. Math. Ann., 12:375, 1877. 

C.J.L. Doran, A.N. Lasenby, and S.F. Gull. Gravity as a gauge theory in 
the spacetime algebra. In F. Brackx, R. Delanghe, and H. Serras, editors, 
Clifford Algebras and their Applications in Mathematical Physics (1993), 
page 375. Kluwer Academic, Dordrecht, 1993. 

M. Nakahara. Geometry, Topology and Physics. Adam Hilger, Bristol, 
1990. 

M. Gockeler and T. Schuckcr. Differential Geometry, Gauge Theories, 
and Gravity. Cambridge University Press, 1987. 

R.T. Rauch. Equivalence of an R + R 2 theory of gravity to Einstein- 
Cartan-Sciama-Kibble theory in the presence of matter. Phys. Rev. D, 
26(4):931, 1982. 

R.D. Hccht, J. Lemke, and R.P. Wallncr. Can Poincare gauge theory be 
saved? Phys. Rev. D, 44(8):2442, 1991. 

A.V. Khodunov and V.V. Zhytnikov. Gravitational equations in space- 
time with torsion. J. Math. Phys., 33(10):3509, 1992. 

C.J.L Doran, A.N. Lasenby, A.D. Challinor, and S.F Gull. Effects of spin- 
torsion in gauge theory gravity. J. Math. Phys., 39(6):3303-3321, 1998. 

R.P. Fcynman, F.B. Morinigo, and W.G. Wagner. Feynman Lectures on 
Gravitation. Addison-Wesley, Reading MA, 1995. 



109 



[51] R. Penrose and W. Rindler. Spinors and space-time, Volume I: two-spinor 
calculus and relativistic fields. Cambridge University Press, 1984. 

[52] D. Kramer, H. Stephani, M. MacCallum, and E. Herlt. Exact Solutions 
of Einstein's Field Equations. Cambridge University Press, 1980. 

[53] C.J.L. Doran. Geometric Algebra and its Application to Mathematical 
Physics. PhD thesis, Cambridge University, 1994. 

[54] E.T. Newman and A.I. Janis. Note on the Kerr spinning-particle metric. 
J. Math. Phys., 6(4):915, 1965. 

[55] C.J.L. Doran and A.N. Lasenby. Integral equations, Kerr-Schild fields 
and gravitational sources gr-qc/0404081, 2004. 

[56] F.B. Estabrook and H.D. Wahlquist. Dyadic analysis of spacetime con- 
gruences. J. Math. Phys., 5:1629, 1965. 

[57] H.D. Wahlquist. The problem of exact interior solutions for rotating rigid 
bodies in general relativity. J. Math. Phys., 33(1):304, 1992. 

[58] A.N. Lasenby, C.J.L. Doran, Y. Dabrowski, and A.D. Challinor. Rotat- 
ing astrophysical systems and a gauge theory approach to gravity. In 
N. Sanchez and A. Zichichi, editors, Current Topics in Astrofundamental 
Physics, Erice 1996, page 380. World Scientific, Singapore, 1997. 

[59] A.N. Lasenby, C.J.L. Doran, and S.F. Gull. Astrophysical and cosmo- 
logical consequences of a gauge theory of gravity. In N. Sanchez and 
A. Zichichi, editors, Advances in Astrofundamental Physics, Erice 1994, 
page 359. World Scientific, Singapore, 1995. 

[60] C.J.L. Doran and A.N. Lasenby. New techniques for analysing axisym- 
mctric gravitational systems I. Vacuum fields. Class. Quant. Gravity,. 
20(6):1077-1101, 2003. 

[61] C.J.L. Doran, A.N. Lasenby, and S.F. Gull. Physics of rotating cylindrical 
strings. Phys. Rev. D, 54(10):6021, 1996. 

[62] R. Gautreau. Curvature coordinates in cosmology. Phys. Rev. D, 
29(2):186, 1984. 

[63] R. Gautreau and J.M. Cohen. Gravitational collapse in a single coordinate 
system. Am. J. Phys., 63(11):991, 1995. 

[64] K. Martel and E. Poisson. Regular coordinate systems for Schwarzschild 
and other spherical spacetimes. Am. J. Phys., 69(4):476, 2001. 

[65] M.W. Choptuik. Universality and scaling in gravitational collapse of a 
massless scalar field. Phys. Rev. Lett., 70(1):9, 1993. 

[66] A.M. Abrahams and C.R. Evans. Critical behaviour and scaling in vacuum 
axisymmetric gravitational collapse. Phys. Rev. Lett., 70(20):2980, 1993. 

[67] R. Gautreau. Light cones inside the Schwarzschild radius. Am. J. Phys., 
63(5):431, 1995. 



110 



[68] H. Stephani. General Relativity. Cambridge University Press, 1982. 

[69] G. Lcmaitrc. Spherical condensations in the expanding universe. Acad. 
ScL, Paris, Comptes Rend., 196:903, 1933. 

[70] M.D. Kruskal. Maximal extension of the Schwarzschild metric. Phys. 
Rev., 119:1743, 1960. 

[71] J.R. Oppenheimer and H. Snyder. On continued gravitational contraction. 
Phys. Rev., 56:455, 1939. 

[72] C.W. Misner and D.H. Sharp. Relativistic equations for adiabatic, spher- 
ically symmetric gravitational collapse. Phys. Rev., 136(2B):571, 1964. 

[73] C.W. Misner, K.S. Thorne, and J. A. Wheeler. Gravitation. W.H. Freeman 
and Company, San Francisco, 1973. 

[74] M. Panck. Cosmic background radiation anisotropics from cosmic struc- 
tures: Models based on the Tolman solution. Ap.J., 388:225, 1992. 

[75] R.C. Tolman. Effect of inhomogeneity on cosmological models. Proc. Nat. 
Acad. Sci. U.S., 20:169, 1934. 

[76] H. Bondi. Spherically symmetrical models in general relativity. Mon. Not. 
R. Astron. Soc, 107:410, 1947. 

[77] A.N. Lasenby, C.J.L Doran, M.P. Hobson, Y. Dabrowski, and A.D. Challi- 
nor. Microwave background anisotropics and nonlinear structures I. Im- 
proved theoretical models. Mon. Not. R. Astron. Soc, 302:748, 1999. 

[78] J. D. Bckcnstein. Generalized second law of thermodynamics in black-hole 
physics. Phys. Rev. D, 9(12):3292, 1974. 

[79] I.D. Novikov and V.P. Frolov. Physics of Black Holes. Kluwer Academic, 
Dordrecht, 1989. 

[80] J.V. Narlikar. Introduction to Cosmology. Cambridge University Press, 
1993. 

[81] G.F.R. Ellis and T. Rothman. Lost horizons. Am. J. Phys., 61(10):883, 
1993. 

[82] E.T. Copson. On electrostatics in a gravitational field. Proc. R. Soc. A, 
118:184, 1928. 

[83] B. Linet. Electrostatics and magnetostatics in the Schwarzschild metric. 
J. Phys. A, 9(7):1081, 1976. 

[84] R. S. Hanni and R. Rufiini. Lines of force of a point charge near a Schwarz- 
schild black hole. Phys. Rev. D., 8(10):3259, 1973. 

[85] K.S. Thorne, R.H. Price, and D.A. Macdonald. Black Holes: The Mem- 
brane Paradigm. Yale University Press, 1986. 

[86] A.G. Smith and CM. Will. Force on a static charge outside a Schwarz- 
schild black hole. Phys. Rev. D, 22(6):1276, 1980. 



Ill 



[87] C.J.L. Doran. New form of the Kerr solution. Phys. Rev. D, 61(6):067503, 
2000. 

[88] A.N. Lasenby et al. Bound states and decay times of fermions in a Schwarz- 
schild black hole background, gr-qc/0209090, 2002. 

[89] J.D. Bjorken and S.D. Droll. Relativistic Quantum Mechanics, vol 1. 
McGraw-Hill, New York, 1964. 

[90] T. Damour and R. Rufhni. Black-hole evaporation in the Klein-Sauter- 
Hciscnberg-Eulcr formalism. Phys. Rev. D, 14(2):332, 1976. 

[91] Zhao Zheng, Guci Yuan-xing, and Liu Liao. Hawking evaporation of Dirac 
particles. Chin. Phys., 1(4):934, 1981. 

[92] Zhao Zheng, Guei Yuan-xing, and Liu Liao. On the Hawking evaporation 
of Dirac particles in Kerr-Newman space-time. Chin. Phys., 2(2):386, 
1982. 

[93] M. Martellini and A. Treves. Comment on the Damour-Rufhni treatment 
of black-hole evaporation. Phys. Rev. D, 15(8):2415, 1977. 

[94] S. W. Hawking. Black hole explosion? Nature, 248:30, 1974. 

[95] A.N. Lasenby, C.J.L. Doran, and S.F. Gull. Cosmological consequences of 
a flat-space theory of gravity. In F. Brackx, R. Delanghe, and H. Serras, 
editors, Clifford Algebras and their Applications in Mathematical Physics 
(1993), page 387. Kluwer Academic, Dordrecht, 1993. 

[96] A.D. Challinor, A.N. Lasenby, C.J.L Doran, and S.F Gull. Massive, non- 
ghost solutions for the Dirac field coupled self-consistcntly to gravity. Gen- 
eral Rel. Grav., 29:1527, 1997. 

[97] C.J. Isham and J.E. Nelson. Quantization of a coupled Fermi field and 
Robertson- Walker metric. Phys. Rev. D, 10(10):3226, 1974. 

[98] E.P. Tryon. Is the Universe a vacuum fluctuation? Nature, 246:396, 1973. 

[99] E.P. Tryon. What made the world? New Sci., 8 March: 14, 1984. 

[100] E.A. Milne and W.H. McCrca. Newtonian Universes and the curvature of 
space. Q. J. Maths., 5:73, 1934. 

[101] G.F.R. Ellis. The covariant and gauge invariant approach to perturbations 
in cosmology. In N. Sanchez and A. Zichichi, editors, Current Topics in 
Astro fundamental Physics: The Early Universe, page 1. Kluwer Academic, 
Dordrecht, 1995. 

[102] C. Itzykson and J-B. Zuber. Quantum Field Theory. McGraw-Hill, New 
York, 1980. 

[103] S.F. Gull. Charged particles at potential steps. In A. Weingartshofer and 
D. Hestenes, editors, The Electron, page 37. Kluwer Academic, Dordrecht, 
1991. 



112 



